bzoj4873: [Shoi2017]寿司餐厅（最大权闭合子图）

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

　 跑网络流。。。

蒟蒻的代码（丑的一匹，慎点）：

 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<cmath>
 #include<algorithm>
 #define qread(x)x=read();
 using namespace std;
 typedef long long LL;
 inline int read()
 {
     int f=,x=;char ch;
     while(ch<'' || ch>''){if(ch=='-')f=-;ch=getchar();}
     while(ch>='' && ch<=''){x=x*+ch-'';ch=getchar();}
     return f*x;
 }
 struct node
 {
     int x,y,c,next,other;
 }a[];int last[],len;
 int st,ed,n,m;
 void ins(int x,int y,int c)
 {
     int k1,k2;
     len++;k1=len;
     a[len].x=x;a[len].y=y;a[len].c=c;
     a[len].next=last[x];last[x]=len;

     len++;k2=len;
     a[len].x=y;a[len].y=x;a[len].c=;
     a[len].next=last[y];last[y]=len;

     a[k1].other=k2;
     a[k2].other=k1;
 }
 int list[],h[];
 int head,tail;
 bool bfs()
 {
     memset(h,,sizeof(h));h[st]=;
     list[]=st;head=;tail=;
     while(head!=tail)
     {
         int x=list[head];
         for(int k=last[x];k;k=a[k].next)
         {
             int y=a[k].y;
             if(h[y]== && a[k].c>)
             {
                 h[y]=h[x]+;
                 list[tail++]=y;
             }
         }
         head++;
     }
     if(h[ed]>)return true;
     return false;
 }
 int findflow(int x,int flow)
 {
     if(x==ed)return flow;
     int s=,t;
     for(int k=last[x];k;k=a[k].next)
     {
         int y=a[k].y;
         if(h[y]==h[x]+ && a[k].c> && flow>s)
         {
             t=findflow(y,min(a[k].c,flow-s));
             s+=t;
             a[k].c-=t;a[a[k].other].c+=t;
         }
     }
     if(s==)h[x]=;
     return s;
 }
 int e[];
 int d[][];
 int id1[][],id2[];
 bool v[];
 int jb()
 {
     memset(v,false,sizeof(v));
     LL sum=,cnt=;
     for(int i=;i<=n;i++)
         for(int j=i;j<=n;j++)
             id1[i][j]=++cnt;

     for(int i=;i<=n;i++)
         if(v[e[i]]==)
         {
             v[e[i]]=true;
             id2[e[i]]=++cnt;
         }

     ed=cnt+n+;
     memset(v,false,sizeof(v));
     for(int i=;i<=n;i++)
         if(v[e[i]]==)
         {
             v[e[i]]=true;
             ins(id2[e[i]],ed,m*e[i]*e[i]);//种类的花费，和ed连边
         }
     for(int i=;i<=n;i++)
     {
         ins(cnt+i,id2[e[i]],);
         ins(cnt+i,ed,e[i]);
     }
     for(int i=;i<=n;i++)
         for(int j=i;j<=n;j++)
         {
             if(d[i][j]>)
             {
                 sum+=d[i][j];
                 ins(st,id1[i][j],d[i][j]);
                 ins(id1[i][j],cnt+i,);
                 ins(id1[i][j],cnt+j,);
             }
             else if(d[i][j]<)
             {
                 ins(id1[i][j],ed,-d[i][j]);
                 ins(id1[i][j],cnt+i,);
                 ins(id1[i][j],cnt+j,);
             }
             if(i!=j)
             {
                 ins(id1[i][j],id1[i][j-],);
                 ins(id1[i][j],id1[i+][j],);
             }
         }
     return sum;
 }
 int main()
 {
     qread(n);qread(m);
     for(int i=;i<=n;i++)qread(e[i]);
     for(int i=;i<=n;i++)
         for(int j=i;j<=n;j++)
             qread(d[i][j]);
     st=;
     LL sum=jb();
     /*
     for(int i=1;i<=len;i++)
             if(a[i].x==0)
                 printf("%d %d\n",a[i].y,a[i].c);
     printf("\n");
     for(int i=1;i<=len;i++)
             if(a[i].y==ed)
                 printf("%d %d\n",a[i].x,a[i].c);
     */
     int ans=;
     while(bfs())ans+=findflow(st,);
     printf("%lld\n",sum-ans);
     return ;
 }