BZOJ 2400(網絡流最小割)

思路：
首先比較顯然能夠想到二進制拆分
每一位的貢獻是獨立的
然後可以用最小割解決
S->i 表示這一個位置是0 inf
i-T 表示爲1 inf
然後所有的邊保留 1
這樣只要最後再找一遍t集合裏面的點就可以確定那些不確定的點了
第二問好像沒什麼用? 確定了每個位置後他的值也就確定了
10^9的限制是幹嘛的?

代碼

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
#include<queue>
#include<vector>
using namespace std;
typedef long long LL;
const int imax=500+9;
const int dmax=imax;
const int bmax=100000+229;
const int inf=100000229;
int n,m;
int a[imax],A[bmax],B[bmax],val[dmax],ind[dmax];
int S,T,num,head[dmax],inext[bmax],to[bmax],re[bmax];

void iread()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    for(int i=1;i<=m;i++) scanf("%d%d",&A[i],&B[i]);    
}

void findset(int x)
{
    ind[x]=1;
    for(int i=head[x];i!=-1;i=inext[i]) 
        if(!ind[to[i]] && re[i^1]) findset(to[i]); 
        //反向邊有流量說明這條邊有貢獻 
}

void iadd(int u,int v,int flow){ to[num]=v; re[num]=flow; inext[num]=head[u]; head[u]=num++;}
void add(int u,int v,int flow) { iadd(u,v,flow); iadd(v,u,0); }

int d[dmax]; 
queue<int> q; 
bool BFS()
{
    memset(d,0,sizeof(d));
    d[S]=1; q.push(S);
    while(!q.empty())
    {
        int u=q.front(); q.pop();
        for(int i=head[u];i!=-1;i=inext[i])
            if(re[i] && !d[to[i]])
            {
                d[to[i]]=d[u]+1;
                q.push(to[i]);  
            }
    }
    return d[T]!=0; 
}

int DFS(int x,int c)
{
    if(x==T || c==0) return c;
    int r=c;
    for(int i=head[x];i!=-1;i=inext[i])
        if(d[to[i]]==d[x]+1)
        {
            int ff=DFS(to[i],min(re[i],r));
            r-=ff; re[i]-=ff; re[i^1]+=ff;
            if(!r) break;   
        }   
    if(r==c) d[x]=0;
    return c-r;
}

int dinic()
{
    int nowans=0;
    while(BFS()) nowans+=DFS(S,inf);
//  printf("%d\n",nowans);
    return nowans;  
} 

void build(int k)
{
    S=0; T=n+1; num=0;
    for(int i=S;i<=T;i++) head[i]=-1;
    for(int i=1;i<=n;i++)
        if(a[i]>=0)
        {
            if(a[i]&(1<<(k-1))) add(i,T,inf);
            else add(S,i,inf);
        }
    for(int i=1;i<=m;i++)
    {
        add(A[i],B[i],1);
        add(B[i],A[i],1);
    }
}

void iwork()
{
    int Max=30;
    LL ans1=0; LL ans2=0;
    for(int i=1;i<=Max;i++) 
    {
        build(i);
        ans1+=(1LL<<(i-1))*dinic();

        memset(ind,false,sizeof(ind));
        findset(T);
        for(int j=1;j<=n;j++) 
            if(ind[j]) val[j]+=(1<<(i-1));  
    }
    for(int i=1;i<=n;i++)
    {
        if(a[i]>0) ans2+=a[i];
        else ans2+=val[i];
    }
    cout<<ans1<<endl<<ans2<<endl;
}

int main()
{
    iread();
    iwork();
    return 0;
}