【BZOJ1061】【NOI2008】志愿者招募
【BZOJ1061】【NOI2008】志愿者招募
题面
BZOJ
题解
我们设每类志愿者分别招募了\(B[i]\)个
那么,我们可以得到一系列的方程
\[\sum_{S[i]\leq x\leq T[i]}B[i]\geq A[x]\]
一共\(n\)天,所以是\(n\)个方程
因为有不等式让我们很不爽,所以换成等式。
这里,设\(y[i]\)表示第\(i\)天多招募的人数
这样,方程变成了
\[(\sum_{S[i]\leq x\leq T[i]}B[i])-y[i]=A[x]\]
所以我们得到了一个方程组
\[\begin{cases}0&=0\\\sum_{S[i]\leq 1\leq T[i]}B[i]-y[1]&=A[1]\\ \sum_{S[i]\leq 2\leq T[i]}B[i]-y[2]&=A[2]\\...\\\sum_{s[i]\leq n\leq T[i]}B[i]-y[n]&=A[n]\\0&=0\end{cases}\]
对于每一类志愿者,他一定出现在连续的一段
如果我们对于相邻的两个方程做差,得到\(n-1\)个方程,
这样每一类志愿者就只会在两个方程中出现过
所以,方程组变成了
\[(\sum_{S[i]=x}B[i]-\sum_{T[i]+1=x}B[i])-y[x]+y[x-1]=A[x]-A[x-1]\]
这样一个等式很像流量平衡
所以我们考虑如何连边?
我们把\(+\)看成流入,\(-\)看成流出
这样子相当于每天就是一个点
前面两个\(\sum\)的差,相当于从\(S[i]\)流入,向\(T[i]+1\)流出
也就是如果这个点作为一类志愿者的\(S[i]\)那么他就会流入一定的流
如果作为一个志愿者的\(T[i]+1\),那么它就要流出去,显然是流到对应的\(S[i]\)去
而\(-y[x]+y[x-1]\)相当于从前一天多出来的志愿者可以流入,然后可以向后一天流出
右边如果非零的话,可以看做从源(汇)点流入(出)
如果大于零则证明流入多于流出,则多出来的流流向汇点
如果小于零则证明流出多于流入,则多出来的流从源点流入
这样建图就可以跑费用流了
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 1111
#define inf 1000000000
inline int read()
{RG int x=0,t=1;RG char ch=getchar();while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();if(ch=='-')t=-1,ch=getchar();while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();return x*t;
}
int n,m,a[MAX];
struct Line{int v,next,w,fy;}e[50000];
int h[MAX],cnt=2;
inline void Add(int u,int v,int w,int fy)
{e[cnt]=(Line){v,h[u],w,fy};h[u]=cnt++;e[cnt]=(Line){u,h[v],0,-fy};h[v]=cnt++;
}
bool vis[MAX];
int S,T;
ll Cost,dis[MAX];
bool SPFA(int S,int T)
{memset(vis,0,sizeof(vis));memset(dis,63,sizeof(dis));queue<int> Q;Q.push(S);dis[S]=0;vis[S]=true;while(!Q.empty()){int u=Q.front();Q.pop();for(int i=h[u];i;i=e[i].next){if(!e[i^1].w)continue;int v=e[i].v;if(dis[v]>dis[u]-e[i].fy){dis[v]=dis[u]-e[i].fy;if(!vis[v])vis[v]=true,Q.push(v);}}vis[u]=false;}if(dis[T]>=1e18)return false;return true;
}
int dfs(int u,int flow)
{if(u==T||!flow)return flow;int ret=0;vis[u]=true;for(int i=h[u];i;i=e[i].next){int v=e[i].v;if(!vis[v]&&e[i].w&&dis[v]==dis[u]-e[i].fy){int d=dfs(v,min(flow,e[i].w));e[i].w-=d;e[i^1].w+=d;flow-=d;ret+=d;if(!flow)break;}}return ret;
}
int main()
{n=read();m=read();S=0;T=n+2;for(int i=1;i<=n;++i)a[i]=read();for(int i=1;i<=n+1;++i){if(a[i]-a[i-1]>=0)Add(i,T,a[i]-a[i-1],0);else Add(S,i,a[i-1]-a[i],0);if(i<=n)Add(i,i+1,inf,0);}for(int i=1;i<=m;++i){int S=read(),T=read(),C=read();Add(T+1,S,inf,C);}while(SPFA(T,S))Cost+=1ll*dis[S]*dfs(S,inf);printf("%lld\n",Cost);return 0;
}
转载于:https://www.cnblogs.com/cjyyb/p/9087319.html
【BZOJ1061】【NOI2008】志愿者招募相关推荐
- [BZOJ1061][Noi2008]志愿者招募
[BZOJ1061][Noi2008]志愿者招募 试题描述 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难 题:为即将启动的奥运新项目招募一批短期志愿 ...
- BZOJ1061: [Noi2008]志愿者招募
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1061 1061: [Noi2008]志愿者招募 Time Limit: 20 Sec M ...
- 【费用流】BZOJ1061: [Noi2008]志愿者招募(这题超好)
1061: [Noi2008]志愿者招募 Time Limit: 20 Sec Memory Limit: 162 MB Submit: 5291 Solved: 3173 [Submit][St ...
- [BZOJ1061] [NOI2008] 志愿者招募 - 最小费用最大流
大部分内容转自: BYVOID - NOI2008 志愿者招募 如果讲道理的话,就是说我们抽象一下这个模型--然后每条费用边就是连接起始日期和结束日期的边,也就是说这条边上的流量增加1,就要增加一个 ...
- BZOJ1061 [NOI2008]志愿者招募
Description 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难 题:为即将启动的奥运新项目招募一批短期志愿者.经过估算,这个项目需要N 天才能 ...
- [bzoj1061] [NOI2008]志愿者招募
Description 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难题:为即将启动的奥运新项目招募一批短期志愿者.经过估算,这个项目需要N 天才能完 ...
- BZOJ1061: [Noi2008]志愿者招募(线性规划)
Time Limit: 20 Sec Memory Limit: 162 MB Submit: 5725 Solved: 3437 [Submit][Status][Discuss] Descri ...
- [Bzoj1061][Noi2008]志愿者招募(费用流)
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1061 一开始疯狂想dp,然后队友走过来瞄一眼就告诉我像费用流,菜的真实,jpg. 一种比 ...
- BZOJ1061 NOI2008 志愿者招募 线性规划、费用流
传送门 一道思路很妙的线性规划网络流 设\(X_i\)表示第\(i\)天需要的人数,\(P_i\)表示第\(i\)种人雇佣的个数 那么我们可以列出一系列式子 比如说样例就可以列出三个式子: \(P_1 ...
- 【BZOJ1061/3265】[Noi2008]志愿者招募/志愿者招募加强版 单纯形法
[BZOJ1061][Noi2008]志愿者招募 Description 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难题:为即将启动的奥运新项目招募 ...
最新文章
- 如何把文件隐藏在一张图片里面
- hdu4067 费用流(混合欧拉的宽展和延伸)
- CrazePony飞行器--相关资料网址
- springMVC请求发生重复路径
- Spread for Windows Forms快速入门(10)---绑定到数据库
- 树莓派IO口驱动代码的编写、微机总线地址、物理地址、虚拟地址、BCM2835芯片手册
- NIO学习–核心概念与基本读写
- 用于连接mysql的java类_【考试】列举Java连接数据库用到的类有哪些。
- 掌握穷变富的12条原则 迅速从普通人变成有钱人
- Django使用mysql连接池_Django db使用MySQL连接池
- 浏览器兼容console对象的简要解决方案
- 洛谷 - P3935 - Calculating - 整除分块
- 移动端click事件延迟300ms到底是怎么回事,该如何解决?
- WinForm TabControl美化
- 请你来说一下数字证书是什么,里面都包含那些内容
- ubuntu修改u盘权限_Ubuntu下U盘只读文件系统,图标上锁,提示无法修改
- .net 后端生成海报
- 面向对象在python中的实践---如何定义类、对象?
- 除了神经网络其他的智能算法,人工神经网络算法与机器算法
- gmssl国密双证书制作
热门文章
- 已经被美团、神策“问烂”了的三道多线程面试题,你能答对几道?
- python 计量做hp滤波_R语言提取时间序列的周期性成分应用EMD,小波滤波器,Baxter过滤器等...
- 微信小程序—for循环包括自定义for循环中的item和index(图文)
- 新宝总结过年红包行情来了
- Python编程基础课后答案
- a标签有关用法以及锚点定位;
- 进程控制——创建,终止
- 堂堂中华网china.com上竟然全是垃圾小广告
- matplotlib绘制多个雷达图(包含多个子图):给出四大国有银行2018年相关财务和监管指标的排名,要求将4家银行各项指标的排名绘制成雷达图并且以2x2子图的形式显示
- 常见人体铰链关节点数据集中的关节点顺序(SMPL,NTU,MPII,Human3.6M)