4873: [Shoi2017]寿司餐厅
大难题啊啊！！！
题目：传送门
题解：一眼题是网络流，但还是不会OTZ，菜啊...
%题解...
最大权闭合子图！！！
好的...开始花式建边：
1、对于每个区间，我们把它看成一个点，按照权值正负连接源点或汇点（最大权闭合子图的套路）
2、对于所有的寿司类型k[i]，各自开一个点，向ed连边，流量为m*k[i]*k[i];
3、对于1~n的每个寿司，向它的所属类型连边，流量为无限大；再向ed连边，流量为e[i]。
4、对于所有的区间i~j，向（i+1，j）和（i，j-1）连边，因为区间具有包含性。
WA了...再%题解...
漏了第五种情况：
5、对于所有区间，向区间内所有的点连边，流量为无限大，表示必须选对应的寿司才能选这个区间QWQ
跑网络流。。。
蒟蒻的代码（丑的一匹，慎点）：
1 #include<cstdio> 2 #include<cstring> 3 #include<cstdlib> 4 #include<cmath> 5 #include<algorithm> 6 #define qread(x)x=read(); 7 using namespace std; 8 typedef long long LL; 9 inline int read() 10 { 11 int f=1,x=0;char ch; 12 while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();} 13 while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();} 14 return f*x; 15 } 16 struct node 17 { 18 int x,y,c,next,other; 19 }a[2110000];int last[1110000],len; 20 int st,ed,n,m; 21 void ins(int x,int y,int c) 22 { 23 int k1,k2; 24 len++;k1=len; 25 a[len].x=x;a[len].y=y;a[len].c=c; 26 a[len].next=last[x];last[x]=len; 27 28 len++;k2=len; 29 a[len].x=y;a[len].y=x;a[len].c=0; 30 a[len].next=last[y];last[y]=len; 31 32 a[k1].other=k2; 33 a[k2].other=k1; 34 } 35 int list[11000],h[11000]; 36 int head,tail; 37 bool bfs() 38 { 39 memset(h,0,sizeof(h));h[st]=1; 40 list[1]=st;head=1;tail=2; 41 while(head!=tail) 42 { 43 int x=list[head]; 44 for(int k=last[x];k;k=a[k].next) 45 { 46 int y=a[k].y; 47 if(h[y]==0 && a[k].c>0) 48 { 49 h[y]=h[x]+1; 50 list[tail++]=y; 51 } 52 } 53 head++; 54 } 55 if(h[ed]>0)return true; 56 return false; 57 } 58 int findflow(int x,int flow) 59 { 60 if(x==ed)return flow; 61 int s=0,t; 62 for(int k=last[x];k;k=a[k].next) 63 { 64 int y=a[k].y; 65 if(h[y]==h[x]+1 && a[k].c>0 && flow>s) 66 { 67 t=findflow(y,min(a[k].c,flow-s)); 68 s+=t; 69 a[k].c-=t;a[a[k].other].c+=t; 70 } 71 } 72 if(s==0)h[x]=0; 73 return s; 74 } 75 int e[110]; 76 int d[110][110]; 77 int id1[110][110],id2[1100]; 78 bool v[1100]; 79 int jb() 80 { 81 memset(v,false,sizeof(v)); 82 LL sum=0,cnt=0; 83 for(int i=1;i<=n;i++) 84 for(int j=i;j<=n;j++) 85 id1[i][j]=++cnt; 86 87 for(int i=1;i<=n;i++) 88 if(v[e[i]]==0) 89 { 90 v[e[i]]=true; 91 id2[e[i]]=++cnt; 92 } 93 94 ed=cnt+n+1; 95 memset(v,false,sizeof(v)); 96 for(int i=1;i<=n;i++) 97 if(v[e[i]]==0) 98 { 99 v[e[i]]=true; 100 ins(id2[e[i]],ed,m*e[i]*e[i]);//种类的花费，和ed连边 101 } 102 for(int i=1;i<=n;i++) 103 { 104 ins(cnt+i,id2[e[i]],999999999); 105 ins(cnt+i,ed,e[i]); 106 } 107 for(int i=1;i<=n;i++) 108 for(int j=i;j<=n;j++) 109 { 110 if(d[i][j]>0) 111 { 112 sum+=d[i][j]; 113 ins(st,id1[i][j],d[i][j]); 114 ins(id1[i][j],cnt+i,999999999); 115 ins(id1[i][j],cnt+j,999999999); 116 } 117 else if(d[i][j]<0) 118 { 119 ins(id1[i][j],ed,-d[i][j]); 120 ins(id1[i][j],cnt+i,999999999); 121 ins(id1[i][j],cnt+j,999999999); 122 } 123 if(i!=j) 124 { 125 ins(id1[i][j],id1[i][j-1],999999999); 126 ins(id1[i][j],id1[i+1][j],999999999); 127 } 128 } 129 return sum; 130 } 131 int main() 132 { 133 qread(n);qread(m); 134 for(int i=1;i<=n;i++)qread(e[i]); 135 for(int i=1;i<=n;i++) 136 for(int j=i;j<=n;j++) 137 qread(d[i][j]); 138 st=0; 139 LL sum=jb(); 140 /* 141 for(int i=1;i<=len;i++) 142 if(a[i].x==0) 143 printf("%d %d ",a[i].y,a[i].c); 144 printf(" "); 145 for(int i=1;i<=len;i++) 146 if(a[i].y==ed) 147 printf("%d %d ",a[i].x,a[i].c); 148 */ 149 int ans=0; 150 while(bfs())ans+=findflow(st,999999999); 151 printf("%lld ",sum-ans); 152 return 0; 153 }