【HNOI2016模拟4.4】Stage（几何&&倍增）

Description

Solution

这题很明显要根据线性期望率来做。
答案明显是所有的节点的格子被包含在凸包中的期望和，因为贡献是为1的，所以就是所有节点被包含在凸包中的概率和。
我们可以发现，要求被包含在凸包中的概率十分的难求，所以我们可以考虑一下求一个节点不被包含在凸包中的概率，随后在用1-这个概率。
一个团伙不被包含在一个凸包上的概率明显就是这个团伙和出现的观察点构成的凸包，团伙在凸包上。
那么我们只需要枚举凸包上的边出现的概率。
这个就是一个很经典的问题了。首先枚举两条与团伙相连的边（一条容易算重），然后这两条边后面的观察点都不能出现，否则不会出现在凸包上，然后就可以计算概率。
这样$n^3$就有了60分。
然后我们可以极角排序一下，枚举的边都是相邻的，所以一开始枚举第一条边的时候求出对应的概率总和。然后当边移动的时候，影响的节点也只有夹角处的地方，所以可以很快的修改。
但是这样做有除法啊，除法是相当克精度的！！！
那么乘法用一个数据结构维护，比如说很好打的倍增。

Code

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fod(i,a,b) for(i=a;i>=b;i--)
using namespace std;
typedef double db;
const int maxn=1007;
int i,j,k,l,t,n,m,num;
db p,ans,up,down,o,w[maxn],pp[maxn],p0,p1,qian[maxn],hou[maxn],sum[maxn][10],P;
struct node{
    db x,y,ji;
    db p;
}a[maxn],b[maxn],c[maxn];
bool bz[maxn];
bool cmp(node x,node y){return x.ji<y.ji;}
db cha(db x,db y,db xx,db yy){return x*yy-y*xx;}
db dian(db x,db y,db xx,db yy){return x*xx+y*yy;}
db he(int l,int r) {
    if(l>r)return 1;
    db o=1;int i;
    fod(i,9,0)if(l+(1<<i)<=r)o*=sum[l][i],l+=(1<<i);
    return o*sum[r][0];
}
int main(){
//  freopen("fan.in","r",stdin);
    scanf("%d%d",&n,&m);
    fo(i,1,n)scanf("%lf%lf",&a[i].x,&a[i].y);
    p0=1.0;
    fo(i,1,m)scanf("%lf%lf%lf",&b[i].x,&b[i].y,&b[i].p),p0*=1.0-b[i].p;
    fo(i,1,n){
        fo(j,1,m)c[j].x=b[j].x-a[i].x,c[j].y=b[j].y-a[i].y,c[j].p=b[j].p,c[j].ji=atan2(c[j].y,c[j].x);
        sort(c+1,c+1+m,cmp);
        qian[0]=1;fo(j,1,m)qian[j]=qian[j-1]*(1-c[j].p),sum[j][0]=1-c[j].p;
        hou[m+1]=1;fod(j,m,1)hou[j]=hou[j+1]*(1-c[j].p);
        fo(k,1,9)fo(j,1,m-(1<<(k-1)))sum[j][k]=sum[j][k-1]*sum[j+(1<<(k-1))][k-1];
        fo(j,1,m)w[i]+=qian[j-1]*hou[j+1]*c[j].p;w[i]+=qian[m];
        k=1;P=0;
        while(k<m&&cha(c[1].x,c[1].y,c[k+1].x,c[k+1].y)>=0)
        k++,P+=hou[k+1]*c[k].p;
        w[i]+=P*c[1].p;if(k==1)k++;
        fo(j,2,m) {
            if(cha(c[j].x,c[j].y,c[j-1].x,c[j-1].y)<0)P-=hou[j+1]*qian[j-2]*c[j].p;
            P*=1-c[j-1].p;
            while(k%m+1!=j&&cha(c[j].x,c[j].y,c[k%m+1].x,c[k%m+1].y)>=0){
                k=k%m+1;
                if(j<k)P+=qian[j-1]*hou[k+1]*c[k].p;else P+=he(k+1,j-1)*c[k].p;
            }
            w[i]+=P*c[j].p;
            if(j==k)k++;
        }
    }
    fo(i,1,n)ans+=1-w[i];
    printf("%.8lf\n",(db)ans);
}