【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;
}

【BZOJ1061】【NOI2008】志愿者招募的更多相关文章

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

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

  2. 网络流解线性规划问题 BZOJ1061: [Noi2008]志愿者招募

    线性规划定义: 在给定有限的资源和竞争约束情况下,很多问题都可以表述为最大化或最小化某个目标.如果可以把目标指定为某些变量的线性函数,而且如果可以将资源约束指定为这些变量的等式或不等式,则得到了一个线 ...

  3. 【费用流】BZOJ1061: [Noi2008]志愿者招募(这题超好)

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

  4. [BZOJ1061] [Noi2008] 志愿者招募 (费用流)

    Description 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难 题:为即将启动的奥运新项目招募一批短期志愿者.经过估算,这个项目需要N 天才能 ...

  5. BZOJ1061 [Noi2008]志愿者招募 【单纯形】

    题目链接 BZOJ1061 题解 今天终于用正宗的线性规划\(A\)了这道题 题目可以看做有\(N\)个限制和\(M\)个变量 变量\(x_i\)表示第\(i\)种志愿者的人数,对于第\(i\)种志愿 ...

  6. 【费用流】BZOJ1061[NOI2008]-志愿者招募

    [题目大意] 一个项目需要n天完成,其中第i天至少需要Ai个人.共有m类人可以招募,其中第i类可以从第Si天做到第Ti天,每人的招募费用为Ci元.求最小招募费用. [思路] byvoid神犇的建图详解 ...

  7. BZOJ1061: [Noi2008]志愿者招募(线性规划)

    Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 5725  Solved: 3437[Submit][Status][Discuss] Descript ...

  8. BZOJ1061 NOI2008 志愿者招募 线性规划、费用流

    传送门 一道思路很妙的线性规划网络流 设\(X_i\)表示第\(i\)天需要的人数,\(P_i\)表示第\(i\)种人雇佣的个数 那么我们可以列出一系列式子 比如说样例就可以列出三个式子: \(P_1 ...

  9. 线性规划费用流解法(Bzoj1061: [Noi2008]志愿者招募)

    题面 传送门 Sol 线性规划费用流解法用与求解未知数为非负数的问题 这道题可以列出一堆形如 \(x[i]+x[j]+x[k]+...>=a[p]\) 的不等式 我们强行给每个式子减去一个东西, ...

  10. [BZOJ1061][Noi2008]志愿者招募 线性规划+费用流

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1061 根据题意列方程,然后用网络流解线性规划. 题解直接贴ByVoid的吧,太神了:htt ...

随机推荐

  1. Oracle TRCA 工具(转)

    本篇文章主要介绍了"Oracle TRCA 工具 说明 ",主要涉及到Oracle TRCA 工具 说明 方面的内容,对于Oracle TRCA 工具 说明 感兴趣的同学可以参考一 ...

  2. SQL Server 使用 Hierarchyid 操作层次结构数据

    层次结构数据定义为一组通过层次结构关系互相关联的数据项. 在层次结构关系中,一个数据项是另一个项的父级或子级. sql server2008开始内置的 hierarchyid 数据类型使存储和查询层次 ...

  3. PKU 1379 Run Away(模拟退火算法)

    题目大意:原题链接 给出指定的区域,以及平面内的点集,求出一个该区域内一个点的坐标到点集中所有点的最小距离最大. 解题思路:一开始想到用随机化算法解决,但是不知道如何实现.最后看了题解才知道原来是要用 ...

  4. vue我的总结+转原理

    我的总结 vue:1 mvvm模型,model,view,viewmodel,自底层向上,逐渐增加的模式2 .vue文件 包含html css js 有最近设计原则,将自己需要的放到最近,2 组件化 ...

  5. 2017浙江省赛 H - Binary Tree Restoring ZOJ - 3965

    地址:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3965 题目: iven two depth-first-search ...

  6. Jmeter中解决中文乱码问题

    在使用JMeter过程中,执行结果响应断言总是提示失败,但是使用英文字母可以成功,表示逻辑和其它地方没有问题,问题可能出在编码上,细看了响应数据和日志,其中中文参数被编码成了类似URL编码格式,如下图 ...

  7. 如何在idea中设置 jsp 内容修改以后,立即生效而不用重新启动服务?

    点击 run---->edit configuration--->

  8. C++与C混编

    C++与C混编 本案例通过实现一个简单的UDP服务器来说明C++与C的混合编程问题 C代码 通过C代码来对UDP服务器的创建,监听进行封装 udp.c文件 #include <sys/types ...

  9. Netty资料

    netty 资料  转自   http://calvin1978.blogcn.com/articles/netty-info.html Netty资料皆阵列在前 Posted on 2016-08- ...

  10. Python面试题目之(针对dict或者set数据类型)边遍历 边修改 报错dictionary changed size during iteration

    # result 是一个字典, 把里面属性值是None的属性删除 for key in result: if not result[key]: del result[key] continue 但是报 ...