Description###

Kiana最近喜欢到一家非常美味的寿司餐厅用餐。每天晚上,这家餐厅都会按顺序提供n种寿司,第i种寿司有一个

代号ai和美味度di,i,不同种类的寿司有可能使用相同的代号。每种寿司的份数都是无限的,Kiana也可以无限次

取寿司来吃,但每种寿司每次只能取一份,且每次取走的寿司必须是按餐厅提供寿司的顺序连续的一段,即Kiana

可以一次取走第1,2种寿司各一份,也可以一次取走第2,3种寿司各一份,但不可以一次取走第1,3种寿司。由于餐

厅提供的寿司种类繁多,而不同种类的寿司之间相互会有影响:三文鱼寿司和鱿鱼寿司一起吃或许会很棒,但和水

果寿司一起吃就可能会肚子痛。因此,Kiana定义了一个综合美味度\(d_{i,j}\)

(i<j),表示在一次取的寿司中,如果包含

了餐厅提供的从第i份到第j份的所有寿司,吃掉这次取的所有寿司后将获得的额外美味度。由于取寿司需要花费一

些时间,所以我们认为分两次取来的寿司之间相互不会影响。注意在吃一次取的寿司时,不止一个综合美味度会被

累加,比如若Kiana一次取走了第1,2,3种寿司各一份,除了d1,3以外,d1,2,d2,3也会被累加进总美味度中。神奇

的是,Kiana的美食评判标准是有记忆性的,无论是单种寿司的美味度,还是多种寿司组合起来的综合美味度,在

计入Kiana的总美味度时都只会被累加一次。比如,若Kiana某一次取走了第1,2种寿司各一份,另一次取走了第2,3

种寿司各一份,那么这两次取寿司的总美味度为d1,1+d2,2+d3,3+d1,2+d2,3,其中d2,2只会计算一次。奇怪的是,

这家寿司餐厅的收费标准很不同寻常。具体来说,如果Kiana一共吃过了c(c>0)种代号为x的寿司,则她需要为这些

寿司付出mx^2+cx元钱,其中m是餐厅给出的一个常数。现在Kiana想知道,在这家餐厅吃寿司,自己能获得的总美

味度(包括所有吃掉的单种寿司的美味度和所有被累加的综合美味度)减去花费的总钱数的最大值是多少。由于她

不会算,所以希望由你告诉她。

Input###

第一行包含两个正整数n,m,分别表示这家餐厅提供的寿司总数和计算寿司价格中使用的常数。

第二行包含n个正整数,其中第k个数ak表示第k份寿司的代号。

接下来n行,第i行包含n-i+1个整数,其中第j个数di,i+j-1表示吃掉寿司能

获得的相应的美味度,具体含义见问题描述。

N<=100,Ai<=1000

Output###

输出共一行包含一个正整数,表示Kiana能获得的总美味度减去花费的总钱数的最大值。

Sample Input###

3 1

2 3 2

5 -10 15

-10 15

15

Sample Output###

12

【样例1说明】

在这组样例中,餐厅一共提供了3份寿司,它们的代号依次为a1=2,a2=3,a3=2,计算价格时的常数m=1。在保证每次取寿司都能获得新的美味度的前提下,Kiana一共有14种不同的吃寿司方案:

1.Kiana一个寿司也不吃,这样她获得的总美味度和花费的总钱数都是0,两者相减也是0;

2.Kiana只取1次寿司,且只取第1个寿司,即她取寿司的情况为{[1,1]},这样获得的总美味度为5,花费的总钱数

为1-2^2+1*2=6,两者相减为-1;

3.Kiana只取1次寿司,且只取第2个寿司,即她取寿司的情况为{[2,2]},这样获得的总美味度为-10,花费的总钱

数为1-3^2+1*3=12,两者相减为-22;

4.Kiana只取1次寿司,且只取第3个寿司,即她取寿司的情况为{[3,3]},这样获得的总美味度为15,花费的总钱数

为12^2+12=6,两者相减为9;

5.Kiana只取1次寿司,且取第1,2个寿司,即她取寿司的情况为{[1,2]},这样获得的总美味度为5+(-10)+(-10)=-1

5,花费的总钱数为(1-22+1*2)+(1-32+1*3)=18,两者相减为-33;

6.Kiana只取1次寿司,且取第2,3个寿司,即她取寿司的情况为{[2,3]},这样获得的总美味度为(-10)+15+15=20,

花费的总钱数为(1-22+1*2)+(1*32+1*3)=18,两者相减为2;

7.Kiana只取1次寿司,且取第1,2,3个寿司,即她取寿司的情况为{[1,3]},这样获得的总美味度为5+(-10)+15+(-1

0)+15+15=30,花费的总钱数为(122+2*2)+(1*32+13)=20,两者相减为10。

8.Kiana取2次寿司,第一次取第1个寿司,第二次取第2个寿司,即她取寿司的情况为{[1,1],[2,2]},这样获得的

总美味度为5+(-10)=-5,花费的总钱数为(122+1*2)+(1*32+13)=18,两者相减为-23;

9.Kiana取2次寿司,第一次取第1个寿司,第二次取第3个寿司,即她取寿司的情况为{[1,1],[3,3]},这样获得的

总美味度为5+15=20,花费的总钱数为12^2+22=8,两者相减为12;

10.Kiana取2次寿司,第一次取第2个寿司,第二次取第3个寿司,即她取寿司的情况为{[2,2],[3,3]},这样获得的

总美味度为(-10)+15=5,花费的总钱数为(122+1*2)+(1*32+13)=18,两者相减为-13;

11.Kiana取2次寿司,第一次取第1,2个寿司,第二次取第3个寿司,即她取寿司的情况为{[1,2],[3,3]},这样获得

的总美味度为5+(-10)+(-10)+15=0,花费的总钱数为(122+2*2)+(1*32+13)=20,两者相减为-20;

12.Kiana取2次寿司,第一次取第1个寿司,第二次取第2,3个寿司,即她取寿司的情况为{[1,1],[2,3]},这样获得

的总美味度为5+(-10)+15+15=25,花费的总钱数为(1-22+2-2)+(1-32+1-3)=20,两者相减为5;

13.Kiana取2次寿司,第一次取第1,2个寿司,第二次取第2,3个寿司,即她取寿司的情况为{[1,2],[2,3]},这样获

得的总美味度为5+(-10)+15+(-10)+15=15,花费的总钱数为(122+2*2)+(1*32+13)=20,两者相减为-5;

14.Kiana取3次寿司,第一次取第1个寿司,第二次取第2个寿司,第三次取第3个寿司,即她取寿司的情况为{[1,1]

,[2,2],[3,3]},这样获得的总美味度为5+(-10)+15=10,花费的总钱数为(122+2*2)+(1*32+13)=20,两者相减

为-10。

所以Kiana会选择方案9,这时她获得的总美味度减去花费的总钱数的值最大为12。


想法##

我得表示看这种巨长无比的题面心中很不爽。。。

费了好半天看明白题后,发现这是最大权闭合图模型。

当Kiana取走[l,r]的寿司时,她会获得输入数据中对应的一个小直角三角形的美味度之和

如果Kiana想获得d[i,j],那么她也必须获得d[i+1,j]与d[i,j-1]

如果Kiana想获得d[i,i],那么她必须交买第i种寿司的钱x

如果Kiana交了买第i中寿司的钱x,她必须交代号为x的寿司的钱\(mx^2\)

可以把所有d[i,j],每种寿司,每种代号作为点

表示d[i,j]的点,点权为d[i,j]

表示每种寿司的点,点权为-x

表示每种代号的点,点权为\(-mx^2\)

这样我们可以根据这些点之间的依赖关系建一张图了。

然后按照最大权闭合图的套路,点权为正的与S相连,点权为负的与T相连,跑一遍最大流,用所有与S连边的点权之和减去最大流就好了。


代码##

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<map> #define INF 2100000000 using namespace std; const int N = 6000; struct node{
int v,f;
node *next,*rev;
}pool[N*10],*h[N];
int cnt;
void addedge(int u,int v,int f){
node *p=&pool[++cnt],*q=&pool[++cnt];
p->v=v;p->next=h[u];h[u]=p;p->f=f;p->rev=q;
q->v=u;q->next=h[v];h[v]=q;q->f=0;q->rev=p;
} int S,T;
int dis[N],que[N],head,tail;
bool bfs(){
int u,v;
head=tail=0;
for(int i=S;i<=T;i++) dis[i]=-1;
que[tail++]=S; dis[S]=0;
while(head<tail){
u=que[head++];
for(node *p=h[u];p;p=p->next)
if(p->f && dis[v=p->v]==-1){
dis[v]=dis[u]+1;
que[tail++]=v;
}
if(dis[T]!=-1) return true;
}
return false;
}
int find(int u,int f){
int v,s=0,t;
if(u==T) return f; /**/
for(node *p=h[u];p;p=p->next)
if(p->f && s<f && dis[v=p->v]==dis[u]+1){
t=find(v,min(p->f,f-s));
if(t){
p->f-=t;
p->rev->f+=t;
s+=t;
}
}
if(s==0) dis[u]=-1;
return s;
}
int dinic(){
int f=0;
while(bfs()) f+=find(S,INF);
return f;
} int n,m;
int a[105],d[105][105],id[105][105],tot,num;
map<int,int> mp1,mp2; int main()
{
int ans=0;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) {
scanf("%d",&a[i]);
if(!mp1.count(a[i])) mp1[a[i]]=++num,mp2[num]=a[i];
}
for(int i=1;i<=n;i++)
for(int j=i;j<=n;j++)
scanf("%d",&d[i][j]),id[i][j]=++tot; S=0; T=tot+n+num+1;
for(int i=1;i<=n;i++)
for(int j=i;j<=n;j++){
if(d[i][j]>=0) addedge(S,id[i][j],d[i][j]),ans+=d[i][j];
else addedge(id[i][j],T,-d[i][j]);
if(i<n) addedge(id[i][j],id[i+1][j],INF);
if(j>i) addedge(id[i][j],id[i][j-1],INF);
}
for(int i=1;i<=n;i++) {
addedge(tot+i,T,a[i]);
addedge(id[i][i],tot+i,INF);
addedge(tot+i,mp1[a[i]]+tot+n,INF);
}
for(int i=1;i<=num;i++)
addedge(tot+n+i,T,mp2[i]*mp2[i]*m); ans-=dinic();
printf("%d\n",ans); return 0;
}

[bzoj4873] [洛谷P3749] [Shoi2017] 寿司餐厅的更多相关文章

  1. 【BZOJ4873】[Shoi2017]寿司餐厅 最大权闭合图

    [BZOJ4873][Shoi2017]寿司餐厅 Description Kiana最近喜欢到一家非常美味的寿司餐厅用餐.每天晚上,这家餐厅都会按顺序提供n种寿司,第i种寿司有一个代号ai和美味度di ...

  2. 【最大权闭合子图】bzoj4873 [Shoi2017]寿司餐厅

    4873: [Shoi2017]寿司餐厅 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 369  Solved: 256[Submit][Status ...

  3. bzoj4873: [Shoi2017]寿司餐厅(最大权闭合子图)

    4873: [Shoi2017]寿司餐厅 大难题啊啊!!! 题目:传送门 题解:一眼题是网络流,但还是不会OTZ,菜啊... %题解... 最大权闭合子图!!! 好的...开始花式建边: 1.对于每个 ...

  4. BZOJ:4873: [Shoi2017]寿司餐厅

    4873: [Shoi2017]寿司餐厅 首先很开心在膜你赛的时候做了出来. 看到数据范围,看到不能dp,看到贡献去重后计算,咦,流? 那就容易了,转最大权闭合子图,每个区间建一个点,取了就一定要取他 ...

  5. bzoj 4873: [Shoi2017]寿司餐厅 [最小割]

    4873: [Shoi2017]寿司餐厅 题意:略 唯一会做的... 一眼最小割 就是最大权闭合子图呀 \(s\rightarrow d_{positive} \rightarrow -d_{negt ...

  6. BZOJ_4873_[Shoi2017]寿司餐厅_最大权闭合子图

    BZOJ_4873_[Shoi2017]寿司餐厅_最大权闭合子图 题意:http://www.lydsy.com/JudgeOnline/problem.php?id=4873 分析:我们发现分数正负 ...

  7. [LOJ 2146][BZOJ 4873][Shoi2017]寿司餐厅

    [LOJ 2146][BZOJ 4873][Shoi2017]寿司餐厅 题意 比较复杂放LOJ题面好了qaq... Kiana 最近喜欢到一家非常美味的寿司餐厅用餐. 每天晚上,这家餐厅都会按顺序提供 ...

  8. 【洛谷P3749】[六省联考2017]寿司餐厅(网络流)

    洛谷 题意: 给出\(n\)份寿司,现可以选取任意多次连续区间内的寿司,对于区间\([l,r]\),那么贡献为\(\sum_{i=l}^r \sum_{j=i}^rd_{i,j}\)(对于相同的\(d ...

  9. BZOJ4873[Shoi2017]寿司餐厅——最大权闭合子图

    题目描述 Kiana最近喜欢到一家非常美味的寿司餐厅用餐.每天晚上,这家餐厅都会按顺序提供n种寿司,第i种寿司有一个 代号ai和美味度di,i,不同种类的寿司有可能使用相同的代号.每种寿司的份数都是无 ...

随机推荐

  1. MFC入门

    目录 001.MFC_应用程序类型    002.MFC_对话框_静态文本_编辑框  003.MFC_对话框_访问控件_7种方法_A   004.MFC_对话框_访问控件_7种方法_B   005.M ...

  2. .Net Core 3.0 发布为单独的可执行文件

    dotnet publish -r win-x64 -c Release /p:PublishSingleFile=true /p:PublishTrimmed=true

  3. 专业版12.0试用,打开演示账套,提示&ldquo;列名FPlatver无效&rdquo;

    你好,我代表研发来说明一下这个问题:一.产生的原因:由于KIS产品今年陆续都增加了应用平台,前期平台兼容性还有些问题,如果一台电脑有多个不同版本的平台,就会产生冲突,因此报错.1.可能您的机器装了多个 ...

  4. 构建锁与同步组件的基石AQS:深入AQS的实现原理与源码分析

    Java并发包(JUC)中提供了很多并发工具,这其中,很多我们耳熟能详的并发工具,譬如ReentrangLock.Semaphore,它们的实现都用到了一个共同的基类--AbstractQueuedS ...

  5. 010 Ceph RGW对象存储

    一.对象存储 1.1 介绍 通过对象存储,将数据存储为对象,每个对象除了包含数据,还包含数据自身的元数据 对象通过Object ID来检索,无法通过普通文件系统操作来直接访问对象,只能通过API来访问 ...

  6. Linux 学习笔记 4 创建、复制、移动、文件的基本操作

    写在前面 通过上一节的学习,我们基本的了解到在Linux 里面对于设备的挂载.卸载以及设备存在的目录.挂载目录.都有了一个基本的了解 本节主要了解文件.以及目录的相关操作,比如文件.目录的创建.以及删 ...

  7. linux-head、tail、sort、uniq、pstree、ps

    1.head 默认查看文件前10行 head -7 /etc/yum.conf -n num -n7代表查看文件前7行 2.tail 默认查看文件的后10行 tail -7 /etc/yum.conf ...

  8. [UWP]XAML中的响应式布局技术

    响应式布局的概念是一个页面适配多个终端及不同分辨率.在针对特定屏幕宽度优化应用 UI 时,我们将此称为创建响应式设计.WPF设计之初响应式设计的概念并不流行,那时候大部分网页设计师都按着宽度960像素 ...

  9. JavaScript数据类型 - Symbol

    ES5:对象的属性名只能是字符串,当给对象添加新属性时,很容易造成属性名冲突,从而覆盖了原有的属性. ES6:所以ES6中引入了symbol数据类型,他表示独一无二的值,避免了属性名的冲突,此时对象的 ...

  10. 【Python3爬虫】突破反爬之应对前端反调试手段

    一.前言 在我们爬取某些网站的时候,会想要打开 DevTools 查看元素或者抓包分析,但按下 F12 的时候,却出现了下面这一幕:   此时网页暂停加载,自动跳转到 Source 页面并打开了一个 ...