【洛谷】P3980 [NOI2008]志愿者招募
【洛谷】P3980 [NOI2008]志愿者招募
我居然现在才会用费用流解线性规划……
当然这里解决的一类问题比较特殊
以式子作为点,变量作为边,然后要求就是变量在不同的式子里出现了两次,系数一次为+1,一次为-1
这样的话就作为了一个出度和一个入度,和边正好对应了
我们设每种志愿者选择人数是\(x_{i}\)
我们的限制是
x_{i} \geq 0\\
x_{1} + x_{2} \geq A_{1}\\
x_{1} + x_{3} \geq A_{2}\\
x_{3} + x_{4} \geq A_{3}
\end{matrix}\right.
\]
假如数据是这样的吧~
但是我们需要流量平衡,于是我们新加几个变量
x_{i} \geq 0\\
y_{i} \geq 0\\
x_{1} + x_{2} - y_{1} - A_{1} = 0\\
x_{1} + x_{3} - y_{2} - A_{2} = 0\\
x_{3} + x_{4} - y_{3} - A_{3} = 0
\end{matrix}\right.
\]
然后呢我们在前面和后面补上0 = 0
差分后就是
x_{i} \geq 0\\
y_{i} \geq 0\\
x_{1} + x_{2} - y_{1} - A_{1} = 0\\
x_{1} + x_{3} - y_{2} - A_{2} = 0\\
x_{3} + x_{4} - y_{3} - A_{3} = 0
\end{matrix}\right.
\]
x_{i} \geq 0\\
y_{i} \geq 0\\
x_{1} + x_{2} - y_{1} - A_{1} = 0\\
-x_{2} + x_{3} - y_{2} + y_{1} + A_{1} - A_{2} = 0\\
-x_{1} + x_{4} - y_{3} + y_{2} + A_{2} - A_{3} = 0 \\
-x_{3} - x_{4} + y_{3} + A_{3} = 0
\end{matrix}\right.
\]
这样的话就满足我们的限制了
我们对于每个式子建一个点,我们把\(-x_{1}\)的式子向\(x_{1}\)连一条边,容量为正无穷,费用为志愿者的费用
\(y\)同理,不过费用是0
然后再看常数项,如果常数项是负的,则向汇点流一条容量为常数项绝对值的边,如果是正的,则源点向这个式子流一条容量为常数项的边
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define ba 47
#define MAXN 200005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 +c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
struct node {
int to,next,cap;int64 val;
}E[1000005];
int head[1005],sumE = 1,Ncnt,S,T,cur[1005];
int64 ans;
int N,M;
int st[10005],ed[10005],cs[10005],a[1005];
bool vis[1005];
void add(int u,int v,int c,int a) {
E[++sumE].to = v;
E[sumE].next = head[u];
E[sumE].cap = c;
E[sumE].val = a;
head[u] = sumE;
}
void addtwo(int u,int v,int c,int a) {
add(u,v,c,a);add(v,u,0,-a);
}
int64 dis[1005];
bool inq[1005];
queue<int> Q;
bool SPFA() {
for(int i = 1 ; i <= Ncnt ; ++i) dis[i] = 1e18;
dis[S] = 0;Q.push(S);
memset(inq,0,sizeof(inq));
while(!Q.empty()) {
int u = Q.front();Q.pop();
inq[u] = 0;
for(int i = head[u] ; i ; i = E[i].next) {
int v = E[i].to;
if(E[i].cap) {
if(dis[u] + E[i].val < dis[v]) {
dis[v] = dis[u] + E[i].val;
if(!inq[v]) Q.push(v);
}
}
}
}
return dis[T] < 1e18;
}
int dfs(int u,int aug) {
if(u == T) {
ans += aug * dis[T];
return aug;
}
vis[u] = 1;
for(int i = head[u] ; i ; i = E[i].next) {
int v = E[i].to;
if(!vis[v] && E[i].cap && dis[v] == dis[u] + E[i].val) {
int t = dfs(v,min(E[i].cap,aug));
if(t) {
E[i].cap -= t;
E[i ^ 1].cap += t;
return t;
}
}
}
return 0;
}
void MCMF() {
while(SPFA()) {
do {
memset(vis,0,sizeof(vis));
}while(dfs(S,0x7fffffff));
}
}
void Solve() {
read(N);read(M);
S = N + 2;T = N + 3;Ncnt = N + 3;
for(int i = 1 ; i <= N ; ++i) {
read(a[i]);
if(a[i] > a[i - 1]) addtwo(i,T,a[i] - a[i - 1],0);
else addtwo(S,i,a[i - 1] - a[i],0);
addtwo(i,i + 1,0x7fffffff,0);
}
addtwo(S,N + 1,a[N],0);
for(int i = 1 ; i <= M ; ++i) {
read(st[i]);read(ed[i]);read(cs[i]);
addtwo(ed[i] + 1,st[i],0x7fffffff,cs[i]);
}
MCMF();
out(ans);enter;
}
int main(){
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
return 0;
}
【洛谷】P3980 [NOI2008]志愿者招募的更多相关文章
- 洛谷P3980 [NOI2008]志愿者招募
题解 最小费用最大流 每一天是一条边\((inf-a[i], 0)\) 然后对于一类志愿者, 区间两端连一条\((inf, c[i])\) \(S\)向第一个点连\((inf, 0)\) 最后一个点向 ...
- Solution -「NOI 2008」「洛谷 P3980」志愿者招募
\(\mathcal{Description}\) Link. 一项持续 \(n\) 天的任务,第 \(i\) 天需要至少 \(a_i\) 人工作.还有 \(m\) 种雇佣方式,第 \(i\) ...
- P3980 [NOI2008]志愿者招募 费用流 (人有多大胆地有多大产
https://www.luogu.org/problemnew/show/P3980 感觉费用流比网络流的图更难想到,要更大胆.首先由于日期是连续的,所以图中的点是横向排列的. 这道题有点绕道走的意 ...
- P3980 [NOI2008]志愿者招募
思路 巧妙的建图 因为每个志愿者有工作的时段,所以考虑让一个志愿者的流量能够从S流到T产生贡献 所以每个i向i+1连INF-a[x]的边(类似于k可重区间集),每个si向ti连边cap=INF,cos ...
- luogu P3980 [NOI2008]志愿者招募
传送门 网络流又一神仙套路应用 首先考虑列不等式,设\(x_i\)为第i种人的个数,记\(b_{i,j}\)为第i种人第j天是否能工作,那么可以列出n个不等式,第j个为\(\sum_{i=1}^{m} ...
- P3980 [NOI2008]志愿者招募 (费用流)
题意:最多1000天 每天需要至少ai个工人施工 有10000种工人可以雇佣 每种工人可以工作si到ti天 雇佣一个的花费是ci 问怎样安排使得施工花费最少 思考:最直白的建模方式 就是每种工人可以和 ...
- BZOJ 1061: [Noi2008]志愿者招募
1061: [Noi2008]志愿者招募 Time Limit: 20 Sec Memory Limit: 162 MBSubmit: 4064 Solved: 2476[Submit][Stat ...
- BZOJ 1061: [Noi2008]志愿者招募 [单纯形法]【学习笔记】
1061: [Noi2008]志愿者招募 Time Limit: 20 Sec Memory Limit: 162 MBSubmit: 3975 Solved: 2421[Submit][Stat ...
- [BZOJ1061][Noi2008]志愿者招募
[BZOJ1061][Noi2008]志愿者招募 试题描述 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难 题:为即将启动的奥运新项目招募一批短期志愿 ...
随机推荐
- bochs调试命令
Bochs几条基本指令: 通过物理地址查看内存时,可以不加参数'/nuf': 其中n指定显示的单元数,默认是1: u 指定每个显示单元的大小(b表示字节.h表示字(2字节).w表示双字(4字节)),默 ...
- yarn 强制孙依赖的版本
今天博主遇到一个棘手的问题,@vue/cli-service 依赖了一个包 portfiner@^1.0.20,但是 2 天前,这个包更新到了1.0.22,带来了一些问题. 博主第一反应就是想 yar ...
- P5022 旅行
原题链接 https://www.luogu.org/problem/P5022 本着快csp了,做点往年的NOIp的题试试水来着,没想到水这么深 难度还挺大的,耗了我一天的时间(可能是我太菜了) ...
- c 字符串的结束标志
//char carray[]="nihao"; char carray[]={'a','b','c','d','\0'}; printf("array=%s" ...
- 监控zabbix 3.4.11异常通过邮件报警步骤
监控的目的一个是可以查看历史状态,可以对比零晨和工作区间数据的对比,以便后期进行优化指导.还有一个是报警,总不能等到服务器出现异常了才去从头查是什么问题吧.所以这篇主要介绍报警中最基础的一个 配置邮件 ...
- [转] node.js express的安装与部署,以及pm2下的运行启动。
node.js 下的express 安装可以参考官方网站的 http://www.expressjs.com.cn/ 这里主要讲一下一些不同的地方,将我们生成的dist文件夹里的文件放在public里 ...
- Java语言发展史
Java语言发展史 詹姆斯·高斯林(James Gosling)1977年获得了加拿大卡尔加里大学计算机科学学士学位,1983年获得了美国卡内基梅隆大学计算机科学博士学位,毕业后到IBM工作,设计IB ...
- JV默认是如何处理异常
main函数收到这个问题时,有两种处理方式: a:自己将该问题处理,然后继续运行 b:自己没有针对的处理方式,只有交给调用main的jvm来处理 jvm有一个默认的异常处理机制,就将该异常进行处理. ...
- python方法未绑定错误
相信 Python 程序员多多少少都和我一样遇到过 Method Unbound Error,直译过来就是 “方法未绑定错误”,虽然搜索之后知道了使用 @classmethod 这样的装饰起后就可以解 ...
- @Component默认是单例还是多例?
@Component默认是单例还是多例? 答: @Component注解默认实例化的对象是单例,如果想声明成多例对象可以使用@Scope("prototype") @Repos ...