Delight for a Cat
Time Limit: 1000 ms Memory Limit: 512 MB
Description
从前,有一只懒猫叫CJB。每个小时,这只猫要么在睡觉,要么在吃东西,但不能一边睡觉一边吃东西,并且这只猫会在一整个小时干同一件事情。
对于接下来的n个小时,CJB知道他在那n个小时睡觉和吃东西的快乐值。
为了健♂康♀地生活,在任意的连续 k个整小时内,CJB要有至少\(m_s\)小时睡觉,至少\(m_e\)个小时在吃东西。也就是说一共有 n−k+1段的 k小时需要满足上述条件。
你的任务是告诉CJB他接下来n个小时能获得的最大快乐值是多少。
Input
第一行四个整数 \(n,k,ms,me(1≤k≤n≤1000,0≤m_s,m_e≤k,m_s+m_e≤k)\),含义如题目所述。
第二行包含 n个整数 \(s_1,s_2,...,s_n(0≤s_i≤10^9)\),表示每个小时CJB睡觉能获得的快乐值。
第三行包含 n个整数 \(e_1,e_2,...,e_n(0≤e_i≤10^9)\),表示每个小时CJB吃东西能获得的快乐值。
Output
输出一个整数,表示CJB能获得的最大快乐值。
Sample Input
10 4 1 2
1 2 3 4 5 6 7 8 9 10
10 9 8 7 6 5 4 3 2 1
Sample Output
69
HINT
【样例解释】
CJB在接下来的10个小时是这样的(S表示睡觉,E表示吃东西):EEESESEESS
【数据范围与约定】
30%的数据中, 1≤n≤10
100%的数据中, 1≤n≤1000
吐槽
这题好神啊哈哈哈qwq(好吧虽然好像印象中之前做过类似的费用流的题目。。但是当时并没有搞太懂所以看到这题的时候根本没有往那边去想qwq菜醒)
然而为什么网上说是。。单纯形裸题??好吧听起来很棒但是并不会单纯形qwq
所以还是写费用流吧嗯qwq
正题
这里要用到一种特殊的建模方式!十分巧妙啊ovo
首先假设这只懒猫一直在睡觉,然后我们再决策哪些天数需要改成吃东西
用\(x_i =0\)表示这只猫第\(i\)天在睡觉,\(x_i=1\)表示在吃东西,那么我们可以列出这样的不等式组
&m_e<=x_1+x_2+x_3+...+x_k<=k-m_s\\
&m_e<=x_2+x_3+x_4+...+x_{k+1}<=k-ms\\
&...\\
&m_e<=x_{n-k+1}+x_{n-k+2}+x_{n-k+3}+...+x_n<=k-m_s\\
\end{aligned}
\]
用同样的方法将这些不等式组变成等式,其中辅助变量\(0<=y_i<=k-m_s-m_e\)
&x_1+x_2+x_3+...+x_k+y_1=k-m_s\\
&x_2+x_3+x_4+...+x_{k+1}+y_2=k-m_s\\
&...\\
&x_{n-k+1}+x_{n-k+2}+...+x_n=k-m_s\\
\end{aligned}
\]
同样的套路,在一头一尾加上两个0=0等式,用上面的减去下面的得到:
&x_1+x_2+x_3+...+x_k+y_1=k-m_s\\
&x_1-x_{k+1}+y_1-y_2=0\\
&x_2-x_{k+2}+y_2-y_3=0\\
&...\\
&x_{n-k}-x_{n}+y_{n-k}-y_{n-k+1}=0\\
&x_{n_k+1}+x_{n-k+2}+...+x_{n}+y_{n-k+1}=k-m_s\\
\end{aligned}
\]
(啊?想用下面的减去下面的?额一样的反正都可以)
然后按照套路:
&x_1-x_{k+1}+y_1-y_2=0\\
&x_2-x_{k+2}+y_2-y_3=0\\
&...\\
&x_{n-k}-x_{n}+y_{n-k}-y_{n-k+1}=0\\
&(k-m_s)-x_1-x_2-x_3-...-x_k-y_1=0\\
&x_{n_k+1}+x_{n-k+2}+...+x_{n}+y_{n-k+1}-(k-m_s)=0\\
\end{aligned}
\]
(其实倒数第二条式子是第1条式子和第0条式子相减得到,的但是为了和后面的代码中的编号对应我们稍微改变一下位置)
观察一下,变成了这个样子之后所有式子加起来变量和常数项就都抵消掉了,满足流量平衡
那么就用和bzoj1061同样的方式把每个等式看成一个点,负的流出正的流入就好啦
跑完一遍最小费用最大流之后,最终的\(ans\)就是$\sum\limits_{i=1}^{n}s_i - $最小费用
然而这并没有结束
首先是,这题的\(x_i\)和\(y_i\)是有范围的!
\(y_i\)的范围上面写了就是\(0<=y_i<=k-m_s-m_e\)
而\(x_i\)的话要么是\(0\)要么是\(1\)
所以两个变量对应的边的容量应该分别为\(k-m_s-m_e\)和\(1\)
再有就是,如果稍微想一下等式之间的连边情况,会发现。。\(k<n-k\)和\(k>n-k\)这两种情况是不一样的!
所以说,这题还要分类讨论
然而具体的连边方式。。写起来十分混乱不如自己去推qwq,详情可以参考代码
(注意,代码中将上面的第一条等式放在了倒数第二条,所以编号为\(n-k+1\),最后一条等式的编号为\(n-k+2\),其他的等式依次为\(1\)到\(n-k\))
这里只讲一下大概要连哪几类边,下面所说的编号顺序按照最终给出的等式组来(也就是代码中的顺序):
- 第一种\(y\)边。从\(i+1\)连到\(i\)
- 第二种\(y\)边。从最后一条等式连到第\(1\)条等式,从\(n-k\)条等式连到倒数第二条等式
- 源点到\(n-k+1\)的边,汇点到\(n-k+2\)的边
- 第\(1\)到\(n-k\)条等式中每个正的\(x\)的边。
- 第\(1\)到\(n-k\)条等式中每个负的\(x\)的边。
稍微提一下,前3类是固定的,第4和第5类又可以分第\(1\)到\(n-k\)条连到最后两条,和\(1\)到\(n-k\)条里面自己乱连两小类,具体就\(k\)和\(n\)的大小关系而定
自己写的时候把\(k=n-k\)的情况也分出来了,可能是。。可以合并的不过不管了反正这么写也不影响正确性
然后就很愉快(个鬼啊分类讨论这种东西)滴做完啦ovo
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define ll long long
#define F 1
using namespace std;
const int N=10010;
const int M=1000010;
const ll inf=1e15;
struct xxx{
int y,next,op;
ll c,r;
}a[M*10];
queue<int> q;
int h[N];
ll e[N],s[N];
ll pre[M],mn[N],pred[M];
ll d[N];
ll ans,sum;
bool vis[N];
int n,m,vs,vt,tot;
ll k,ms,me;
int add(int x,int y,int r,ll c);
int spfa();
int build();
int main(){
#ifndef ONLINE_JUDGE
freopen("a.in","r",stdin);
// freopen("a.out","w",stdout);
#endif
scanf("%d%lld%lld%lld",&n,&k,&ms,&me);
vs=0; vt=n-k+3;
sum=0;
for (int i=1;i<=n;++i) scanf("%lld",s+i),sum+=s[i];
for (int i=1;i<=n;++i) scanf("%lld",e+i);
build();
ans=0;
while (spfa());
printf("%lld\n",sum+ans);
}
int add(int x,int y,int r,ll c){
a[++tot].y=y; a[tot].next=h[x]; h[x]=tot; a[tot].r=r; a[tot].op=tot+1; a[tot].c=c;
a[++tot].y=x; a[tot].next=h[y]; h[y]=tot; a[tot].r=0; a[tot].op=tot-1; a[tot].c=-c;
}
int spfa(){
while (!q.empty()) q.pop();
for (int i=vs;i<=vt;++i) d[i]=-inf,vis[i]=false,mn[i]=inf;
q.push(vs); vis[vs]=1; d[vs]=0;
mn[vs]=inf; pre[vs]=-1; pred[vs]=-1;
int v,u;
while (!q.empty()){
v=q.front(); q.pop();
for (int i=h[v];i!=-1;i=a[i].next){
u=a[i].y;
if (!a[i].r) continue;
if (d[u]<d[v]+a[i].c){
pre[u]=i; pred[u]=v;
mn[u]=min(mn[v],a[i].r);
d[u]=d[v]+a[i].c;
if (!vis[u])
vis[u]=1,q.push(u);
}
}
vis[v]=false;
}
if (d[vt]==-inf) return false;
ll flow=mn[vt];
ans+=flow*d[vt];
u=vt;
while (pre[u]!=-1){
a[pre[u]].r-=flow;
a[a[pre[u]].op].r+=flow;
u=pred[u];
}
return true;
}
int build(){
memset(h,-1,sizeof(h));
if (k<n-k){
/*X*/
for (int i=1;i<=k;++i) add(n-k+1,i,F,e[i]-s[i]);
for (int i=n-2*k+1;i<=n-k;++i) add(i,n-k+2,F,e[i+k]-s[i+k]);
for (int i=1;i<=n-2*k;++i) add(i,i+k,F,e[i+k]-s[i+k]);
}
else if (k==n-k){
/*X*/
for (int i=1;i<=n-k;++i) add(n-k+1,i,F,e[i]-s[i]);
for (int i=1;i<=n-k;++i) add(i,n-k+2,F,e[i+k]-s[i+k]);
}
else{
/*X*/
for (int i=1;i<=n-k;++i) add(n-k+1,i,F,e[i]-s[i]);
for (int i=1;i<=n-k;++i) add(i,n-k+2,F,e[i+k]-s[i+k]);
for (int i=n-k+1;i<=k;++i) add(n-k+1,n-k+2,F,e[i]-s[i]);
}
/*Y*/for (int i=1;i<=n-k-1;++i) add(i,i+1,k-ms-me,0);
/*Y*/add(n-k+1,1,k-ms-me,0); add(n-k,n-k+2,k-ms-me,0);
/*X*/add(vs,n-k+1,k-ms,0); add(n-k+2,vt,k-ms,0);
}
Delight for a Cat的更多相关文章
- 【BZOJ4842】[Neerc2016]Delight for a Cat 线性规划+费用流
[BZOJ4842][Neerc2016]Delight for a Cat Description ls是一个特别堕落的小朋友,对于n个连续的小时,他将要么睡觉要么打隔膜,一个小时内他不能既睡觉也打 ...
- [bzoj4842][bzoj1283][Neerc2016]Delight for a Cat/序列_线性规划_费用流
4842: [Neerc2016]Delight for a Cat_1283: 序列 题目大意:ls是一个特别堕落的小朋友,对于n个连续的小时,他将要么睡觉要么打隔膜,一个小时内他不能既睡觉也打隔膜 ...
- 【xsy2115】Delight for a Cat
Time Limit: 1000 ms Memory Limit: 512 MB Description 从前,有一只懒猫叫CJB.每个小时,这只猫要么在睡觉,要么在吃东西,但不能一边睡觉一边吃东 ...
- BZOJ.4842.[NEERC2016]Delight for a Cat(费用流)
BZOJ 参考这儿. 首先如果一个活动的时间满足条件,那么另一个活动也一定满足.还有就是这题就是费用流没有为什么.不妨假设最初所有时间都用来睡觉,那么我们要对每个\(k\)大小区间选出\([t2,k- ...
- bzoj 4842 [Neerc2016]Delight for a Cat 最小费用最大流,线性规划
题意:有n个小时,对于第i个小时,睡觉的愉悦值为si,打隔膜的愉悦值为ei,同时对于任意一段连续的k小时,必须至少有t1时间在睡觉,t2时间在打隔膜.如果要获得的愉悦值尽 量大,求最大的愉悦值和睡觉还 ...
- bzoj4842 Delight for a Cat
题意:n天内你每天可以s或者e,分别有一定的收益. 每连续k天中s的天数要大于ds,e的天数要大于de,求最大收益. 解:费用流解线性规划. 先假设全部选e,然后一天s的收益为si - ei ai表示 ...
- 题解-bzoj1283序列 & bzoj4842 [Neerc2016]Delight for a Cat
因为这两题有递进关系,所以放一起写 Problem bzoj1283 序列 题意概要:一个长度为 \(n\) 的序列\(\{c_i\}\),求一个子集,使得原序列中任意长度为 \(m\) 的子串中被选 ...
- [BZOJ4842]Delight for a Cat[费用流]
题意 题目链接 分析 类似 最长k可重区间集 一题. 由于本题区间长度相同,首先可以将点的影响看成区间,区间看成点. 先默认所有位置选择事件2,选择区间看做改选事件1 .于是问题变成了求收益最大的方案 ...
- bzoj 4842: [Neerc2016]Delight for a Cat
Description ls是一个特别堕落的小朋友,对于n个连续的小时,他将要么睡觉要么打隔膜,一个小时内他不能既睡觉也打隔膜 ,因此一个小时内他只能选择睡觉或者打隔膜,当然他也必须选择睡觉或打隔膜, ...
- bzoj4842: [Neerc2016]Delight for a Cat
bzoj4842 这是一道网络流的题(大家都看出来了吧) 首先我们简化一下题目,选出最关键的部分(就是知道什么和要求什么,还有条件) 我们在这里把睡觉设为0,至少有t0时间在睡觉,把打隔膜设为1,至少 ...
随机推荐
- hdu2066一个人的旅行(floyd优化)
一个人的旅行 Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Subm ...
- 安装SQLSEVER与MySQL
昨天装了一整填的SQLSEVER,今天是把昨天遗留的问题给重新整合一下,上午安装MySQL的时候,是在网上找的帖子通过压缩包安装的,捣鼓了一上午,下午花不到一个小时, 去安装好了,我觉得通过压缩包安装 ...
- Linux命令应用大词典-第9章 数字计算
9.1 bc:任意精度的计算器 9.2 dc:一个任意精度的计算器 9.3 expr:将表达式的值打印到标准输出 9.1 bc:任意精度的计算器 9.2 dc:一个任意精度的计算器 9.3 expr: ...
- 【Set jsonObj = toJson( jsonString )】创建JSON实例
创建JSON实例: 原型: toJson( jsonString ) 说明: 创建JSON实例 返回: [JSON] 参数: jsonString [可选] 可以用json格式字符串创建实例 示例: ...
- 孤荷凌寒自学python第八十一天学习爬取图片1
孤荷凌寒自学python第八十一天学习爬取图片1 (完整学习过程屏幕记录视频地址在文末) 通过前面十天的学习,我已经基本了解了通过requests模块来与网站服务器进行交互的方法,也知道了Beauti ...
- 孤荷凌寒自学python第七十七天开始写Python的第一个爬虫7
孤荷凌寒自学python第七十七天开始写Python的第一个爬虫7 (完整学习过程屏幕记录视频地址在文末) 今天在上一天的基础上继续完成对我的第一个代码程序的书写. 今天的学习仍然是在纯粹对docx模 ...
- 【MySQL解惑笔记】Centos7下卸载彻底MySQL数据库
彻底卸载Yum安装的MySQL数据库 在我第二章MySQL数据库基于Centos7.3-部署过程中,因为以前安装过其它的版本所以没有卸载干净影响后期安装 一.卸载Centos7自带的Maridb数据库 ...
- 使用eclipse创建maven项目出现的一个问题
错误信息 This error occurs when you employ a plugin that Maven could not download. Possible causes for t ...
- Python3 小工具-TCP半连接扫描
from scapy.all import * import optparse import threading def scan(ip,port): pkt=IP(dst=ip)/TCP(dport ...
- vmware中三种网络连接方式
原文来自http://note.youdao.com/share/web/file.html?id=236896997b6ffbaa8e0d92eacd13abbf&type=note vmw ...