[Wc2007]剪刀石头布

http://www.lydsy.com/JudgeOnline/problem.php?id=2597

Time Limit: 20 Sec Memory Limit: 128 MBSec Special Judge

Description

在一些一对一游戏的比赛（如下棋、乒乓球和羽毛球的单打）中，我们经常会遇到A胜过B，B胜过C而C又胜过A的有趣情况，不妨形象的称之为剪刀石头布情况。有的时候，无聊的人们会津津乐道于统计有多少这样的剪刀石头布情况发生，即有多少对无序三元组(A, B, C)，满足其中的一个人在比赛中赢了另一个人，另一个人赢了第三个人而第三个人又胜过了第一个人。注意这里无序的意思是说三元组中元素的顺序并不重要，将(A, B, C)、(A, C, B)、(B, A, C)、(B, C, A)、(C, A, B)和(C, B, A)视为相同的情况。

有N个人参加一场这样的游戏的比赛，赛程规定任意两个人之间都要进行一场比赛：这样总共有场比赛。比赛已经进行了一部分，我们想知道在极端情况下，比赛结束后最多会发生多少剪刀石头布情况。即给出已经发生的比赛结果，而你可以任意安排剩下的比赛的结果，以得到尽量多的剪刀石头布情况。

Input

输入文件的第1行是一个整数N，表示参加比赛的人数。

之后是一个N行N列的数字矩阵：一共N行，每行N列，数字间用空格隔开。

在第(i+1)行的第j列的数字如果是1，则表示i在已经发生的比赛中赢了j；该数字若是0，则表示在已经发生的比赛中i败于j；该数字是2，表示i和j之间的比赛尚未发生。数字矩阵对角线上的数字，即第(i+1)行第i列的数字都是0，它们仅仅是占位符号，没有任何意义。

输入文件保证合法，不会发生矛盾，当i≠j时，第(i+1)行第j列和第(j+1)行第i列的两个数字要么都是2，要么一个是0一个是1。

Output

输出文件的第1行是一个整数，表示在你安排的比赛结果中，出现了多少剪刀石头布情况。

输出文件的第2行开始有一个和输入文件中格式相同的N行N列的数字矩阵。第(i+1)行第j个数字描述了i和j之间的比赛结果，1表示i赢了j，0表示i负于j，与输入矩阵不同的是，在这个矩阵中没有表示比赛尚未进行的数字2；对角线上的数字都是0。输出矩阵要保证合法，不能发生矛盾。

Sample Input

3
0 1 2
0 0 2
2 2 0

Sample Output

1
0 1 0
0 0 1
1 0 0

首先最多的合法三元组转化为最少的不合法三元组

每个三元组中，不合法的情况是胜利者有2条出边

所以每个人对不合法三元组的贡献是C（win[i],2）

win[i]表示i的胜利次数

ans=C（n，3）- Σ （win[i],2）

所以最小化 Σ （win[i],2）

考虑最小费用流

C（1,2）=0

C（2,2）=1

C（3,2）=3

C（4,2）=6

每次对答案的贡献一次+1、+2、+3、+4……

所以构图方式：

首先累加已知的Σ C（win[i],2）

源点向每一场比赛连流量为1、费用为0的边

这场比赛如果已知结果，向胜利者连流量为1、费用为0的边

如果不知道结果，向两人都连流量为1、费用为0的边

对于每个人，都向汇点连n-1条边，前win[i]条流量为1，费用为0

其余的流量为1，费用为win[i]、win[i]+1、win[i]+2……

#include<queue>
#include<cstdio>
#include<cstring>
#define N 5500
#define M 50000
using namespace std;
int win[101],cnt;
int src,decc;
int front[N],nxt[M],to[M],tot=1,from[M];
int cost[M],cap[M],pre[N];
int dis[N],ans[101][101];
bool v[N];
queue<int>q;
void add(int u,int v,int w,int val)
{
    to[++tot]=v; nxt[tot]=front[u]; front[u]=tot; cap[tot]=w; cost[tot]=val; from[tot]=u;
    to[++tot]=u; nxt[tot]=front[v]; front[v]=tot; cap[tot]=0; cost[tot]=-val; from[tot]=v;
}
bool spfa()
{
    memset(dis,127,sizeof(dis));
    memset(v,0,sizeof(v));
    dis[src]=0; q.push(src); v[src]=true;
    int now;
    while(!q.empty())
    {
        now=q.front();
        q.pop(); v[now]=false;
        for(int i=front[now];i;i=nxt[i])
        {
            if(cap[i]>0&&dis[to[i]]>dis[now]+cost[i])
            {
                dis[to[i]]=dis[now]+cost[i];
                pre[to[i]]=i;
                if(!v[to[i]]) 
                {
                    v[to[i]]=true;
                    q.push(to[i]);
                }
            }
        }
    }
    return dis[decc]<2e9;
}
int main()
{
    int n,x;
    scanf("%d",&n);
    cnt=n+1;decc=n+1;
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      {
          scanf("%d",&x);
          if(j>=i) continue;
        if(!x) 
        {
            win[j]++;
            ans[j][i]=1;
            cnt++;
            add(src,cnt,1,0);
            add(cnt,j,1,0);
        }
          else if(x==1) 
        {
              win[i]++;
              ans[i][j]=1;
              cnt++;
              add(src,cnt,1,0);
              add(cnt,i,1,0);
        }
          else 
          {
              cnt++;
              add(src,cnt,1,0);
              add(cnt,i,1,0);
              add(cnt,j,1,0);
        }
      }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=win[i];j++) add(i,decc,1,0);
        for(int j=win[i]+1;j<n;j++) add(i,decc,1,j-1);
    }
    int tmp=0;
    while(spfa()) 
    {
        tmp+=dis[decc];
        for(int i=pre[decc];i;i=pre[from[i]])
        {
            cap[i]--; cap[i^1]++;
        }
    }
    for(int i=front[src];i;i=nxt[i])
     {
         int j1=front[to[i]],j2=nxt[j1];
        if(cap[j1]) ans[to[j2]][to[j1]]=1;
        else ans[to[j1]][to[j2]]=1;
     }
    int sum=n*(n-1)*(n-2)/6;
    for(int i=1;i<=n;i++) 
     sum-=win[i]*(win[i]-1)/2;
    printf("%d
",sum-tmp);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<n;j++)
         printf("%d ",ans[i][j]);
        printf("
");
    }
}

错误：一边spfa一边更新答案

原因：费用流也有退流