BZOJ2597: [Wc2007]剪刀石头布（费用流）

传送门

题解：
从反方向考虑找不满足的个数$cnt$，那么答案是$\begin{pmatrix}{3\\n} \end{pmatrix}-cnt$。

从每个不满足的环考虑，因为原图是竞赛图，所以肯定有一个点连出了两条边，那么有:

$$cnt=\sum_i{\begin{pmatrix}{2\\d[i]}\end{pmatrix}}$$
$d[i]$为点 $i$出度。

转化可得

$$\begin{aligned}cnt&=&\sum_i\frac{(d[i]^2-d[i])}{2}\\cnt*2&=&-\begin{pmatrix}2\\n\end{pmatrix}+\sum_id[i]^2\end{aligned}$$

现在要最小化

$$\sum_id[i]^2$$
每个边向两点连边，每个点按照等差数列向汇点连边做费用流即可。

#include<bits/stdc++.h>
using namespace std;
inline int rd(){
    char ch=getchar();int i=0,f=1;
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){i=(i<<1)+(i<<3)+ch-'0';ch=getchar();}
    return i*f;
}
const int Maxn=1e4+50,Maxm=1e6+50,INF=0x3f3f3f3f;
int n,G[110][110],ans,d[110],tot;
int g[Maxn],nt[Maxm],c[Maxm],v[Maxm],dis[Maxn],w[Maxm],walk[Maxn],cur[Maxn],vt,ec=1;
const int src=0,des=10020;
inline void add(int x,int y,int o,int wi){
    ++ec;nt[ec]=g[x];g[x]=ec;v[ec]=y;w[ec]=wi;c[ec]=o;
    ++ec;nt[ec]=g[y];g[y]=ec;v[ec]=x;w[ec]=-wi;c[ec]=0;
}
inline bool spfa(){
    queue<int>q;
    for(int i=1;i<=tot;i++)dis[i]=INF;
    dis[des]=INF;dis[src]=0;q.push(src);
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int j=g[u];j;j=nt[j]){
            if(!c[j]||dis[v[j]]<=dis[u]+w[j])continue;
            dis[v[j]]=dis[u]+w[j];q.push(v[j]);
        }
    }
    return dis[des]<INF;
}
inline int dinic(const int &x,int f,int co){
    if(x==des)return ans+=1ll*f*co,f;
    int rs=0;walk[x]=vt;
    for(int &j=cur[x];j;j=nt[j]){
        if(!c[j]||(walk[v[j]]==vt)||(dis[v[j]]<dis[x]+w[j]||(dis[x]<=-INF)))continue;
        int o=dinic(v[j],min(f-rs,c[j]),co+w[j]);
        rs+=o;c[j]-=o;c[j^1]+=o;
        if(rs==f)return rs;
    }
    return dis[x]=-INF,rs;
}
inline int mincostflow(){
    while(spfa()){
        ++vt,memcpy(cur,g,sizeof(cur));
        while(dinic(src,INF,0))++vt,memcpy(cur,g,sizeof(cur));
    }
    return ans;
}
int main(){
    n=rd();tot=n;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++){
            G[i][j]=rd();if(i==j)continue;
            if(G[i][j]==1){
                d[i]++;ans+=d[i]*2-1;
            }else if(G[i][j]==2&&i<j){
                ++tot;add(src,tot,1,0);
                add(tot,i,1,0);add(tot,j,1,0);
            }
        }
    for(int i=1;i<=n;i++){
        for(int j=d[i]+1;j<n;j++){
            add(i,des,1,j*2-1);
        }
    }
    mincostflow();printf("%d\n",n*(n-1)*(n-2)/6-(ans-n*(n-1)/2)/2);
    for(int i=n+1;i<=tot;i++){
        int a=v[g[i]],b=v[nt[g[i]]];
        if(c[g[i]])swap(a,b);
        G[a][b]=1,G[b][a]=0;
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            printf("%d ",G[i][j]);
        }
        printf("\n");
    }
}