[BZOJ]1324: Exca王者之剑

题目大意:在一个N*M的网格之中，每个格子上都有一定价值的宝石，Amber可以任意选则一个起始点。并且Amber可以瞬间取走当前所在格子中的宝石，Amber每秒可以走一步，在偶数秒时他周围4格中的宝石会消失（包括0秒），问最多能取到宝石的价值和。
这题就是要我们求这个图的最大点权独立集，因为最大点权独立集的点权和=所有点权和-最小点权覆盖的点权和，所以这题转换为求最小点权覆盖的点权和，先把图染色，保证相邻的2点颜色不同（一黑一白），然后从源点向每个白点连一条容量等于白点权值的边，从每个黑点向汇点连一条一条容量等于黑点权值的边，然后从白点向相邻的黑点连一条容量无限大的边，跑一边最大流就可以了。
最小点权覆盖的解法参考：http://blog.youkuaiyun.com/tanmlh/article/details/41979505

#include <cstdio>
#include <cstring>
#include <queue>
#include <vector>
using namespace std;
const int INF=2e9;
int n,m,tot=-1,all,d[10003],cur[10003],ans;
bool vis[10003];
vector<int> o[10003];
queue<int> p;

 struct Edge
{
    int to,cap,flow;
    Edge(int a=0,int b=0,int c=0):to(a),cap(b),flow(c){}
}edge[200002];

inline int min(int a,int b)
{
    return a<b?a:b;
}

bool BFS(void)
{
    memset(vis,0,sizeof(vis));
    vis[0]=1,p.push(0),d[0]=0;
    while(!p.empty())
    {
        int now=p.front();
        p.pop();
        for (int i=0;i<o[now].size();++i)
            if(!vis[edge[o[now][i]].to]&&edge[o[now][i]].flow<edge[o[now][i]].cap)
            {
                vis[edge[o[now][i]].to]=1;
                p.push(edge[o[now][i]].to);
                d[edge[o[now][i]].to]=d[now]+1;
            }
    }
    return vis[all];
}

int DFS(int x,int a)
{
    if(x==all||a==0) return a;
    int flow=0,f;
    for (int &i=cur[x];i<o[x].size();++i)
    {
        Edge&e=edge[o[x][i]];
        if(d[x]+1==d[e.to]&&(f=DFS(e.to,min(a,e.cap-e.flow)))>0)
        {
            e.flow+=f;
            edge[o[x][i]^1].flow-=f;
            a-=f;
            flow+=f;
            if(a==0)
                break;
        }
    }
    return flow;
}

int main(void)
{
    register int i,j,x;
    scanf("%d%d",&n,&m);
    all=n*m+1;
    for (i=1;i<=n;++i)
        for (j=1;j<=m;++j)
        {
            scanf("%d",&x);
            ans+=x;
            if(i+j&1) 
                edge[++tot]=Edge((i-1)*m+j,x,0),o[0].push_back(tot),
                edge[++tot]=Edge(0,0,0),o[(i-1)*m+j].push_back(tot);
            else
            {
                edge[++tot]=Edge(all,x,0),o[(i-1)*m+j].push_back(tot), 
                edge[++tot]=Edge((i-1)*m+j,0,0),o[all].push_back(tot);
                if(i^1)
                    edge[++tot]=Edge((i-1)*m+j,INF,0),o[(i-2)*m+j].push_back(tot), 
                    edge[++tot]=Edge((i-2)*m+j,0,0),o[(i-1)*m+j].push_back(tot);
                if(j^1)
                    edge[++tot]=Edge((i-1)*m+j,INF,0),o[(i-1)*m+j-1].push_back(tot), 
                    edge[++tot]=Edge((i-1)*m+j-1,0,0),o[(i-1)*m+j].push_back(tot);  
                if(i^n)
                    edge[++tot]=Edge((i-1)*m+j,INF,0),o[i*m+j].push_back(tot), 
                    edge[++tot]=Edge(i*m+j,0,0),o[(i-1)*m+j].push_back(tot);
                if(j^m)
                    edge[++tot]=Edge((i-1)*m+j,INF,0),o[(i-1)*m+j+1].push_back(tot), 
                    edge[++tot]=Edge((i-1)*m+j+1,0,0),o[(i-1)*m+j].push_back(tot);
            }
        }
    while(BFS())
    {
        memset(cur,0,sizeof(cur));
        ans-=DFS(0,INF);
    }
    printf("%d",ans);
    return 0;
}