题目:

doc 最近太忙了, 每天都有课. 这不怕, doc 可以请假不去上课.
偏偏学校又有规定, 任意连续 n 天中, 不得请假超过 k 天.

doc 很忧伤, 因为他还要陪学姐去逛街呢.

后来, doc发现, 如果自己哪一天智商更高一些, 陪学姐逛街会得到更多的好感度.
现在 doc 决定做一个实验来验证自己的猜想, 他拜托 小岛 预测出了 自己 未来 3n 天中, 每一天的智商.
doc 希望在之后的 3n 天中选出一些日子来陪学姐逛街, 要求在不违反校规的情况下, 陪学姐逛街的日子自己智商的总和最大.

可是, 究竟这个和最大能是多少呢?

格式

输入格式

第一行给出两个整数, n 和 k, 表示我们需要设计之后 3n 天的逛街计划, 且任意连续 n 天中不能请假超过 k 天.
第二行给出 3n 个整数, 依次表示 doc 每一天的智商有多少. 所有数据均为64位无符号整数

输出格式

输出只有一个整数, 表示可以取到的最大智商和.

样例1

样例输入1

5 3
14 21 9 30 11 8 1 20 29 23 17 27 7 8 35

样例输出1

195

限制

对于 20% 的数据, 1 <= n <= 12 , k = 3.
对于 70% 的数据, 1 <= n <= 40 .
对于 100% 的数据, 1 <= n <= 200 , 1 <= k <= 10.

题解:

直接引用官网题解吧····第一次遇到根据等式转化成网络流的····涨知识了··

该题是最大流最大费用流!!!不是最大流最小费用流!!

设第i天是否去逛街为a[i],c[i]表示第i天的智商,a[i]=1表示去逛街,a[i]=0表示不去
则可得2n个不等式
a[1]+a[2]+...+a[n]<=k
a[2]+a[3]+...+a[n+1]<=k
...
a[2n+1]+....+a[3n]<=k
求c[1]*a[1]+c[2]*a[2]+...+c[3n]*a[3n]的最大值
添加一个辅助变量
a[1]+a[2]+...+a[n]+y[1]=k
a[2]+a[3]+...+a[n+1]+y[2]=k
...
a[2n+1]+....+a[3n]+y[2n+1]=k
0<=y[i]<=k
将上述不等式相邻两个相减
y[1]+a[1]=a[n+1]+y[2]。。。。。。。。。1
y[2]+a[2]=a[n+2]+y[3]。。。。。。。。。2
......
y[n+1]+a[n+1]=a[2n+1]+y[n+2]。。。。。。。n+1
......
y[2n]+a[2n]=a[3n]+y[2n+1]。。。。。。。。2n
根据网络中每个节点流入量等于流出量的性质
将上述等式编号并抽象成网络中的点,变量a[i]和y[i]抽象为网络中的有向边(弧)
问题等价于求最大费用最大流
以a[n+1]为例 可以看成是节点1部分流出量和节点n+1的部分流入量于是可以建边从n+1到1
故根据这些等式可以建图
设源点为0,汇点为2n+3
i到n+i连一条弧,流量上限为1,费用为c[n+i] 1<=i<=n 
i到i+1连一条弧,流量上限为k,费用为0(即为辅助变量y)1<=i<=2n-1
这时发现题目的k还没用上,
于是发现上述等式成立必需满足这两个等式
a[1]+a[2]+...+a[n]+y[1]=k。。。。。2n+1
a[2n+1]+....+a[3n]+y[2n+1]=k。。。。。2n+2
于是建一个节点2n+1
为了满足2n+1式
则由源点向节点2n+1连一条流量上限为k的边,费用为0。
由节点2n+1向i连一条流量上限为1的边,费用为c[i](即为变量a[i])1<=i<=n
由2n+1向1连一条流量上限为k的边,费用为0(即y[1])
同理建一个节点2n+2
为了满足2n+2式则由节点2n+2向汇点连一条流量上限为k的边,费用为0。
由节点i向2n+2连一条流量上限为1的边,费用为c[i+n]  n+1<=i<=2n;

由2n向2n+1连一条流量上限为k的边,费用为0(及y[2n+1])

代码:

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<cctype>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
using namespace std;
const int N=;
const int M=;
queue<int>que;
int n,k;
int tot=,fst[N*+],nxt[M],go[M],rest[M],src,des;
long long c[N*],cost[M],ans=,dis[N];
bool visit[N],work[N];
inline void comb(int a,int b,int c,long long d)
{
nxt[++tot]=fst[a],fst[a]=tot,go[tot]=b,rest[tot]=c,cost[tot]=d;
nxt[++tot]=fst[b],fst[b]=tot,go[tot]=a,rest[tot]=,cost[tot]=-d;
}
inline void build()
{
for(int i=;i<=n;i++)
comb(i,n+i,,c[n+i]);
for(int i=;i<=*n-;i++)
comb(i,i+,k,); comb(src,*n+,k,);
for(int i=;i<=n;i++)
comb(*n+,i,,c[i]);
comb(*n+,,k,); comb(*n+,des,k,);
for(int i=n+;i<=*n;i++)
comb(i,*n+,,c[i+n]);
comb(*n,*n+,k,);
}
inline bool SPFA()
{
for(int i=src;i<=des;i++) dis[i]= -1e+, work[i]=false;
que.push(src),dis[src]=;
int u,v;
while(!que.empty())
{
u=que.front();
que.pop();
visit[u]=false;
for(int e=fst[u];e;e=nxt[e])
{
if(dis[v=go[e]] < dis[u]+cost[e]&&rest[e])
{
dis[v]=dis[u]+cost[e];
if(!visit[v])
{
visit[v]=true;
que.push(v);
}
}
}
}
return dis[des]!= -1e+;
}
inline int dinic(int u,int flow)
{
if(u==des)
{
ans+=dis[des]*flow;
return flow;
}
work[u]=true;
int delta=,res=,v;
for(int e=fst[u];e;e=nxt[e])
{
if(dis[v=go[e]]==dis[u]+cost[e]&&rest[e]&&!work[v])
{
delta=dinic(v,min(flow-res,rest[e]));
if(delta)
{
rest[e]-=delta;
rest[e^]+=delta;
res+=delta;
if(res==flow) break;
}
}
}
return res;
}
void maxflow()
{
while(SPFA())
dinic(src,);
}
int main()
{
//freopen("a.in","r",stdin);
scanf("%d%d",&n,&k);
for(int i=;i<=n*;i++)
scanf("%d",&c[i]);
src=,des=*n+;
build();
maxflow();
cout<<ans<<endl;
return ;
}

刷题总结——学姐的逛街计划(vijos1891费用流)的更多相关文章

  1. Vijos1891 学姐的逛街计划 【费用流】*

    Vijos1891 学姐的逛街计划 描述 doc 最近太忙了, 每天都有课. 这不怕, doc 可以请假不去上课. 偏偏学校又有规定, 任意连续 n 天中, 不得请假超过 k 天. doc 很忧伤, ...

  2. vijos1891 学姐的逛街计划(线性规划)

    P1891学姐的逛街计划 描述 doc 最近太忙了, 每天都有课. 这不怕, doc 可以请假不去上课.偏偏学校又有规定, 任意连续 n 天中, 不得请假超过 k 天. doc 很忧伤, 因为他还要陪 ...

  3. [vijos1891]学姐的逛街计划

                                                                     学姐的逛街计划 描述 doc 最近太忙了, 每天都有课. 这不怕, d ...

  4. Vijos 学姐的逛街计划

    传送门 题解传送门 线性规划,最小费用最大流. 神奇的操作. //Achen #include<algorithm> #include<iostream> #include&l ...

  5. python部落刷题宝学到的内置函数

    最近加入了python部落,感觉里面的刷题宝很有意思,玩了一下,知道了许多以前并不清楚的内置函数,然后感觉到快要记不住了,所以开始陈列一下 1.divmod(a,b):取a除以b的商和余数,功效等价于 ...

  6. python部落刷题宝学到的内置函数(二)

    感觉到刷题宝有一个好处,也许也不是好处,它的答案必须是真正输出的值,也就是说应该输出字符串aaaa的时候,答案必须写成界面上返回的值,即'aaaa'.有利于真正记忆返回值类型,但是....太繁琐了 1 ...

  7. 牛客网Java刷题知识点之输入流、输出流、字节流、字符流、字节流的抽象基类(InputStream、OutputStream)、字符流的抽象基类(Reader、Writer)、FileWriter、FileReader

    不多说,直接上干货! IO流用来处理设备之间的数据传输. java对数据的操作是通过流的方式. java用于操作流的对象都在IO包中. IO流按操作数据分为两种:字节流和字符流. IO流按流向分为:输 ...

  8. [网络流24题] 最长k可重区间集问题 (费用流)

    洛谷传送门 LOJ传送门 很巧妙的建图啊...刚了$1h$也没想出来,最后看的题解 发现这道题并不类似于我们平时做的网络流题,它是在序列上的,且很难建出来二分图的形. 那就让它在序列上待着吧= = 对 ...

  9. 【网络流24题】No.18 分配问题 (二分图最佳匹配 费用流|KM)

    [题意] 有 n 件工作要分配给 n 个人做.第 i 个人做第 j 件工作产生的效益为 cij . 试设计一个将n 件工作分配给 n 个人做的分配方案, 使产生的总效益最大. 输入文件示例input. ...

随机推荐

  1. Java、Node.js、PHP还是.Net? 无论你选谁,我都能教你一招!

    七夕如期而至,不该来的终究还是来了.再傲娇的单身贵族恐怕也难免在今天会感觉一丝丝的空虚.还好你关注了我,因为接下来我准备了三大招教你一个人…..也可以优雅地过七夕. 招式一:移形幻影,无中生有 七夕当 ...

  2. eclipse 在写XML时 包类名自动提醒的问题

    需要加一个STS插件 配置很简单 参考了  https://blog.csdn.net/HH775313602/article/details/70176531 在 https://spring.io ...

  3. 响应式Web设计- Viewport

    什么是Viewport? viewport是用户网页的可视区域, 翻译为中文可以叫做"视区". 设置Viewport 一个常用的针对移动网页优化过的页面的Viewport meta ...

  4. React初识整理(一)

    一.React的特点 1.自动化的UI状态管理:自动完成数据变化与界面效果的更新. 2.虚拟DOM:创建1个虚拟的dom节点树,放在内存里(内存修改数据效率高),数据变化时先修改内存里的虚拟DOM,然 ...

  5. Ubuntu创建应用快捷方式

    Ubuntu创建应用快捷方式 新建一个.desktop文件 vi eclipse.desktop 然后又进行编辑 [Desktop Entry] Encoding=UTF-8 Name=eclipse ...

  6. 如何优化sql查询

    借鉴https://www.cnblogs.com/ssrstm/p/5753068.html和https://www.cnblogs.com/exe19/p/5786806.html 1. 对查询进 ...

  7. APIO2019&&THUSC2019游记

    APIO2019懵十三记: day0: 早上和ljx从沈阳出发,下午一点到的首师大附. 由于工作人员中午十二点就散了,我们就先去试机了. 下午三点接到狗牌和T恤,晚上买麦当劳回如意吃. 晚上还有场模拟 ...

  8. golang 实现冒泡排序

    package main import ( "fmt" ) func main(){ a := [...] int{2,5,9,6,8} fmt.Println(a) num := ...

  9. CSS3-弹性盒模型-FlexBox

    Flex容器属性 display 定义一个Flex容器,根据其取的值来决定是内联还是块.Flex容器会为其内容建立新的伸缩格式化上下文. .container { display: flex; /* ...

  10. perl中foreach(一)

    perl中的foreach结构  首先语法 foreach $rock(qw /bedrock slate lava/){        $rock="\t$rock";      ...