BZOJ 3698: XWW的难题 [有源汇上下界最大流]
2017-03-28 21:44
381 查看
3698: XWW的难题
题意:(1)A=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 , s, t, tot, x, c, a ; struct edge{int v, c, f, ne, lower;}e[M]; int cnt=1, h ; 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 , head, tail, vis , d , cur ; 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] =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 [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); } }
相关文章推荐
- bzoj3698 XWW的难题(有源汇有上下界最大流)
- 【bzoj3698】XWW的难题 有源汇上下界网络流最大流
- bzoj3698: XWW的难题 //有上下界有源汇最大流
- bzoj 3698: XWW的难题 (有源汇有上下界的最大流)
- [BZOJ3698]XWW的难题(有源汇上下界最大流+讲解)
- [BZOJ3698]XWW的难题(有源汇有上下界的最大流)
- [BZOJ3698]XWW的难题解题报告|上下界网络流|有源汇最大流
- bzoj 3698 XWW的难题(有源汇的上下界最大流)
- [bzoj3698]XWW的难题 有源汇的上下界最大流
- BZOJ 3698: XWW的难题|有上下界的网络流之最大流
- 【BZOJ3698】XWW的难题 有上下界的最大流
- 3698: XWW的难题[有源汇上下界最大流]
- 3698: XWW的难题 有源汇上下界最大流
- 【bzoj3698】XWW的难题 有上下界最大流
- 【bzoj3698】【XWW的难题】【有上下界的网络流】
- [上下界网络流 二分] BZOJ 3698 XWW的难题
- 【bzoj3698】XWW的难题 有上下界的网络流
- BZOJ.3698.XWW的难题(有源汇上下界最大流ISAP)
- [BZOJ3698] XWW的难题 网络流
- BZOJ 3698 XWW的难题