Fishes题解

最新推荐文章于 2020-07-07 18:50:20 发布

转载最新推荐文章于 2020-07-07 18:50:20 发布 · 148 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/ljk123-de-bo-ke/p/11318974.html

Fishes题解

我们考虑每一条鱼的贡献，

显然是能被选中的矩形的个数/矩形总数

fish

可知，我们普通的鱼能被选中的矩形个数为r*r[从上至下r,从左至右r],

但是有时会受到矩形大小的限制，

此时x-r+1<1或者x+r-1>n或者y-r+1<1或者y+r-1>m，

能选矩形的下界为min(n-r+1,x),上界为max(1,n-r+1),上下长度为min(n-r+1,x)-max(1,x-r+1)+1=min(n-r+1,x)-max(0,x-r),

同理，左右长度为min(m-r+1,y)-max(0,y-r)，个数即为（min(n-r+1,x)-max(0,x-r)）×（min(m-r+1,y)-max(0,y-r)）。

矩形总数为：(n-r+1)*(n-r+1);

选出k条鱼的方案使期望最大，即是选出前k大期望的正方形格子，越在中间，肯定越不受矩形边界的限制，我们不妨从最中间起，向四周拓展。

代码：

#include<bits/stdc++.h>
using namespace std;
const int N=100006,f[4][2]={{1,0},{-1,0},{0,1},{0,-1}};
map<int,int> ha[N];
int n,m,r,k;
double ans=0.000;
struct xd{
   int x,y;
   double z;
   bool operator < (const xd &a) const {return a.z>z;}
}tmp,nw;
priority_queue<xd> q;
inline int read(){
   int T=0,F=1; char ch=getchar();
   while(ch<'0'||ch>'9'){if(ch=='-') F=-1; ch=getchar();}
   while(ch>='0'&&ch<='9') T=(T<<3)+(T<<1)+(ch-48),ch=getchar();
   return F*T;
}
int main(){
    n=read(),m=read(),r=read(),k=read();
    tmp.x=n/2+1,tmp.y=m/2+1,ha[tmp.x][tmp.y]=1,tmp.z=(min(n-r+1,tmp.x)-max(0,tmp.x-r))*(min(m-r+1,tmp.y)-max(0,tmp.y-r))*1.0000000000000/(n-r+1)*1.0000000000000000/(m-r+1),q.push(tmp);
    for(int i=1;i<=k;++i){
        tmp=q.top(),ans+=tmp.z,q.pop();
        for(int j=0;j<4;++j){
            nw.x=tmp.x+f[j][0],nw.y=tmp.y+f[j][1];
            if(nw.x<1||nw.x>n||nw.y<1||nw.y>m||ha[nw.x].find(nw.y)!=ha[nw.x].end()) continue;
            nw.z=(min(n-r+1,nw.x)-max(0,nw.x-r))*(min(m-r+1,nw.y)-max(0,nw.y-r))*1.0000000000000/(n-r+1)*1.0000000000000000/(m-r+1);
            q.push(nw),ha[nw.x][nw.y]=1;
        }
    }
    printf("%.10lf",ans);
    return 0; 
}

转载于:https://www.cnblogs.com/ljk123-de-bo-ke/p/11318974.html