【洛谷】P3980 [NOI2008]志愿者招募

我居然现在才会用费用流解线性规划……

当然这里解决的一类问题比较特殊

以式子作为点,变量作为边,然后要求就是变量在不同的式子里出现了两次,系数一次为+1,一次为-1

这样的话就作为了一个出度和一个入度,和边正好对应了

我们设每种志愿者选择人数是\(x_{i}\)

我们的限制是

\[\left\{\begin{matrix}
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.
\]

假如数据是这样的吧~

但是我们需要流量平衡,于是我们新加几个变量

\[\left\{\begin{matrix}
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

差分后就是

\[\left\{\begin{matrix}
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.
\]

\[\left\{\begin{matrix}
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]志愿者招募的更多相关文章

  1. 洛谷P3980 [NOI2008]志愿者招募

    题解 最小费用最大流 每一天是一条边\((inf-a[i], 0)\) 然后对于一类志愿者, 区间两端连一条\((inf, c[i])\) \(S\)向第一个点连\((inf, 0)\) 最后一个点向 ...

  2. Solution -「NOI 2008」「洛谷 P3980」志愿者招募

    \(\mathcal{Description}\)   Link.   一项持续 \(n\) 天的任务,第 \(i\) 天需要至少 \(a_i\) 人工作.还有 \(m\) 种雇佣方式,第 \(i\) ...

  3. P3980 [NOI2008]志愿者招募 费用流 (人有多大胆地有多大产

    https://www.luogu.org/problemnew/show/P3980 感觉费用流比网络流的图更难想到,要更大胆.首先由于日期是连续的,所以图中的点是横向排列的. 这道题有点绕道走的意 ...

  4. P3980 [NOI2008]志愿者招募

    思路 巧妙的建图 因为每个志愿者有工作的时段,所以考虑让一个志愿者的流量能够从S流到T产生贡献 所以每个i向i+1连INF-a[x]的边(类似于k可重区间集),每个si向ti连边cap=INF,cos ...

  5. luogu P3980 [NOI2008]志愿者招募

    传送门 网络流又一神仙套路应用 首先考虑列不等式,设\(x_i\)为第i种人的个数,记\(b_{i,j}\)为第i种人第j天是否能工作,那么可以列出n个不等式,第j个为\(\sum_{i=1}^{m} ...

  6. P3980 [NOI2008]志愿者招募 (费用流)

    题意:最多1000天 每天需要至少ai个工人施工 有10000种工人可以雇佣 每种工人可以工作si到ti天 雇佣一个的花费是ci 问怎样安排使得施工花费最少 思考:最直白的建模方式 就是每种工人可以和 ...

  7. BZOJ 1061: [Noi2008]志愿者招募

    1061: [Noi2008]志愿者招募 Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 4064  Solved: 2476[Submit][Stat ...

  8. BZOJ 1061: [Noi2008]志愿者招募 [单纯形法]【学习笔记】

    1061: [Noi2008]志愿者招募 Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 3975  Solved: 2421[Submit][Stat ...

  9. [BZOJ1061][Noi2008]志愿者招募

    [BZOJ1061][Noi2008]志愿者招募 试题描述 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难 题:为即将启动的奥运新项目招募一批短期志愿 ...

随机推荐

  1. AtCoder Grand Contest 001 题解

    传送门 \(A\) 咕咕咕 const int N=505; int a[N],n,res; int main(){ scanf("%d",&n); fp(i,1,n< ...

  2. 数据结构实验之图论三:判断可达性(SDUT 2138)(简单DFS)

    #include <bits/stdc++.h> using namespace std; int gra[1002][1005]; int vis[1002]; int n,m; voi ...

  3. 二分算法题目训练(四)——Robin Hood详解

    codeforces672D——Robin Hood详解 Robin Hood 问题描述(google翻译) 我们都知道罗宾汉令人印象深刻的故事.罗宾汉利用他的射箭技巧和他的智慧从富人那里偷钱,然后把 ...

  4. PHP 连接本地mysql

    <?php echo microtime(true); ?> <?php $servername = "localhost"; $username = " ...

  5. 2017杭电ACM集训队单人排位赛 - 2 题解

    1001,水题,直接模拟即可.比赛中开局连wa三发,因为把int写成了bool.. 1002,积分题,比赛中找到了下面这个积分公式, 但是并没什么用,,因为带入以后存在误差,估计是展开了以后出现了误差 ...

  6. zabbix(10)自动发现规则(low level discovery)

    1.概念 在配置Iterms的过程中,有时候需要对类似的Iterms进行添加,这些Iterms具有共同的特征,表现为某些特定的参数是变量,而其他设置都是一样的,例如:一个程序有多个端口,而需要对端口配 ...

  7. PHP 之查找字符串位置函数封装

    /** * 正数查找字符串n次出现的位置 * @param $str * @param $find * @param $n * @return bool|int */ function str_n_p ...

  8. python3编程基础之一:操作

    基本操作有:读数据.写数据.运算.控制.输入.输出.语句块 1.读取数据: num1 = 50 num2 = num1 //通过num2取得num1的值,这就是逻辑上的读取 测试数据:print(nu ...

  9. GO- 使用JSON

    1 json.Marshal  把对象转换为JSON的方法 原型如下 func Marshal(v interface{}) ([]byte, error)这个函数接收任意类型的数据 v,并转换为字节 ...

  10. Flask上下文源码分析(二)

    前面第一篇主要记录了Flask框架,从http请求发起,到返回响应,发生在server和app直接的过程. 里面有说到,Flask框架有设计了两种上下文,即应用上下文和请求上下文 官方文档里是说先理解 ...