bzoj [Usaco2005 mar]Ombrophobic Bovines 发抖的牛 最大流
2017-06-05 09:50
489 查看
很好的题目= =并没有想出正解= =
首先二分最短的路径x
然后每次二分的时候跑最大流。
每个牛拆成i和i’,起点向i连容量为数量的边,i’向终点连容量为容量的边。
如果i,j最短路小于二分的答案x就连i和j’容量为inf。
首先二分最短的路径x
然后每次二分的时候跑最大流。
每个牛拆成i和i’,起点向i连容量为数量的边,i’向终点连容量为容量的边。
如果i,j最短路小于二分的答案x就连i和j’容量为inf。
#include<cstdio> #include<cstring> #include<algorithm> #define fo(i,a,b) for(int i=a;i<=b;i++) #define fd(i,a,b) for(int i=a;i>=b;i--) #define inf 1000000000 using namespace std; const int N=1e5+5; int n,m; typedef long long ll; int a ,b ; int head ,next ,go ; int val ; int q ; ll f[1005][1005]; int dis ; int s,T,ans,tot,sum; inline void add(int x,int y,int z) { go[++tot]=y; next[tot]=head[x]; val[tot]=z; head[x]=tot; } inline void ins(int x,int y,int z) { add(x,y,z); add(y,x,0); } inline bool bfs() { int t=0,w=1; memset(dis,-1,sizeof(dis)); q[1]=s,dis[s]=0; while (t<w) { int x=q[++t]; for(int i=head[x];i;i=next[i]) { int v=go[i]; if (val[i]&&dis[v]==-1) { dis[v]=dis[x]+1; q[++w]=v; } } } return dis[T]!=-1; } inline int dfs(int x,int f) { if (x==T)return f; int w,used=0; for(int i=head[x];i;i=next[i]) { int v=go[i]; if (dis[v]==dis[x]+1) { w=f-used; w=dfs(v,min(w,val[i])); val[i]-=w; val[i^1]+=w; used+=w; if (used==f)return f; } } if (!used)dis[x]=-1; return used; } inline void dinic() { while (bfs())ans+=dfs(s,inf); } inline bool check(ll x) { s=0,T=2*n+1; tot=1;ans=0; memset(head,0,sizeof(head)); fo(i,1,n)ins(s,i,a[i]),ins(n+i,T,b[i]); fo(i,1,n) fo(j,1,n) if (f[i][j]<=x)ins(i,j+n,inf); dinic(); if (ans==sum)return 1; else return 0; } int main() { scanf("%d%d",&n,&m); fo(i,1,n) { scanf("%d%d",&a[i],&b[i]); sum+=a[i]; } fo(i,1,n) fo(j,1,n)f[i][j]=inf*200ll; fo(i,1,m) { int x,y; ll z; scanf("%d%d%lld",&x,&y,&z); f[x][y]=f[y][x]=min(f[x][y],z); } fo(i,1,n)f[i][i]=0; fo(k,1,n) fo(i,1,n) fo(j,1,n) f[i][j]=min(f[i][j],f[i][k]+f[k][j]); ll l=0,r=inf*200ll-1; ll ans=-1; while (l<=r) { ll mid=(l+r)>>1; if (check(mid)) { ans=mid; r=mid-1; } else l=mid+1; } printf("%lld\n",ans); }
相关文章推荐
- BZOJ 1738: [Usaco2005 mar]Ombrophobic Bovines 发抖的牛( floyd + 二分答案 + 最大流 )
- 【bzoj1738】[Usaco2005 mar]Ombrophobic Bovines 发抖的牛 二分答案+最大流
- 【BZOJ1738】[Usaco2005 mar]Ombrophobic Bovines 发抖的牛【二分】【Floyd】【最大流】
- BZOJ 1738: [Usaco2005 mar]Ombrophobic Bovines 发抖的牛 二分 网络流
- BZOJ1738: [Usaco2005 mar]Ombrophobic Bovines 发抖的牛
- 【bzoj1738】[Usaco2005 mar]Ombrophobic Bovines 发抖的牛 Floyd+二分+网络流最大流
- 【BZOJ1738】【POJ2391】【Usaco2005 mar】 Ombrophobic Bovines 发抖的牛 网络流
- BZOJ 1738: [Usaco2005 mar]Ombrophobic Bovines 发抖的牛
- BZOJ1738 [Usaco2005 mar]Ombrophobic Bovines 发抖的牛
- bzoj 1675: [Usaco2005 Feb]Rigging the Bovine Election 竞选划区(暴力)
- BZOJ1681 [Usaco2005 Mar]Checking an Alibi 不在场的证明
- bzoj1681[Usaco2005 Mar]Checking an Alibi 不在场的证明
- bzoj1681[Usaco2005 Mar]Checking an Alibi 不在场的证明
- BZOJ 1734: [Usaco2005 feb]Aggressive cows 愤怒的牛( 二分答案 )
- bzoj1682 [Usaco2005 Mar]Out of Hay 干草危机(kruskal裸题)
- bzoj 1741: [Usaco2005 nov]Asteroids 穿越小行星群【最大点覆盖】
- bzoj:1675 [Usaco2005 Feb]Rigging the Bovine Election 竞选划区
- bzoj1680[Usaco2005 Mar]Yogurt factory
- BZOJ_P1682 [Usaco2005 Mar]Out of Hay 干草危机(最小生成树)
- bzoj1680[Usaco2005 Mar]Yogurt factory