您的位置:首页 > 其它

hdu 4635(tarjan)

2017-08-07 22:42 330 查看
题意:给出一个有向图,如果是联通直接输出-1,不是的话求最多可以增加多少条边使得图还是不联通的。

思路:先跑一遍他tarjan进行缩点,然后选取一个联通块作为与其他的联通块不联通的那个块,然后其他的联通块都联通上,这个能作为选取的那个联通块的话只有入度或者只有出度,因为当你把其他的点该连上的都连上的话就会使得这个联通块与其他的联通块相连,破坏连通性。

最大的点的话就是将两块都分别变成完全图,然后在他们的中间变成与入度或者出度同向的边的个数。是a * b,a,b分别为两个联通块的点个数,之后直接减去原来的边数就行了。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
typedef set<int>::iterator ITER;
#define fi first
#define se second
#define INF 0x3f3f3f3f
#define clr(x,y) memset(x,y,sizeof x)

const int maxn = 1e5 + 10;
const int Mod = 1e9 + 7;
const int N = 2;
int n,m;
struct Edge{int to,next;}edge[maxn * 2];
int edge_num;
int head[maxn];
int low[maxn],dfn[maxn];
int col[maxn],sizes[maxn];
int num[maxn];
int top;
int st[maxn];
int lens;
bool vis[maxn];
int in[maxn],out[maxn];
void Init()
{
edge_num = 0;clr(head,-1);clr(col,0);clr(dfn,0);clr(low,0);
clr(sizes,0);clr(vis,false);lens = 0;clr(num,0);clr(in,0);clr(out,0);
}
void add_edge(int x,int y){edge[edge_num] = (Edge){y,head[x]};head[x] = edge_num ++;}
void tarjan(int u,int &len)
{
dfn[u] = low[u] = len;
vis[u] = true;st[++ top] = u;
for(int i = head[u];i != -1; i = edge[i].next)
{
int v = edge[i].to;
if(!dfn[v])
{
tarjan(v,++ len);low[u] = min(low[u],low[v]);
}
else if(vis[v])low[u] = min(low[u],dfn[v]);
}
if(low[u] == dfn[u])
{
++ lens;vis[u] = false;col[u] = lens;sizes[lens] = 1;
// cout << lens << ": " << u <<" ";
while(st[top --] != u && top >= 1)
{
// cout << st[top + 1] << " ";
vis[st[top + 1]] = false;col[st[top + 1]] = lens;sizes[lens] ++;
}
// cout << endl;
}
}
int main()
{
int Tcase;
scanf("%d",&Tcase);
for(int ii = 1; ii <= Tcase; ii ++)
{
Init();
scanf("%d%d",&n,&m);
for(int i = 1; i <= m; i ++)
{
int x,y;scanf("%d%d",&x,&y);add_edge(x,y);
}
int top = 0;
for(int i = 1; i <= n; i ++)
{
if(!dfn[i])
{
int len = 1;tarjan(i,len);
}
}
if(lens == 1)
{
printf("Case %d: -1\n",ii);continue;
}
for(int u = 1; u <= n; u ++)
{
for(int i = head[u];i != -1; i = edge[i].next)
{
int v = edge[i].to;int x = col[u],y = col[v];
if(x != y)in[y] ++,out[x] ++;
}
}
// for(int i = 1; i <= lens; i ++)
// cout << sizes[i] << " ";
// cout << endl;
ll ans = 0;
for(int i = 1; i <= lens; i ++)
{
if(in[i] == 0 || out[i] == 0)
{
ll t = sizes[i] * (ll)(n - sizes[i]) +
sizes[i] * (ll)(sizes[i] - 1) + (n - sizes[i]) * (ll)(n - sizes[i] - 1) - m;
// cout << sizes[i] << " " << t <<endl;
ans = max(ans,t);
}
}
printf("Case %d: %I64d\n",ii,ans);
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: