BZOJ 2597: [Wc2007]剪刀石头布|费用流

神题吓死人- -！费用流神题！！！
n个点考虑最多有C（n，3）个题意中描述的3元环
考虑每个人赢的次数
赢一次对这种3圆环的破坏个数为0，两次为1，三次再加2，四次再加3……
ans=C（n,3）-破坏个数
显然要让破坏个数最少，那就是费用流的模型了
原图中的每条边抽象为点
S连每条边流量为1费用0
每条边若胜负不定（2） 连两个人流量为1费用为0
若胜负确定（0,1）连胜利的人流量为1费用为0
然后每个人连到汇点n-1条边（其实这样显然多了，可以在前面计算一下连边的数目）流量都为1，费用依次为0,1,2,3,……

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<queue>
#include<cmath>
#include<vector>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
#define MX 1e9
#define T 111
#define E 11111
using namespace std;
int sc()
{
    int i=0; char c=getchar();
    while(c>'9'||c<'0')c=getchar();
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i;
}
int dis[E],vis[E],q[E],pre[E],from[E];
int head[E],lst[E*100],nxt[E*100],v[E*100],c[E*100];
int s[T][T],d[T];
int n,m,tot=1,num,ans;
int S,W;
void insert(int x,int y,int _c,int _v)
{
    lst[++tot]=y; nxt[tot]=head[x]; c[tot]=_c; v[tot]=_v; head[x]=tot;
    lst[++tot]=x; nxt[tot]=head[y]; c[tot]=0; v[tot]=-_v; head[y]=tot;
}
bool spfa()
{
    for(int i=1;i<=W;i++) dis[i]=MX; dis[S]=0;
    int l=1,r=2; q[1]=S;
    while(l!=r)
    {
        int x=q[l++]; l%=E; vis[x]=0;
        for(int i=head[x];i;i=nxt[i])
            if(c[i]&&dis[lst[i]]>dis[x]+v[i])
            {
                dis[lst[i]]=dis[x]+v[i];
                from[lst[i]]=x;
                pre[lst[i]]=i;
                if(!vis[lst[i]])
                {
                    vis[lst[i]]=1;
                    q[r++]=lst[i];
                    r%=E;
                }
            }
    }
    return dis[W]!=MX;
}
void mcf()
{
    int cc=MX;
    for(int i=W;i!=S;i=from[i])cc=min(cc,c[pre[i]]);
    for(int i=W;i!=S;i=from[i])
    {
        int k=pre[i];
        ans+=v[k]*cc;
        c[k]-=cc;
        c[k^1]+=cc;
    }
}
int find(int x)
{
    for(int i=head[x];i;i=nxt[i])
        if(!c[i])return lst[i];
}
int main()
{
    n=sc();m=n*(n-1)/2;
    S=n+m+1,W=S+1;num=n;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        {
            s[i][j]=sc();
            if(i>=j)continue;
            num++; insert(S,num,1,0);
            if(s[i][j]==1)insert(num,i,1,0),d[i]++;
            else if(s[i][j]==0)insert(num,j,1,0),d[j]++;
            else insert(num,i,1,0),insert(num,j,1,0),d[i]++,d[j]++;
        }
    for(int i=1;i<=n;i++)
        for(int j=0;j<d[i];j++)
            insert(i,W,1,j);
    while(spfa())mcf();
    cout<<n*(n-1)*(n-2)/6-ans<<endl;
    num=n;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            if(i==j){printf("0 ");continue;}
            if(i>j){printf("%d ",!s[j][i]);continue;}
            num++; int x=find(num);
            if(x==i)s[i][j]=1;else s[i][j]=0;
            printf("%d ",s[i][j]);
        }
        puts("");
    }       
    return 0;
}