【bzoj3698】XWW的难题 有源汇上下界网络流最大流

Description

XWW是个影响力很大的人，他有很多的追随者。这些追随者都想要加入XWW教成为XWW的教徒。但是这并不容易，需要通过XWW的考核。
XWW给你出了这么一个难题：XWW给你一个N*N的正实数矩阵A，满足XWW性。
称一个N*N的矩阵满足XWW性当且仅当：（1）A[N][N]=0；（2）矩阵中每行的最后一个元素等于该行前N-1个数的和；（3）矩阵中每列的最后一个元素等于该列前N-1个数的和。
现在你要给A中的数进行取整操作（可以是上取整或者下取整），使得最后的A矩阵仍然满足XWW性。同时XWW还要求A中的元素之和尽量大。

Input

第一行一个整数N,N ≤ 100。
接下来N行每行包含N个绝对值小于等于1000的实数，最多一位小数。

Output

输出一行，即取整后A矩阵的元素之和的最大值。无解输出No。

Sample Input

4

3.1 6.8 7.3 17.2

9.6 2.4 0.7 12.7

3.6 1.2 6.5 11.3

16.3 10.4 14.5 0

Sample Output

129

HINT

【数据规模与约定】

有10组数据，n的大小分别为10,20,30…100。

【样例说明】

样例中取整后满足XWW性的和最大的矩阵为：

3 7 8 18

10 3 0 13

4 1 7 12

17 11 15 0

题解
有上下界的网络流。
学习文章：http://blog.csdn.net/clove_unique/article/details/54884437

代码

#include<cstdio>
#include<iostream>
#include<queue>
#define ll long long
#define inf 1000000000
#define mod 1000000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
queue<int>q;
int S,T,SS,TT;
int n,cnt=1,tot,ans;
int dis[205],in[205],cur[205];
double a[205][205];
int Head[205],Next[100005],ret[100005],len[100005];
void ins(int u,int v,int w)
{
    ret[++cnt]=v;Next[cnt]=Head[u];Head[u]=cnt;len[cnt]=w;
}
void insert(int u,int v,int w)
{ins(u,v,w);ins(v,u,0);}
void build()
{
    for (int i=1;i<n;i++)
    {
        if (a[i][n]!=(int)a[i][n])
            insert(S,i,1);
        in[S]-=(int)a[i][n];in[i]+=(int)a[i][n];
    }
    for (int i=1;i<n;i++)
    {
        if (a[n][i]!=(int)a[n][i])
            insert(i+n,T,1);
        in[i+n]-=(int)a[n][i];in[T]+=(int)a[n][i];
    }
    for (int i=1;i<n;i++)
        for (int j=1;j<n;j++)
        {
            if (a[i][j]!=(int)a[i][j])
                insert(i,j+n,1);
            in[i]-=(int)a[i][j];in[j+n]+=(int)a[i][j];
        }
    for (int i=1;i<=TT;i++)
        if (in[i]>0){tot+=in[i];insert(SS,i,in[i]);}
        else insert(i,TT,-in[i]);
}
bool bfs(int x,int y)
{
    for (int i=1;i<=TT;i++) dis[i]=-1;
    dis[x]=0;q.push(x);
    while (!q.empty())
    {
        int now=q.front();q.pop();
        for (int i=Head[now];i;i=Next[i])
        {
            if (len[i]&&dis[ret[i]]==-1)
            {
                dis[ret[i]]=dis[now]+1;
                q.push(ret[i]);
            }
        }
    }
    return dis[y]!=-1;
}
int dfs(int x,int y,int aug)
{
    if (x==y||aug==0) return aug;
    int flow=0;
    for (int i=cur[x];i;i=Next[i])if (dis[ret[i]]==dis[x]+1)
    {
        int t=dfs(ret[i],y,min(aug-flow,len[i]));flow+=t;
        len[i]-=t;len[i^1]+=t;if (len[i]) cur[x]=i;
        if (flow==aug) return flow;
    }
    if (flow==0) dis[x]=-1;
    return flow;
}
void dinic(int x,int y)
{
    while(bfs(x,y))
    {
        for(int i=1;i<=TT;i++)cur[i]=Head[i];
        ans+=dfs(x,y,inf);
    }
}
int main()
{
    n=read();S=2*n+1;T=S+1;SS=T+1;TT=SS+1;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            scanf("%lf",&a[i][j]);
    build();
    insert(T,S,inf);
    dinic(SS,TT);
    if (tot!=ans)return puts("No"),0;
    ans=0;
    dinic(S,T);
    printf("%d",ans*3);
    return 0;
}