BZOJ4930: 棋盘

BZOJ4930: 棋盘

https://lydsy.com/JudgeOnline/problem.php?id=4930

分析:

• 基本上就是游戏那道题加上费用流了，所以没啥好说的。
• 记得两边都是拆边。

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <cmath>
using namespace std;
#define inf 0x3f3f3f3f
#define N 2550
#define M 500050
typedef double f2;
namespace EK {
    const int S=N-1,T=N-2;
    int head[N],to[M],nxt[M],flow[M],val[M],cnt=1;
    int dis[N],Q[N],path[N],vis[N];
    inline void add(int u,int v,int f,int c) {
        to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt; flow[cnt]=f; val[cnt]=c;
        to[++cnt]=u; nxt[cnt]=head[v]; head[v]=cnt; flow[cnt]=0; val[cnt]=-c;
    }
    bool spfa() {
        memset(path,0,sizeof(path));
        memset(dis,0x3f,sizeof(dis));
        int l=0,r=0;
        Q[r++]=S; dis[S]=0;
        while(l!=r) {
            int x=Q[l++]; if(l==S) l=0;
            int i; vis[x]=0;
            for(i=head[x];i;i=nxt[i]) if(dis[to[i]]>dis[x]+val[i]&&flow[i]) {
                dis[to[i]]=dis[x]+val[i]; path[to[i]]=i^1;
                if(!vis[to[i]]) {
                    vis[to[i]]=1; Q[r++]=to[i]; if(r==S) r=0;
                }
            }
        }
        return path[T]!=0;
    }
    pair<int,int> ek() {
        int minc=0,maxf=0;
        while(spfa()) {
            int nf=inf;
            int i;
            for(i=T;i!=S;i=to[path[i]]) {
                nf=min(nf,flow[path[i]^1]);
            }maxf+=nf;
            for(i=T;i!=S;i=to[path[i]]) {
                flow[path[i]^1]-=nf;
                flow[path[i]]+=nf;
                minc+=nf*val[path[i]^1];
            }
        }
        return make_pair(minc,maxf);
    }
}
int n,bl[2][55][55],m,ans[2550];
char mp[55][55];
int main() {
    using namespace EK;
    scanf("%d",&n);
    int i,j,k;
    for(i=1;i<=n;i++) scanf("%s",mp[i]+1);
    for(i=1;i<=n;i++) {
        for(j=1;j<=n;j++) {
            if(mp[i][j]=='#') {continue;}
            int st=j;
            for(;j<n&&mp[i][j+1]=='.';j++) ;
            m++;
            for(k=st;k<=j;k++) bl[0][i][k]=m;
            for(k=1;k<=n;k++) add(S,m,1,k-1);
        }
    }
    for(j=1;j<=n;j++) {
        for(i=1;i<=n;i++) {
            if(mp[i][j]=='#') {continue;}
            int st=i;
            for(;i<n&&mp[i+1][j]=='.';i++) ;
            m++;
            for(k=st;k<=i;k++) bl[1][k][j]=m;
            for(k=1;k<=n;k++) {
                add(m,T,1,k-1);
            }
        }
    }
    for(i=1;i<=n;i++) for(j=1;j<=n;j++) if(mp[i][j]=='.') {
        add(bl[0][i][j],bl[1][i][j],1,0);
    }
    int now=1;
    int minc=0;
    while(spfa()) {
        int nf=inf;
        int i;
        for(i=T;i!=S;i=to[path[i]]) {
            nf=min(nf,flow[path[i]^1]);
        }
        for(i=T;i!=S;i=to[path[i]]) {
            flow[path[i]^1]-=nf;
            flow[path[i]]+=nf;
            minc+=nf*val[path[i]^1];
        }
        ans[now]=minc; now+=1;
    }
    int q;
    scanf("%d",&q);
    while(q--) {
        scanf("%d",&k);
        printf("%d\n",ans[k]);
    }
}

转载于:https://www.cnblogs.com/suika/p/10205913.html