CCF-高速公路

问题描述

　　某国有n个城市，为了使得城市间的交通更便利，该国国王打算在城市之间修一些高速公路，由于经费限制，国王打算第一阶段先在部分城市之间修一些单向的高速公路。
　　现在，大臣们帮国王拟了一个修高速公路的计划。看了计划后，国王发现，有些城市之间可以通过高速公路直接（不经过其他城市）或间接（经过一个或多个其他城市）到达，而有的却不能。如果城市A可以通过高速公路到达城市B，而且城市B也可以通过高速公路到达城市A，则这两个城市被称为便利城市对。
　　国王想知道，在大臣们给他的计划中，有多少个便利城市对。

输入格式

　　输入的第一行包含两个整数n, m，分别表示城市和单向高速公路的数量。
　　接下来m行，每行两个整数a, b，表示城市a有一条单向的高速公路连向城市b。

输出格式

　　输出一行，包含一个整数，表示便利城市对的数量。

样例输入

5 5
1 2
2 3
3 4
4 2
3 5

样例输出

3

样例说明

　　城市间的连接如图所示。有3个便利城市对，它们分别是(2, 3), (2, 4), (3, 4)，请注意(2, 3)和(3, 2)看成同一个便利城市对。

评测用例规模与约定

　　前30%的评测用例满足1 ≤ n ≤ 100, 1 ≤ m ≤ 1000；
　　前60%的评测用例满足1 ≤ n ≤ 1000, 1 ≤ m ≤ 10000；
　　所有评测用例满足1 ≤ n ≤ 10000, 1 ≤ m ≤ 100000。

强联通分量的题目，对于每一个强联通分量求出点的个数num对应对数就是(num)*(num+1)/2。

敲了强联通分量的板子，1A，（装个B，嘿嘿嘿）

#include<iostream>#include<cmath>#include<cstring>#include<cstdio>#include<vector>#include<stack>#include<algorithm>#include<sstream>#define inf 0x3f3f3f3fusing namespace std;const int MAXN=10010;struct Edge{    int to,nextx;}edge[100010];stack<int> s;vector<int> g[MAXN];int dfn[MAXN];int low[MAXN];int is_stack[MAXN];int head[MAXN];int belong[MAXN];int numx[MAXN];int n,m,tp,p,num;void tarjan(int u){    dfn[u]=low[u]=++tp;    s.push(u);    is_stack[u]=1;    for(int i=head[u];i!=-1;i=edge[i].nextx)    {        int v=edge[i].to;        if(!dfn[v])        {            tarjan(v);            low[u]=min(low[u],low[v]);        }        else if(is_stack[v])            low[u]=min(low[u],dfn[v]);    }    if(dfn[u]==low[u])    {        num++;        while(!s.empty())        {            int v=s.top();            s.pop();            is_stack[v]=0;            belong[v]=num;            if(dfn[v]==low[u])                break;        }    }}int main(){    cin>>n>>m;    memset(dfn,0,sizeof(dfn));    memset(low,0,sizeof(low));    memset(head,-1,sizeof(head));    memset(is_stack,0,sizeof(is_stack));    memset(numx,0,sizeof(numx));    tp=0,p=0,num=0;    for(int i=1;i<=m;i++)    {        int u,v;        cin>>u>>v;        edge[i].to=v;        edge[i].nextx=head[u];        head[u]=i;    }    for(int i=1;i<=n;i++)        if(!dfn[i])        tarjan(i);    int sum=0;    for(int i=1;i<=n;i++)        numx[belong[i]]++;    for(int i=1;i<=num;i++)        if(numx[i]!=1)        sum+=(numx[i]-1)*(numx[i])/2;    cout<<sum<<endl;    return 0;}