博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
小R的烦恼 BZOJ3280
阅读量:7079 次
发布时间:2019-06-28

本文共 1982 字,大约阅读时间需要 6 分钟。

分析:

一开始一直Wa,发现是建图建错了,必须得拆点。

S连i,流量为a[i],费用为0,i+n连T,流量同上,费用为0,之后i连i+1费用为0,流量为inf,之后S连n*2+i,流量为li,费用为0,之后枚举j从1到n,n*2+i连接j+n,费用为p[i],之后i连接i+d[i]+1,费用为q[i],流量为inf,之后跑费用流就可以了。

附上代码:

#include 
#include
#include
#include
#include
#include
#include
using namespace std;#define N 155#define S 0#define T 151struct node{ int to,next,val,flow,from;}e[N*200];int head[N],cnt,dis[N],vis[N],from[N],ans,n,m,k;void add(int x,int y,int z,int v){ e[cnt].to=y,e[cnt].next=head[x],e[cnt].val=v; e[cnt].flow=z,e[cnt].from=x,head[x]=cnt++;}void insert(int x,int y,int z,int v){ add(x,y,z,v);add(y,x,0,-v);}int spfa(){ memset(dis,0x3f,sizeof(dis));memset(from,-1,sizeof(from)); queue
q;q.push(S);dis[S]=0; while(!q.empty()) { int x=q.front();q.pop();vis[x]=0; for(int i=head[x];i!=-1;i=e[i].next) { int to1=e[i].to; if(e[i].flow&&dis[to1]>dis[x]+e[i].val) { from[to1]=i; dis[to1]=dis[x]+e[i].val; if(!vis[to1])vis[to1]=1,q.push(to1); } } } return dis[T]==0x3f3f3f3f?0:1;}void mcf(){ int x,i=from[T]; while(i!=-1){x=min(e[i].flow,x);i=from[e[i].from];} i=from[T]; while(i!=-1){e[i].flow-=x,e[i^1].flow+=x,ans+=x*e[i].val,i=from[e[i].from];} return ;}int a[N];int main(){ int t; scanf("%d",&t); for(int cas=1;cas<=t;cas++) { printf("Case %d: ",cas); int v=0;ans=0; memset(head,-1,sizeof(head));cnt=0; scanf("%d%d%d",&n,&m,&k); for(int i=1;i
<<30,0); for(int i=1;i<=n;i++) { int x; scanf("%d",&x);a[i]=x; insert(S,i,x,0); insert(i+n,T,x,0); } for(int i=1;i<=m;i++) { int x,y; scanf("%d%d",&x,&y); insert(S,2*n+i,x,0); for(int j=1;j<=n;j++)insert(2*n+i,j+n,x,y); } for(int i=1;i<=k;i++) { int x,y; scanf("%d%d",&x,&y); for(int j=1;j+x+1<=n;j++) { insert(j,j+n+x+1,1<<30,y); } } while(spfa())mcf(); for(int i=head[T];i!=-1;i=e[i].next) { if(e[i^1].flow!=0) { v=1; puts("impossible"); break; } } if(!v)printf("%d\n",ans); } return 0;}

  

转载于:https://www.cnblogs.com/Winniechen/p/9118888.html

你可能感兴趣的文章
Mysql学习总结(6)——MySql之ALTER命令用法详细解读
查看>>
大型网站技术架构(五)网站高可用架构
查看>>
SVN学习总结(2)——SVN冲突解决
查看>>
MySQL基本
查看>>
<org manual>翻译--3.2 列的宽度与对齐
查看>>
我的友情链接
查看>>
《未测试》如何使用自己已经编译过的lamp安装cacti nagios zabbix
查看>>
SCVMM 2012 R2之为云分配虚拟机
查看>>
中国企业的等级制度
查看>>
pfSense安装和配置pfBlockerNg
查看>>
Ubuntu 14.04 X64 下安装Appcelerator Titanium Studio。
查看>>
R.java was modified manually! Reverting to generated version!(R文件丢失异常原因汇总)
查看>>
ASP.NET MVC Model验证(一)
查看>>
opensuse12.1安装VirtualBox
查看>>
获取系统特殊文件
查看>>
C语言中的for命令
查看>>
Hash算法
查看>>
urlWithString、fileURLWithPath的区别
查看>>
Elasticsearch
查看>>
【ASM内部原理】_asm_kill_unresponsive_clients & _asm_healthcheck_timeout
查看>>