BZOJ 2400: Spoj 839 Optimal Marks|最小割

bzoj评测姬已经被群攻了..
可以发现每一位都是独立的
然后枚举每一位,每一个点要么是1要么是0
然后可以划分为两个集合,0集合和1集合,只有连接两个集合的边才会对答案有贡献,然后跑最小割就可以了
建图:
枚举每一位

  1. S连这一位上是1的点,流量inf
  2. 这一位上是0的点连T,流量inf
  3. 再加上原图中的边,流量为1

这样跑出来的最小割就是这一位对答案1的贡献
然后再看属于S集的点的个数就是对答案2的贡献

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
#define mx 1e9
#define lowbit(x) (x&(-x))
#define N 505
using namespace std;
int sc()
{
    int i=0,f=1; char c=getchar();
    while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i*f;
}
long long sum,Sum;
int dis[505],q[505],mark[505];
int head[505],nxt[9444],lst[9444],c[9444];
int v[505],a[2002],b[2002];
int n,m,S,T,tot;
void insert(int x,int y,int a,int b)
{
    lst[++tot]=y; nxt[tot]=head[x]; head[x]=tot; c[tot]=a;
    lst[++tot]=x; nxt[tot]=head[y]; head[y]=tot; c[tot]=b;
}
void cal(int x,int y,int k)
{
    if(v[y]<0)
        insert(x,y,1,1);
    else 
    {
        if(v[y]&k) insert(y,x,1,0);
        else       insert(x,y,1,0);
    }
}
bool BFS()
{
    for(int i=1;i<=T;i++)dis[i]=0; 
    dis[q[1]=S]=1;
    for(int l=1,r=2;l<r;l++)
    {
        int x=q[l];
        for(int i=head[x];i;i=nxt[i])
            if(c[i]&&dis[lst[i]]==0)
            {
                dis[lst[i]]=dis[x]+1;
                q[r++]=lst[i];
            }
    }
    return dis[T];
}
int dfs(int x,int f)
{
    if(x==T)return f;
    int ww=0,w;
    for(int i=head[x];i;i=nxt[i])
        if(c[i]&&dis[lst[i]]==dis[x]+1)
        {
            w=dfs(lst[i],min(f-ww,c[i]));
            c[i]-=w;c[i^1]+=w; ww+=w;
            if(ww==f)return f;
        }
    if(ww==0)dis[x]=0;
    return ww;
}
int bfs()
{
    memset(mark,0,sizeof(mark));
    int r=2; mark[q[1]=S]=1;
    for(int l=1;l<r;l++)
    {
        int x=q[l];
        for(int i=head[x];i;i=nxt[i])
            if(c[i]&&!mark[lst[i]])
                mark[q[r++]=lst[i]]=1;
    }
    return r-2;
}           
int main()
{
    n=sc(),m=sc();
    for(int i=1;i<=n;i++)v[i]=sc();
    for(int i=1;i<=m;i++)a[i]=sc(),b[i]=sc();
    S=n+1,T=S+1;
    for(int i=1<<30;i;i>>=1)
    {
        memset(head,0,sizeof(head));
        tot=1;
        for(int j=1;j<=m;j++)
        {
            if(v[a[j]]<0||v[b[j]]<0)
            {
                if(v[a[j]]<0)
                    cal(a[j],b[j],i);
                else
                    cal(b[j],a[j],i);
            }
            else if((v[a[j]]&i)&&!(v[b[j]]&i)) insert(a[j],b[j],1,0);
            else if((v[b[j]]&i)&&!(v[a[j]]&i)) insert(b[j],a[j],1,0);
        }

        for(int j=1;j<=n;j++)
            if(v[j]>=0)
            {
                if(v[j]&i)insert(S,j,mx,0);
                else      insert(j,T,mx,0);
            }
        int ans=0;
        while(BFS())ans+=dfs(S,mx);
        sum+=(long long)ans*i;
        Sum+=(long long)i*(bfs());
    }
    cout<<sum<<endl<<Sum;
    return 0;
}

你可能感兴趣的:(最小割)