其实一开始看到题意很长就不想做。。
其实一开始想什么spfa什么DP什么的。。
而这题是网络流（好像之前在竞赛室的时候听到队长骂这题stm网络流很难写，感觉还好啊，一中午写完了，2A）
直接网络流显然搞不定，因为每个点有点值啊，费用流感觉也不对，想试试最小割，前天吃饭的时候突然想到：诶这题有负数？那搞搞最大权闭合子图！
显然就正负分开st集合和ed集合，然后我就发现中间保护建的边要到inf，可最大权闭合子图不是sum-最小割吗？昨天走的时候O老师（%%%）才告诉我是正权sum-最小割
然后就这样建，发现数据有环，那就要用topsort判环，但是有个要注意的，假如x被y保护，千万建边要建成y->x，如果建成x->y，假如x在环里面，那y的入度肯定不是0，但是实际上是y保护x，x无法被吃掉对于y是没有影响的。
好感动。。。第二题最大权闭合子图居然看出来了而且2A。。。
#include<cstdio> #include<iostream> #include<cstring> using namespace std; int n,m,st,ed,sum; int w[610]; struct enode { int x,y,next; }e[410000];int elen,elast[610]; void eins(int x,int y) { elen++; e[elen].x=x;e[elen].y=y; e[elen].next=elast[x];elast[x]=elen; } int ru[610],top,sta[610]; void topsort() { top=0; for(int i=1;i<=n*m+2;i++) if(ru[i]==0)sta[++top]=i; while(top!=0) { int x=sta[top];top--; for(int k=elast[x];k;k=e[k].next) { int y=e[k].y; ru[y]--; if(ru[y]==0)sta[++top]=y; } } sum=0; for(int i=1;i<=n*m;i++) if(ru[i]==0&&w[i]>=0)sum+=w[i]; } struct node { int x,y,c,next,other; }a[810000];int len,last[610]; 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=0; a[len].next=last[y];last[y]=len; a[k1].other=k2; a[k2].other=k1; } void composition() { len=0;memset(last,0,sizeof(last)); for(int i=1;i<=elen;i++) { int x=e[i].x,y=e[i].y; if(ru[x]==0&&ru[y]==0)ins(x,y,999999999); } for(int i=1;i<=n*m;i++) if(ru[i]==0) if(w[i]>=0)ins(i,ed,w[i]); else ins(st,i,-w[i]); } int h[610],list[610]; bool bt_h() { memset(h,0,sizeof(h));h[st]=1; int head=1,tail=2;list[1]=st; 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]==0&&a[k].c>0) { h[y]=h[x]+1; list[tail]=y; tail++; } } head++; } if(h[ed]==0)return false; return true; } int findflow(int x,int f) { if(x==ed)return f; int s=0; for(int k=last[x];k;k=a[k].next) { int y=a[k].y; if(h[y]==h[x]+1&&a[k].c>0&&s<f) { int t=findflow(y,min(a[k].c,f-s)); s+=t;a[k].c-=t;a[a[k].other].c+=t; } } if(s==0)h[x]=0; return s; } int point(int x,int y){return x*m+y+1;} int main() { int T,x,y; scanf("%d%d",&n,&m); st=n*m+1;ed=n*m+2; elen=0;memset(elast,0,sizeof(elast)); memset(ru,0,sizeof(ru)); for(int i=0;i<n;i++) { for(int j=0;j<m;j++) { int p=point(i,j); scanf("%d",&w[p]); //我所能保护的 int pro; if(j!=0) { pro=point(i,j-1); eins(p,pro), ru[pro]++; } scanf("%d",&T); while(T--) { scanf("%d%d",&x,&y); pro=point(x,y); eins(p,pro), ru[pro]++; } } } topsort(); composition(); int ans=0; while(bt_h()==true) { ans+=findflow(st,999999999); } printf("%d ",sum-ans); return 0; }