BZOJ3232: 圈地游戏

最新推荐文章于 2019-11-30 10:17:03 发布

MirrorGray

最新推荐文章于 2019-11-30 10:17:03 发布

阅读量1.1k

点赞数

分类专栏：二分 spfa OIER的世界文章标签：二分 spfa

本文链接：https://blog.csdn.net/mirrorgray/article/details/51103138

版权

OIER的世界同时被 3 个专栏收录

31 篇文章 3 订阅

订阅专栏

二分

1 篇文章 0 订阅

订阅专栏

spfa

1 篇文章 0 订阅

订阅专栏

本文介绍了如何通过二分法结合SPFA算法解决圈地游戏中的负环判断问题。通过建立特定的图结构，利用前缀和计算边权，确保能够反映出C乘以中间值mid与总费用V的差值是否小于0。这种方法巧妙地避免了对特定路径的限制，简化了问题的求解过程。

摘要由CSDN通过智能技术生成

膜拜下whx大爷…好神的做法…
首先比值类型的最大值一般要二分答案转判定。现在要你找到一个环使得C*mid-V<0。这十分类似spfa判负环…
所以我们考虑建图找负环，每个网格线的交点向周围四个方向建边，边权怎么办呢？将格子上的数在列方向做前缀和，a[i][j]表示第j列，前i行格子的权值和，b[i][j]是横向的格线上的费用，那么横向格线的边从左向右方向的权值是b[i][j]*mid+a[i][j]，从右往左的权值是b[i][j]*mid-a[i][j]，表示对于每一列来说，加上第一次经过的那条边向上所有格子的权值，减去第二次经过的那条边向上所有格子的权值，加上第三次，减掉第四次…以此类推，这样如果原图中存在一条负环，就说明存在一条回路，使得C*mid-V<0。
为什么不会出现下图这种情况呢？显然下图这种情况，只选择一块不劣…
这里写图片描述
（论看题的重要性…我还以为必须是左上角到右下角的一条回路，想了好长时间不会…所以我们也可以很简单建图网络流判定的，但是个人感觉不如spfa优美巧妙…）

#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
//by:MirrorGray
using namespace std;
const int Y=53,N=2711,M=2111111;
const double eps=1e-6;
int n,m,cnt,pos[Y][Y];
int a[Y][Y],b[Y][Y],c[Y][Y];
int inq[N],vis[N];
double dis[N],len[M];
int tot=-1,head[N],ver[M],nxt[M];
int mx[]={0,1},my[]={1,0};
queue <int> q;

void add(int x,int y){
    nxt[++tot]=head[x];head[x]=tot;ver[tot]=y;
    nxt[++tot]=head[y];head[y]=tot;ver[tot]=x;
}

bool spfa(){
    while(!q.empty())q.pop();
    for(int i=1;i<=cnt;i++){
        dis[i]=0,inq[i]=true,vis[i]=0;
        q.push(i);
    }
    while(!q.empty()){
        int x=q.front();q.pop();
        inq[x]=false;
        for(int i=head[x];~i;i=nxt[i])if(dis[ver[i]]>dis[x]+len[i]+eps){
            dis[ver[i]]=dis[x]+len[i];
            vis[ver[i]]=vis[x]+1;
            if(vis[ver[i]]>cnt)return true;//attention please
            if(!inq[ver[i]])inq[ver[i]]=true,q.push(ver[i]);
        }
    }
    return false;
}

bool judge(double mid){
    tot=-1;
    for(int i=0;i<=n;i++)for(int j=0;j<=m;j++)for(int k=0;k<2;k++)
    if(i+mx[k]<=n&&j+my[k]<=m){
        if(mx[k])len[++tot]=c[i][j]*mid,len[++tot]=c[i][j]*mid;
        else len[++tot]=b[i][j]*mid+a[i][j],len[++tot]=b[i][j]*mid-a[i][j];
    }
    return spfa();
}

int main(){
    memset(head,-1,sizeof(head));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)for(int j=0;j<m;j++)scanf("%d",&a[i][j]),a[i][j]+=a[i-1][j];
    for(int i=0;i<=n;i++)for(int j=0;j<m;j++)scanf("%d",&b[i][j]);
    for(int i=0;i<n;i++)for(int j=0;j<=m;j++)scanf("%d",&c[i][j]);
    for(int i=0;i<=n;i++)for(int j=0;j<=m;j++)pos[i][j]=++cnt;
    for(int i=0;i<=n;i++)for(int j=0;j<=m;j++)for(int k=0;k<2;k++)
    if(i+mx[k]<=n&&j+my[k]<=m)add(pos[i][j],pos[i+mx[k]][j+my[k]]);
//  if(judge(2))puts("bingo");else puts("fuck");
//  if(judge(1))puts("bingo");else puts("fuck");
//  if(judge(0))puts("bingo");else puts("fuck");
    double l=0,r=n*m*100;
    for(int i=1;i<=40;i++){
        double mid=(l+r)/2;
        if(judge(mid))l=mid;else r=mid;
    }
    printf("%.3f\n",(l+r)/2);
    return 0;
}