BZOJ 3698: XWW的难题 [有源汇上下界最大流]

3698: XWW的难题

题意：（1）A[N][N]=0；（2）矩阵中每行的最后一个元素等于该行前N-1个数的和；（3）矩阵中每列的最后一个元素等于该列前N-1个数的和。给A中的数进行取整操作（可以是上取整或者下取整），使得最后的A矩阵仍然满足。要求A中的元素之和尽量大。

---

我不知花了多少个小时想为什么这样求出的最大流就是原图的最大流，并没有人来解释这件事...

建图太简单了，注意给出的可能本来就是整数
我把每种做法都试了一遍，都可以通过2333

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define fir first
#define sec second
typedef long long ll;
const int N=2005, M=4e5+5, INF=1e9;
inline ll read(){
    char c=getchar();ll x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}

int n, m, extra[N], s, t, tot, x, c, a[N][N];
struct edge{int v, c, f, ne, lower;}e[M];
int cnt=1, h[N];
inline int ins(int u, int v, int c, int b=0) {
	e[++cnt]=(edge){v, c, 0, h[u], b}; h[u]=cnt;
	e[++cnt]=(edge){u, 0, 0, h[v], b}; h[v]=cnt;
	return cnt-1;
}
int q[N], head, tail, vis[N], d[N], cur[N];
bool bfs(int s, int t) {
	memset(vis, 0, sizeof(vis));
	head=tail=1;
	q[tail++]=s; d[s]=0; vis[s]=1;
	while(head!=tail) {
		int u=q[head++];
		for(int i=h[u];i;i=e[i].ne)
			if(!vis[e[i].v] && e[i].c>e[i].f) {
				vis[e[i].v]=1; d[e[i].v]=d[u]+1;
				q[tail++]=e[i].v;
				if(e[i].v == t) return true;
			}
	}
	return false;
}
int dfs(int u, int a, int t) {
	if(u==t || a==0) return a;
	int flow=0, f;
	for(int &i=cur[u];i;i=e[i].ne)
		if(d[e[i].v]==d[u]+1 && (f=dfs(e[i].v, min(a, e[i].c-e[i].f), t))>0) {
			flow+=f;
			e[i].f+=f;
			e[i^1].f-=f;
			a-=f;
			if(a==0) break;
		}
	if(a) d[u]=-1;
	return flow;
}
int dinic(int s, int t) {
	int flow=0;
	while(bfs(s, t)) {
		for(int i=0; i<=tot; i++) cur[i]=h[i];
		flow+=dfs(s, INF, t);
	}
	return flow;
}

double f;
int main() {
	freopen("in","r",stdin);
	n=read(); s=0; t=n+n+1;
	for(int i=1; i<n; i++) {
		for(int j=1; j<n; j++) {
			scanf("%lf",&f); x=int(f); c= x!=f;
			a[i][j]=ins(i, n+j, c, x), extra[i]-=x, extra[n+j]+=x;
		}
		scanf("%lf",&f); x=int(f); c= x!=f;
		a[i][n]=ins(s, i, c, x), extra[s]-=x, extra[i]+=x;
	}
	for(int i=1; i<n; i++) {
		scanf("%lf",&f); x=int(f); c= x!=f;
		a[n][i]=ins(n+i, t, c, x), extra[n+i]-=x, extra[t]+=x;
	}
	int ss=n+n+2, tt=n+n+3, sum=0; tot=tt;
	for(int i=s; i<=t; i++) {
		if(extra[i]>0) ins(ss, i,  extra[i]), sum+=extra[i];
		if(extra[i]<0) ins(i, tt, -extra[i]);
	}
	ins(t, s, INF);
	int flow=dinic(ss, tt);
	if(flow!=sum) puts("No");
	else {
		//int flow=dinic(s, t); int  ans=0;
		//for(int i=1; i<=n; i++)
			//for(int j=1; j<=n; j++) ans += e[a[i][j]].f + e[a[i][j]].lower;
		int _ = e[cnt-1].f; e[cnt].f =e[cnt].c =e[cnt-1].f = e[cnt-1].c=0;
		printf("%d\n", (_ +dinic(s, t))*3);
	}
}