2017Astar资格赛1002 度度熊的王国战略

啊啊全局最小割？

看了一天没看懂

考虑一个最小割，把原图分成两半（废话

左边取一个点为S，右边去一个点为T，跑出来的最小割一定是答案（废话

于是取1为S,随机T，1000组，保证不重复，随机数据随便过

但是有卡点，3000个点被分成了2999+1或2998+2

这两种情况都是可以预处理的

当2997+3时有三个点随机到一个就行，1000次，A掉的概率有,大概每个点70%

想想卡点不多（2999+1，2998+2为主），随机至上（笑），反正过了

#include <ctime>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define rep(j,k,l) for (int j=k;j<=l;++j)
#define red(j,k,l) for (int j=k;j>=l;--j)
#define N 3005
#define M 100005*2
#define inf 1000000007
#define min(a,b) (((a)<(b))?(a):(b))

using namespace std;
int n,m,a[M],b[M],c[M],d[M],e[M],pos;
struct Dinic{
    
    int st[N],to[M],ne[M],mx[M],fl[M],deep[N],que[N],cur[N],S,T,cnt;
	
    void add(int k,int l,int poi){
        
        to[++cnt]=l;
        mx[cnt]=poi;
        ne[cnt]=st[k];
        st[k]=cnt;
        fl[cnt]=0;
        //if (poi>0) cout<<k<<' '<<l<<' '<<poi<<endl;
        
    }
    
    void set(){
    	
    	S=1,T=d[++pos];
        cnt=1;rep(i,1,n) st[i]=0;
        rep(i,1,m) add(a[i],b[i],c[i]),add(b[i],a[i],c[i]);
        
    }
    
    bool bfs(){
        int head=0,tail=1;
        que[1]=S;memset(deep,0,sizeof(deep));
        deep[S]=1;
        while (head!=tail){
            int x=que[++head];
            for (int i=st[x];i;i=ne[i])
                if (deep[to[i]]==0&&fl[i]!=mx[i]){
                    
                    deep[to[i]]=deep[x]+1;
                    que[++tail]=to[i];
                    if (to[i]==T) return 1;
                    
                }
            
        }
        return deep[T]>0;
        
    }
    
    int dfs(int k,int re){
        
        if (k==T||re==0) return re;
        int flow=0,f;
        for (int i=cur[k];i;cur[k]=i=ne[i])
            if (deep[to[i]]==deep[k]+1&&(f=dfs(to[i],min(re,mx[i]-fl[i])))>0){
                
				//printf("--!--%d %d\n",k,f);
                fl[i]+=f;fl[i^1]-=f;
                flow+=f;re-=f;
                if (!re) return flow;
                
            }
        if (re) deep[k]=-1;
        return flow;
        
    }
    
    int maxflow(){
        
        int flow=0;
        while (bfs()){
            
			//puts("ok");
            rep(i,1,n) cur[i]=st[i];  //зЂвт
            flow+=dfs(S,inf);
            
        }
		//printf("----%d\n",flow);
        return flow;
        
    }
    
} F;


int main(){
	
	srand(time(0));
	while (scanf("%d%d",&n,&m)!=EOF){
		
		for (int i=1;i<n;i++) d[i]=i+1;
		random_shuffle(d+1,d+n);
		pos=0;
		memset(e,0,sizeof(e));
		rep(i,1,m)
			scanf("%d%d%d",a+i,b+i,c+i),e[a[i]]+=c[i],e[b[i]]+=c[i];
		int ans=inf,pi;
		rep(i,1,n) ans=min(ans,e[i]);
		rep(i,1,m) ans=min(ans,e[a[i]]+e[b[i]]-2*c[i]);
		rep(i,1,min(500,n-1)) F.set(),pi=F.maxflow(),ans=min(ans,pi);
		printf("%d\n",ans);
		
	}
    
}