大佬们是怎么一眼看出这是一个最大权闭合子图的……大佬好强->这里
1.把所有区间$(i,j)$看成一个点，如果权值大于0，则从$S$向他连边，容量为权值，否则从它向$T$连边，容量为权值的相反数
2.对于区间$(i,j)$，向所有的寿司$i$到$j$连边，表示选这个区间这些寿司必须选
3.对每一个寿司类型$w[i]$，为他们各自开一个点，向$T$连边，容量为$m*w[i]*w[i]$
4.对每一个寿司向他们所属的类型$w[i]$连边，容量为$inf$，向$T$连边，容量为$w[i]$
5.对于所有区间$(i,j)$，向区间$(i-1,j),(i,j-1)$连边，容量$inf$，表示选了大的必须选小的
然后用总收益减去最小割就行了
1 //minamoto 2 #include<iostream> 3 #include<cstdio> 4 #include<cstring> 5 #include<queue> 6 #define inf 0x3f3f3f3f 7 #define ll long long 8 using namespace std; 9 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++) 10 char buf[1<<21],*p1=buf,*p2=buf; 11 inline int read(){ 12 #define num ch-'0' 13 char ch;bool flag=0;int res; 14 while(!isdigit(ch=getc())) 15 (ch=='-')&&(flag=true); 16 for(res=num;isdigit(ch=getc());res=res*10+num); 17 (flag)&&(res=-res); 18 #undef num 19 return res; 20 } 21 const int N=50005,M=2000005; 22 int head[N],Next[M],ver[M],edge[M],tot=1; 23 int dep[N],cur[N],S,T; 24 queue<int> q; 25 inline void add(int u,int v,int e){ 26 ver[++tot]=v,Next[tot]=head[u],head[u]=tot,edge[tot]=e; 27 ver[++tot]=u,Next[tot]=head[v],head[v]=tot,edge[tot]=0; 28 } 29 bool bfs(){ 30 while(!q.empty()) q.pop(); 31 for(int i=S;i<=T;++i) cur[i]=head[i]; 32 memset(dep,-1,sizeof(dep)); 33 q.push(S),dep[S]=0; 34 while(!q.empty()){ 35 int u=q.front();q.pop(); 36 for(int i=head[u];i;i=Next[i]){ 37 int v=ver[i]; 38 if(dep[v]<0&&edge[i]){ 39 dep[v]=dep[u]+1,q.push(v); 40 if(v==T) return true; 41 } 42 } 43 } 44 return false; 45 } 46 int dfs(int u,int limit){ 47 if(u==T||!limit) return limit; 48 int flow=0,f; 49 for(int i=cur[u];i;i=cur[u]=Next[i]){ 50 int v=ver[i]; 51 if(dep[v]==dep[u]+1&&(f=dfs(v,min(limit,edge[i])))){ 52 flow+=f,limit-=f; 53 edge[i]-=f,edge[i^1]+=f; 54 if(!limit) break; 55 } 56 } 57 if(!flow) dep[u]=-1; 58 return flow; 59 } 60 int dinic(){ 61 int flow=0; 62 while(bfs()) flow+=dfs(S,inf); 63 return flow; 64 } 65 int n,m,a[105],mp[105][105],id[105][105],cnt=0; 66 int idw[1005];bool vis[1005];ll sum=0; 67 void build(){ 68 S=0; 69 for(int i=1;i<=n;++i) 70 for(int j=1;j<=n;++j) 71 id[i][j]=++cnt; 72 for(int i=1;i<=n;++i) 73 if(!vis[a[i]]) vis[a[i]]=1,idw[a[i]]=++cnt; 74 T=cnt+n+1; 75 memset(vis,0,sizeof(vis)); 76 for(int i=1;i<=n;++i) 77 if(!vis[a[i]]) vis[a[i]]=1,add(idw[a[i]],T,m*a[i]*a[i]); 78 for(int i=1;i<=n;++i) 79 add(i+cnt,idw[a[i]],inf),add(i+cnt,T,a[i]); 80 for(int i=1;i<=n;++i) 81 for(int j=i;j<=n;++j){ 82 if(mp[i][j]>0){ 83 sum+=mp[i][j]; 84 add(S,id[i][j],mp[i][j]); 85 add(id[i][j],i+cnt,inf); 86 add(id[i][j],j+cnt,inf); 87 } 88 else if(mp[i][j]<0){ 89 add(id[i][j],T,-mp[i][j]); 90 add(id[i][j],i+cnt,inf); 91 add(id[i][j],j+cnt,inf); 92 } 93 if(i!=j){ 94 add(id[i][j],id[i+1][j],inf); 95 add(id[i][j],id[i][j-1],inf); 96 } 97 } 98 } 99 int main(){ 100 //freopen("testdata.in","r",stdin); 101 n=read(),m=read(); 102 for(int i=1;i<=n;++i) a[i]=read(); 103 for(int i=1;i<=n;++i) 104 for(int j=i;j<=n;++j) 105 mp[i][j]=read(); 106 build(); 107 printf("%lld ",sum-dinic()); 108 return 0; 109 }