GWind概率求解 皇星客栈

G-Wind

Problem Description

N棵树和M个蘑菇排成一排. 每棵树有自己的坐标和高度, 每个蘑菇有自己的坐标和权值. 一阵大风刮来, 每棵树有Li%的概率往左边倒, 有Ri%的概率往右边倒, 有(1 – Li% - Ri%)的概率不倒. 如果一棵高度为H坐标为X的树往左边倒下, 那么坐标在区间[X – H, X)内的蘑菇就会被砸到, 如果往右边倒下, 那么坐标在区间(X, X + H]内的蘑菇就会被砸到. 求不被树砸到的所有的蘑菇的权值之和的期望值。

Input

第一行包含2个整数N, M.
接下来N行, 每行4个整数X, H, L, R. 表示这棵树的坐标, 高度和往左右倒的概率. 即这棵树有L%的概率往左倒, R%的概率往右倒.
接下来M行, 每行2个整数Y, W. 表示这个蘑菇的坐标和权值.
[数据范围]
对于20% 的数据, 1 <= N, M <= 1000;
对 于100% 的数据, 1 <= N <= 10^5, 1 <= M <= 10^4, 树和蘑菇的坐标的绝对值不超过10^9, 树的高度在[1, 10^9]范围内, 每个蘑菇的权值在[1, 1000]范围内, 每棵树的L,R 满足 0 <= L, R, L + R <= 100.

Output

输出一个整数, 为最后答案. 保留4位小数.

Sample Input

1 1
2 2 50 50
1 1

Sample Output

0.5000
第一种方法：比较暴力

#include<stdio.h>
#include<stdlib.h>

typedef struct tree{
    int X;
    int H;
    double L;
    double R;
}node;
node T[100001];

typedef struct mogu{
    int Y;
    double W;
}mu;
mu Mo[10001];

int cmp1( const void *a,const void *b )
{
    node *aa = (node *)a;
    node *bb = (node *)b;
    return aa->X - bb->X;
}

int cmp2( const void *a,const void *b )
{
    mu *aa = (mu *)a;
    mu *bb = (mu *)b;
    return aa->Y - bb->Y;
}

int main( )
{
    int N;
    int M;
    int i;
    int j;
    double sum;
    
    while( scanf("%d %d",&N,&M) == 2 )
    {
        
        sum = 0;
        for( i = 0 ; i < N ; i++ )
            scanf("%d %d %lf %lf",&T[i].X,&T[i].H,&T[i].L,&T[i].R);
        for( i = 0 ; i < M ; i++ )
            scanf("%d %lf",&Mo[i].Y,&Mo[i].W);
        
        qsort(T,N,sizeof(T[0]),cmp1);
        qsort(Mo,M,sizeof(Mo[0]),cmp2);
        
        for( i = 0 ; i < N ; i++ )
            for( j = 0 ; j < M ; j++ )
            {
                if( Mo[j].Y < T[i].X-T[i].H )
                    continue;
                else if( Mo[j].Y > T[i].X+T[i].H )
                    continue;
                else{
                    if( Mo[j].Y < T[i].X )
                        Mo[j].W *= (1-T[i].L/100);
                    else
                        Mo[j].W *= (1-T[i].R/100);
                }
            }
            for( j = 0 ; j < M ; j++ )
                sum += Mo[j].W;
            printf("%.4lf\n",sum);
    }
    return 0;
}