[北京省选集训2019]图的难题 题解

题面

首先你会发现一个性质：一个图符合条件当且仅当E<=2*V-2;

因为假设图中有2*V-2条边，为了符合条件，一定是有一颗白色生成树和一颗黑色生成树。可以看出，没有其他的情况；

那么如果再多一条边呢？显然，填在任意位置都会使任意一颗树变成基环树，导致一种颜色成环，不成立；

所以E的最大值就是2*V-2；

所以我们就是判断任意一个子图的E-2*V<=-2;

为了方便，我们只要判定一个E-2*V最大的子图是否<=-2就好了；

怎么求呢？这时候引入最大权闭合子图的概念就可以秒出答案了(可以看我的这篇文章)

 把边化成点，点权为1，原来图中的点的点权是-2
 
 因为选边就一定要选择点，所以所有选择的点的后继一定要都被选择；
 
 然后就构成了最大权闭合子图的模型；
 
 注意一个细节，需要枚举所有点表示强制选这个点，否则答案是不成立的；

然后就AC掉了；

 

#include <bits/stdc++.h>
#define inc(i,a,b) for(register int i=a;i<=b;i++)
using namespace std;
class littlestar{
    public:        
        int to;
        int w;
        int nxt;
        void add(int u,int v,int gg);
}star[100010];
int head[100010],cnt=1;
inline void littlestar::add(int u,int v,int gg)
{
    to=v;
    nxt=head[u];
    head[u]=cnt;
    w=gg;
}
template<class nT>
inline void read(nT &x)
{
    char c;
    while(c=getchar(),!isdigit(c));
    x=c^48;
    while(c=getchar(),isdigit(c)) x=x*10+c-48;
}
int S,T;
int dis[100010];
int n,m;
bool bfs()
{
    memset(dis,0,sizeof(dis));
    queue<int> q;
    q.push(S);
    dis[S]=1;
    while(q.size()){
        int u=q.front();
        q.pop();
        for(int i=head[u];i;i=star[i].nxt){
            int v=star[i].to;
            if(!dis[v]&&star[i].w){
                dis[v]=dis[u]+1;
                q.push(v);
                if(v==T) return 1;
            }
        }
    }
    return 0;
}
int dinic(int u,int flow)
{
    if(u==T) return flow;
    int rest=flow,tmp;
    for(int i=head[u];i&&rest;i=star[i].nxt){
        int v=star[i].to;
        if(dis[v]==dis[u]+1&&star[i].w){
            tmp=dinic(v,min(star[i].w,rest));
            if(!tmp) dis[v]=0;
            star[i].w-=tmp;
            star[i^1].w+=tmp;
            rest-=tmp;
        }        
    }
    return flow-rest;
}
int sum=0;
int tmpu[10000],tmpv[10000];
int myenum(int s)
{
    memset(head,0,sizeof(head));
    cnt=1;
    inc(i,1,m){
        int u,v;
        u=tmpu[i]; v=tmpv[i];
        star[++cnt].add(S,i,1);
        star[++cnt].add(i,S,0);
        star[++cnt].add(i,m+u,INT_MAX);
        star[++cnt].add(m+u,i,0);
        star[++cnt].add(i,m+v,INT_MAX);
        star[++cnt].add(m+v,i,0);
    }
    inc(i,1,n){
        if(s==i) continue;
        star[++cnt].add(i+m,T,2);
        star[++cnt].add(T,i+m,0);
    }
    sum=m-2;
    while(bfs()){
        sum-=dinic(S,INT_MAX);
    }
    return sum;
}
int main()
{
    int tt;
    read(tt);
    while(tt--){
        read(n); read(m);
        S=n+m+1,T=n+m+2;
        int ans=-INT_MAX;
        inc(i,1,m) read(tmpu[i]),read(tmpv[i]);
        for(int i=1;i<=n;i++){
            ans=max(ans,myenum(i));
        }
        if(ans+2>0){
            cout<<"No"<<endl;
        }
        else{
            cout<<"Yes"<<endl;
        }
    }
}