暑末 Day1 Face

2 篇文章 0 订阅
0 篇文章 0 订阅
//伪概率DP
#include<bits/stdc++.h>
using namespace std;
int hm, high, low, n, m, c, a, b, t;
int delta;
double dp[810][16001], sum[16001];

inline int read(){
    int num = 0;
    char c;
    while((c = getchar()) == ' ' || c == '\r' || c == '\n');
    num = c - '0';
    while(isdigit(c = getchar()))   num = num*10 + c - '0';
    return num;
}

inline void getInit(){
    memset(dp, 0, sizeof(dp));
}

int main(){
    t = read();
    while(t--){
        m = read(); c = read();
        n = read(); a = read();
        b = read();
        getInit();
        delta = b - a + 1;
        hm = m + c;
        low = n * a;
        high = n * b;
        if(high < m){
            printf("0.00000\n");
            continue;
        }
        else if(low > hm){
            printf("0.00000\n");
            continue;
        }
        for(int i = a; i <= b; i++) dp[1][i] = 1;
        for(int i = 1; i <= n; i++){
            dp[i][a*i] = 1;
            dp[i][b*i] = 1;
        }   
        for(int i = 2; i <= n; i++){
            memset(sum, 0, sizeof(sum));
            for(int j = i*a-b; j <= i*b-a; j++) sum[j] = sum[j-1] + dp[i-1][j];
            for(int j = i*a+1; j <= i*b-1; j++)
                //for(int k = a; k <= b; k++)
                    //if(j - k >= 0)
                    dp[i][j] += sum[j-a] - sum[j-b-1]; //dp[i-1][j-k];
        }

        double ans = 0;
        long long final = 1;
        for(int i = 1; i <= n; i++) final *= (long long)delta;
        for(int i = m; i <= hm; i++)    ans += dp[n][i];
        ans /= (double)final;
        printf("%.5lf\n", ans);
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值