BZOJ 1061费用流
来源:互联网 发布:胡歌男粉丝知乎 编辑:程序博客网 时间:2024/04/28 15:55
思路:
我们可以列出几个不等式
用y0带进去变成等式
下-上 可以消好多东西
我们发现 等式左边的加起来=0
可以把每个方程看成一个点
正->负 连边
跑费用流即可
//By SiriusRen#include <queue>#include <cstdio>#include <cstring>#include <algorithm>using namespace std;#define int long long#define mem(x,y) memset(x,y,sizeof(x))const int N=24005,M=1005,inf=0x3f3f3f3f;int edge[N],cost[N],v[N],w[N],first[M],next[N],tot,ans;int n,m,xx,yy,zz,least[M],with[M],vis[N],d[M],minn[M];void Add(int x,int y,int C,int E){edge[tot]=E,cost[tot]=C,v[tot]=y,next[tot]=first[x],first[x]=tot++;}void add(int x,int y,int C,int E){Add(x,y,C,E),Add(y,x,-C,0);}bool tell(){ mem(with,0),mem(vis,0),mem(d,0x3f),mem(minn,0x3f); queue<int>q;q.push(0);d[0]=0; while(!q.empty()){ int t=q.front();q.pop();vis[t]=0; for(int i=first[t];~i;i=next[i]){ if(d[v[i]]>d[t]+cost[i]&&edge[i]){ d[v[i]]=d[t]+cost[i],minn[v[i]]=min(minn[t],edge[i]),with[v[i]]=i; if(!vis[v[i]])vis[v[i]]=1,q.push(v[i]); } } }return d[n+1]<inf;}int zeng(){ for(int i=n+1;i;i=v[with[i]^1]) edge[with[i]]-=minn[n+1],edge[with[i]^1]+=minn[n+1]; return minn[n+1]*d[n+1];}signed main(){ mem(first,-1); scanf("%lld%lld",&n,&m); for(int i=1;i<=n;i++){ scanf("%lld",&least[i]); if(least[i]-least[i-1]>0)add(0,i,0,least[i]-least[i-1]); else add(i,n+1,0,least[i-1]-least[i]); if(i!=n)add(i+1,i,0,inf); } for(int i=1;i<=m;i++) scanf("%lld%lld%lld",&xx,&yy,&zz),add(xx,yy+1,zz,inf); while(tell())ans+=zeng(); printf("%lld\n",ans);}
0 0
- BZOJ 1061费用流
- bzoj 1061(线性规划+费用流)
- BZOJ 1061 (费用流,填”坑“ )
- BZOJ 1061: [Noi2008]志愿者招募 费用流
- bzoj-3130 费用流
- bzoj 3130 费用流
- bzoj 2424(费用流)
- bzoj 3171(费用流)
- bzoj 2245(费用流)
- bzoj 1221(费用流)
- bzoj 3280(费用流)
- BZOJ 1070 费用流
- BZOJ 1070 费用流
- BZOJ 3280 费用流
- BZOJ 1283 费用流
- BZOJ 1283: 序列 最大费用费用流
- BZOJ 1283 序列 费用流
- bzoj 2661(最大费用流)
- 解决问题
- 初学Hive数据仓库
- C++中的const
- 3.4 T4
- iOS 图片旋转方法
- BZOJ 1061费用流
- 从TCP三次握手说起——浅析TCP协议中的疑难杂症
- hadoop学习笔记4:hadoop、spark概念
- ubuntu在登录界面添加root用户登录的办法
- git命令记录
- 基于hibernate实现jpa规范
- 在公司的第一个项目总结
- Tomcat日志打印请求响应时间
- BZOJ 1507 splay