[六省联考2017]寿司餐厅

Description

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

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

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

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

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

果寿司一起吃就可能会肚子痛。因此，Kiana定义了一个综合美味度di,j(i 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+12=6，两者相减为-1；
3.Kiana只取1次寿司，且只取第2个寿司，即她取寿司的情况为{[2,2]}，这样获得的总美味度为-10，花费的总钱
数为1-3^2+13=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-2^2+12)+(1-3^2+13)=18，两者相减为-33；
6.Kiana只取1次寿司，且取第2,3个寿司，即她取寿司的情况为{[2,3]}，这样获得的总美味度为(-10)+15+15=20，
花费的总钱数为(1-2^2+12)+(13^2+13)=18，两者相减为2；
7.Kiana只取1次寿司，且取第1,2,3个寿司，即她取寿司的情况为{[1,3]}，这样获得的总美味度为5+(-10)+15+(-1
0)+15+15=30，花费的总钱数为(12^2+22)+(13^2+13)=20，两者相减为10。
8.Kiana取2次寿司，第一次取第1个寿司，第二次取第2个寿司，即她取寿司的情况为{[1,1],[2,2]}，这样获得的
总美味度为5+(-10)=-5，花费的总钱数为(12^2+12)+(13^2+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，花费的总钱数为(12^2+12)+(13^2+13)=18，两者相减为-13；
11.Kiana取2次寿司，第一次取第1,2个寿司，第二次取第3个寿司，即她取寿司的情况为{[1,2],[3,3]}，这样获得
的总美味度为5+(-10)+(-10)+15=0，花费的总钱数为(12^2+22)+(13^2+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，花费的总钱数为(12^2+22)+(13^2+13)=20，两者相减为-5；
14.Kiana取3次寿司，第一次取第1个寿司，第二次取第2个寿司，第三次取第3个寿司，即她取寿司的情况为{[1,1]
,[2,2],[3,3]}，这样获得的总美味度为5+(-10)+15=10，花费的总钱数为(12^2+22)+(13^2+1*3)=20，两者相减
为-10。
所以Kiana会选择方案9，这时她获得的总美味度减去花费的总钱数的值最大为12。

Solution

最近写题的时间确实少了点，去搞网络流的总结了。。

说一下这个题。明显题目这么长。不仅长而且不好懂。。。看着样例都看不出来输入的啥。。。

不过经过随便考虑之后就可以发现这题明显是一个最大权闭合子图的模型。问怎么看出来的？不很明显的选了大的区间就一定要选小区间么？然后自然连边。然后正权区间向S连边，负权区间向T连边。好了，接着解决花费的事就可以了。可以发现的是把那个式子随便拆拆就可以拆成这样一个样子\(mx^2和cx\)。把所有的寿司种类都开一个点，然后把每个寿司连向它所属的种类，然后对于前面那个，把种类连向T，对于后面的，把每个寿司再向T连边。好了搞定。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#define re register
#define inf 400000000
#define MAXN 6055
#define MAXM 100001
#define get(a,b) 1000+(a-1)*(b-a+2)+(a-1)*(a-2)/2
using namespace std;
int n,s,q,dis[2000011],t,l,cur[200051],m;
struct po
{
    int nxt,to,w;
}edge[MAXM];
int head[MAXN],dep[MAXN],num=-1,de[105][105],id[105][105],vis[5000],a[1001],cnt,tot;
inline char gc(){
    static char BUFF[1000000],*S=BUFF,*T=BUFF;
    return S==T&&(T=(S=BUFF)+fread(BUFF,1,1000000,stdin),S==T)?EOF:*S++;
}
inline int read()
{
    int x=0,c=1;
    char ch=' ';
    while((ch>'9'||ch<'0')&&ch!='-')ch=gc();
    while(ch=='-')c*=-1,ch=gc();
    while(ch<='9'&&ch>='0')x=x*10+ch-'0',ch=gc();
    return x*c;
}
inline void add_edge(int from,int to,int w)
{
    edge[++num].nxt=head[from];
    edge[num].to=to;
    edge[num].w=w;
    head[from]=num;
}
inline void add(int from,int to,int w)
{
    add_edge(from,to,w);
    add_edge(to,from,0);
}
inline bool bfs()
{
    memset(dep,0,sizeof(dep));
    queue<int> q;
    while(!q.empty())
    q.pop();
    q.push(s);
    dep[s]=1;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(re int i=head[u];i!=-1;i=edge[i].nxt)
        {
            int v=edge[i].to;
            if(dep[v]==0&&edge[i].w>0)
            {
                dep[v]=dep[u]+1;
                if(v==t)
                return 1;
                q.push(v);
            }
        }
    }
    return 0;
}
inline int dfs(int u,int dis)
{
    if(u==t)
    return dis;
    int diss=0;
    for(re int& i=cur[u];i!=-1;i=edge[i].nxt)
    {
        int v=edge[i].to;
        if(edge[i].w!=0&&dep[v]==dep[u]+1)
        {
            int check=dfs(v,min(dis,edge[i].w));
            if(check>0)
            {
                dis-=check;
                diss+=check;
                edge[i].w-=check;
                edge[i^1].w+=check;
                if(dis==0) break;
            }
        }
    }
    return diss;
}
inline int dinic()
{
    int ans=0;
    while(bfs())
    {
        for(re int i=0;i<=t;i++)
        cur[i]=head[i];
        while(int d=dfs(s,inf))
        ans+=d;
    }
    return ans;
}
int main() 
{
    memset(head,-1,sizeof(head));
    n=read();m=read();
    s=0,t=1000+n*(n+1)/2+1;
    for(re int i=1;i<=n;i++){
        a[i]=read();
    }
    cnt=1000;
    for(re int i=1;i<=n;i++)
        for(re int j=i;j<=n;j++){
            de[i][j]=read();
            id[i][j]=++cnt;
        }
    for(re int i=1;i<=n;i++){
        if(!vis[a[i]]&&m){
            add(a[i],t,a[i]*a[i]*m);
            vis[a[i]]=1;
        }
        add(id[i][i],a[i],inf);
        add(id[i][i],t,a[i]);
    }
    for(re int i=1;i<=n;i++)
        for(re int j=i;j<=n;j++){
            if(de[i][j]>0) add(s,id[i][j],de[i][j]),tot+=de[i][j];
            else add(id[i][j],t,-de[i][j]);
            if(i!=j){
                add(id[i][j],id[i][j-1],inf);
                add(id[i][j],id[i+1][j],inf);
            }
        }
    cout<<tot-dinic();
    return 0;
}//貌似直接暴力开就能过。。。