題目描述
輸入
輸出
僅包含一個整數,表示可以獲得的最大能源收入。注意,你也可以選擇不進行任何攻擊,這樣能源收入爲0。
樣例輸入
3 2
10 0
20 0
-10 0
-5 1 0 0
100 1 2 1
100 0
樣例輸出
25
。。一看就是一個網絡流,首先有兩種植物,一種有收益,一種需花費。
我們可以嘗試着想到最小割,先統計出總收益tot(所有收益之和),在想辦法構造出一個圖,
使得它的每一個割對應:不要的收益+所花的花費 。那麼這時,答案就等於 tot-最小割。
至於怎麼構圖,我們可以吧最小割割開的圖看做S->T集(即和源或匯點分在一個集合中的點集)
S集中的點一定是被喫掉的植物,相反,T中就是沒有被喫的。
於是,我們考慮將收益和花費的植物分別將S,T連邊,有因爲殭屍是從左向右喫的,要想喫掉一棵植物必須喫掉保護他的那棵植物,對於這樣的拓撲關係,很明顯要先將先被喫的那棵植物劃分到S集,也就是將他和T的連邊斷掉,於是對於每一點,我們向在拓撲關係中的他上一層的節點連一條爲INF的邊。
不過這裏要注意,當拓撲關係中存在環時,在流網絡中就形成了一個INF的環路,最小割割掉的一定是環路與S,T的所有邊,那麼這就錯了。
於是我們必須預處理拓撲關係圖中的所有環路,將所有不可能到達的點刪去,再構圖跑網絡流,這裏JeremyGuo提供了一個很方便的解決方案:我們對於原圖中的每一個強聯通分量都從中找出一個點,將其與T連一條INF邊,那麼這個強連通分量就一定屬於T集合,於是就只算了不需要的收益。至於找點向T連邊,這個有很多搞法,而且又不是一個分量只能連一條,所以亂搞搞就行了。
/**************************************************************
Problem: 1565
User: szpszp
Language: C++
Result: Accepted
Time:188 ms
Memory:29440 kb
****************************************************************/
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
using namespace std;
#define MAXN 20
#define MAXM 30
#define INF 0x3f3f3f3f
typedef long long int LL;
template<class T>
void Read(T &x){
x=0;char c=getchar();bool flag=0;
while(c<'0'||'9'<c){if(c=='-')flag=1;c=getchar();}
while('0'<=c&&c<='9'){x=x*10+c-'0';c=getchar();}
if(flag)x=-x;
}
namespace ISAP{
const int MAXP = MAXN*MAXM;
const int MAXE = MAXN*MAXM*MAXN*MAXM*2;
struct node{
int v,c;
node *nxt,*bck;
}*adj[MAXP+10],Edges[MAXE*2+100],*New;
int d[MAXP+10],vd[MAXP+10];
int S,T,P;
void addedges(int u,int v,int c){
node *p=++New;
p->v=v;
p->c=c;
p->nxt=adj[u];
p->bck=New+1;
adj[u]=p;
p=++New;
p->v=u;
p->c=0;
p->nxt=adj[v];
p->bck=New-1;
adj[v]=p;
}
void init(int s,int t,int p){
S=s,T=t,P=p;
memset(adj,0,sizeof(adj));
New=Edges;
}
int aug(int x,int augco){
if(x==T)return augco;
int Dmin=P-1,augc=augco,delta;
for(node *p=adj[x];p!=NULL;p=p->nxt)
if(p->c){
if(d[x]==d[p->v]+1){
delta=min(augc,p->c);
delta=aug(p->v,delta);
p->c-=delta;
p->bck->c+=delta;
augc-=delta;
if(d[S]>=P)
return augco-augc;
if(!augc)
break;
}
Dmin=min(Dmin,d[p->v]);
}
if(augco==augc){
if(!(--vd[d[x]]))
d[S]=P;
++vd[d[x]=Dmin+1];
}
return augco-augc;
}
int isap(){
memset(d,0,sizeof(d));
memset(vd,0,sizeof(vd));
vd[0]=P;
int flow=0;
while(d[S]<P)
flow+=aug(S,INF);
return flow;
}
}
int w[MAXN+10][MAXM+10];
int n,m;
int id[MAXN+10][MAXM+10],tot;
void makeid(int n,int m){
tot=0;
for(int i=0;i<n;++i)
for(int j=0;j<m;++j)
id[i][j]=++tot;
}
struct node{
int v;
node *nxt;
}*adj[MAXN*MAXM+10],Edges[MAXN*MAXM*MAXN*MAXM*2+100],*New=Edges;
void addedge(int u,int v){
node *p=++New;
p->v=v;
p->nxt=adj[u];
adj[u]=p;
}
bool vis[MAXN*MAXM+10],tag[MAXN*MAXM+10];
vector<int>point;
void dfs(int x){
vis[x]=tag[x]=1;
for(node *p=adj[x];p!=NULL;p=p->nxt){
if(!vis[p->v])dfs(p->v);
else if(tag[p->v])point.push_back(p->v);
}
tag[x]=0;
}
int main(){
Read(n),Read(m);
makeid(n,m);
int t,x,y;
for(int i=0;i<n;++i)
for(int j=0;j<m;++j){
Read(w[i][j]);
Read(t);
for(int k=0;k<t;++k){
Read(x),Read(y);
addedge(id[i][j],id[x][y]);
}
}
for(int i=0;i<n;++i)
for(int j=1;j<m;++j)addedge(id[i][j],id[i][j-1]);
for(int i=1;i<=tot;++i)
if(!vis[i])dfs(i);
ISAP::init(tot+1,tot+2,tot+2);
for(int i=0;i<n;++i)
for(int j=0;j<m;++j){
if(w[i][j]>=0){
ISAP::addedges(ISAP::S,id[i][j],w[i][j]);
//ISAP::addedges(id[i][j],ISAP::T,0);
}
else{
//ISAP::addedges(ISAP::S,id[i][j],0);
ISAP::addedges(id[i][j],ISAP::T,-w[i][j]);
}
for(node *p=adj[id[i][j]];p!=NULL;p=p->nxt)
ISAP::addedges(p->v,id[i][j],INF);
}
for(int i=0;i<point.size();++i)
ISAP::addedges(point[i],ISAP::T,INF);
int ans=0;
for(int i=0;i<n;++i)
for(int j=0;j<m;++j)
if(w[i][j]>0)ans+=w[i][j];
printf("%d\n",ans-ISAP::isap());
}