题意
对于一个N个点M条有向边的简单图,问我们最多还能加多少条边进去 使得该图仍然是有向图且非强联通。必须保证加边以后的图不是强联通图
解题思路:
首先缩点,因为环可以不必再考虑。对于两个联通分量,加边使得他们不连通则必须是一个块的所有点指向另一个块的所有点。一个N个点的图最多有N*(N-1)条边,已经连接了M条边,那么还剩下N*(N-1)-M条边 。剩下的点中有部分边是不可取的,因为要保证整个图是非强联通图。那么对于缩点后的图,对于某个度为0的点,他必定不能有任何的连出去/ 其他点连进来 ,我们只需要一个这样子的点即可。那么该点意味着该点内的点数P 不能与其他N-P个点有边。即上述N*(N-1)-M边中的不可取边。
所以只需要最小化P*(N-P)那么我们的答案就是最大化的
#include<cstring>
#include<cstdio>
#include<queue>
#define int long long
#define Re int
const int N=1e5+300,M=3e5+3;
int n,m,x,y,Q_o,ans,A[N],ru[N],gs[N],ip[N],dp[N],cu[N];
inline void in(Re &x){
int f=0;x=0;char c=getchar();
while(c<'0'||c>'9')f|=c=='-',c=getchar();
while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c^48),c=getchar();
x=f?-x:x;
}
inline int min(Re a,Re b){return a<b?a:b;}
inline int max(Re a,Re b){return a>b?a:b;}
struct Tarjan{
int o,t,dfn_o,Q[N],low[N],dfn[N],pan[N],head[N];
struct QAQ{int x,to,next;}a[M];
inline void add(Re x,Re y){a[++o].x=x,a[o].to=y,a[o].next=head[x],head[x]=o;}
inline void cle(){
o=1;
t=0;
dfn_o=0;
memset(head,0,sizeof(head));
memset(low,0,sizeof(low));
memset(dfn,0,sizeof(dfn));
Q_o=0;
}
inline void tarjan(Re x){
dfn[x]=low[x]=++dfn_o,Q[++t]=x,pan[x]=1;
for(Re i=head[x],to;i;i=a[i].next)
if(!dfn[to=a[i].to])tarjan(to),low[x]=min(low[x],low[to]);
else if(pan[to])low[x]=min(low[x],dfn[to]);
if(low[x]==dfn[x]){
++Q_o;
while(1){
ip[Q[t]]=Q_o,gs[Q_o]+=A[Q[t]],pan[Q[t]]=0;
if(x==Q[t--])break;
}
}
}
inline void SuoPoint(){for(Re i=1;i<=n;++i)if(!dfn[i])tarjan(i);}
}T1;
struct Tuopu{
int o,pan[N],head[N];std::queue<int>Q;
struct QAQ{int to,next;}a[M<<1];
inline void add(Re x,Re y){a[++o].to=y,a[o].next=head[x],head[x]=o;}
inline void creat(){
for(Re i=1;i<=T1.o;++i)
if((x=ip[T1.a[i].x])!=(y=ip[T1.a[i].to]))
++ru[y],add(x,y),cu[x]++;
}
inline void tuopu(){
for(Re i=1;i<=Q_o;++i)if(!ru[i])dp[i]=gs[i],Q.push(i);
while(!Q.empty()){
Re x=Q.front();Q.pop();
for(Re i=head[x],to;i;i=a[i].next){
to=a[i].to;
dp[to]=max(dp[to],dp[x]+gs[to]);
if(!(--ru[to]))Q.push(to);
}
}
}
}T2;
signed main(){
int ca=0;
int cab;
scanf("%lld",&cab);
while(cab--)
{
// in(n),in(m);
scanf("%lld%lld",&n,&m);
int edg=m;
ca++;
T1.cle();
memset(cu,0,sizeof(cu));
memset(ru,0,sizeof(ru));
memset(T2.head,0,sizeof(T2.head));
memset(gs,0,sizeof(gs));
for(Re i=1; i<=n; ++i)
A[i]=1;
while(m--)
in(x),in(y),T1.add(x,y);
T1.SuoPoint(),T2.creat();//T2.tuopu();
int ma=1e16;
for(int i=1; i<=Q_o; i++)
{
if(!ru[i]||!cu[i])
ma=min(ma,((n-gs[i])*gs[i]));
}
printf("Case %lld: ",ca);
if(Q_o==1)
{
printf("-1\n");
continue;
}
printf("%lld\n",(n*(n-1)-edg-ma));
}
}