我冒着明天困屎= = 的危险写了这道题和题解
【题解】
最大权闭合子图模型,注意同一排里面的植物可以保护外面的
若A能保护B,就连边 A->B,容量无穷大,将s与所有正权点相连,所有负权点与t相连,容量都为点权的绝对值
还要注意若A,B互相保护,即图上有环,或环能保护一个点,则它们都不可能取!预处理时把它们都删掉。这个预处理只需记录入度,拓扑排序即可。再注意,千万别把s和t也考虑进去了!
然后只考虑有效点,用ISAP做一遍就行了
【代码】
#include<stdio.h>
#include<stdlib.h>
#define INF 10000000
int ord[250][350],val[800],u[800000],v[800000],cap[800000],flow[800000],first[800],next[800000];
int f[800],q[800000],d[800],hash[800],p[800],cur[800],num[800000];
int e=0;
int min(int a,int b)
{
if(a<b) return a;
return b;
}
int fan(int x)
{
if(x&1) return x+1;
return x-1;
}
void tj(int x,int y,int C,int flag)
{
u[++e]=x;
v[e]=y;
cap[e]=C;
next[e]=first[x];
first[x]=e;
if(flag==1&&cap[e]==0) f[y]++;
}
int main()
{
int n,m,s,t,head=0,tail=0,i,j,tot=0,x,y,ans=0,a=0;
scanf("%d%d",&n,&m);
for(i=0;i<n;i++)
for(j=0;j<m;j++)
ord[i][j]=++tot;
s=++tot;
t=++tot;
for(i=0;i<n;i++)
for(j=0;j<m;j++)
{
scanf("%d%d",&val[ord[i][j]],&tot);
if(val[ord[i][j]]>0)
{
tj(s,ord[i][j],val[ord[i][j]],0);
tj(ord[i][j],s,0,0);
}
if(val[ord[i][j]]<0)
{
tj(ord[i][j],t,-val[ord[i][j]],0);
tj(t,ord[i][j],0,0);
}
for(;tot>0;tot--)
{
scanf("%d%d",&x,&y);
tj(ord[x][y],ord[i][j],INF,1);
tj(ord[i][j],ord[x][y],0,1);
}
}
for(i=0;i<n;i++)
for(j=0;j<m-1;j++)
{
tj(ord[i][j],ord[i][j+1],INF,1);
tj(ord[i][j+1],ord[i][j],0,1);
}
for(i=1;i<s;i++)//拓扑排序判环
if(f[i]==0) q[tail++]=i;
while(head<tail)
{
for(i=first[q[head]];i!=0;i=next[i])
if(cap[i]==0&&v[i]<s)
{
f[v[i]]--;
if(f[v[i]]==0) q[tail++]=v[i];
}
head++;
}
for(i=1;i<=t;i++)
d[i]=t;
d[t]=0;
head=0;
tail=1;
q[0]=t;
hash[t]=1;
while(head<tail)
{
for(i=first[q[head]];i!=0;i=next[i])
if(f[v[i]]==0&&d[v[i]]>d[u[i]]+1&&cap[fan(i)]>0)
{
d[v[i]]=d[u[i]]+1;
cur[v[i]]=fan(i);
if(hash[v[i]]==0)
{
q[tail++]=v[i];
hash[v[i]]=1;
}
}
hash[q[head++]]=0;
}
for(i=1;i<=t;i++)
if(f[i]==0)
{
num[d[i]]++;
if(val[i]>0) ans+=val[i];
}
x=s;
while(d[s]<t)
{
if(x==t)
{
a=INF;
for(i=t;i!=s;i=u[p[i]])
a=min(a,cap[p[i]]-flow[p[i]]);
for(i=t;i!=s;i=u[p[i]])
{
flow[p[i]]+=a;
flow[fan(p[i])]-=a;
}
ans-=a;
x=s;
}
for(i=cur[x];i!=0;i=next[i])
if(d[x]==d[v[i]]+1&&cap[i]>flow[i])
{
cur[x]=i;
p[v[i]]=i;
x=v[i];
break;
}
if(i==0)
{
a=t-1;
for(i=first[x];i!=0;i=next[i])
if(cap[i]>flow[i]) a=min(a,d[v[i]]);
if(--num[d[x]]==0) break;
d[x]=a+1;
num[d[x]]++;
cur[x]=first[x];
if(x!=s) x=u[p[x]];
}
}
printf("%d",ans);
return 0;
}