aoj 600 数字游戏

转载 2016年08月28日 15:42:35

题目链接:http://icpc.ahu.edu.cn/OJ/Problem.aspx?id=600


题目描述:

Description
cxlove对数字情有独钟,最近又开始玩一个有趣的数字游戏。
首先我们定义一种Lucky number:最高位为1的数字(10进制)。
接下来,会给你n个区间,[Li,Ri]
随机从每一个区间内取出1个整数。
问取出的这n个数中至少有K%是Lucky number的概率是多少。

Input
一个整数 T,表示T组数据。(1<=T<=50)
一个整数n,表示区间的个数 (1<=n<=1000)
接下来n行,每一行两个整数Li,Ri表示区间[Li,Ri],并且保证(1<=Li<=Ri<=10^18)
最后一行是一个整数k (0<=k<=100)

Output
一个实数,表示至少有k%是Lucky Number的概率,小数点后保留6位。

Sample Input
Original Transformed
2
1
1 2
50
2
1 2
9 11
50

Sample Output
Original Transformed
0.500000
0.833333

Source
安徽大学第五届ACM/ICPC程序设计竞赛 现场赛


解题思路:转自http://blog.csdn.net/j_sure/article/details/41624065

警示1:整数的n次方,再也不要用自带的pow()函数。自己写!心情好还能写个快速幂。因为自带的pow()函数是double类型的,精度损失非常大

警示2:区间右端点减去区间左端点,左端点要先减去1。防止左端点被减掉。

警示3:一定要根据dp数组的状态描述仔细地进行初始化。

设dp[i][j]表示前i个区间选到了j个幸运数字的概率。

那么dp[0][0] = 1(0个区间一定是0个幸运数字)

设每个区间选择到幸运数字的概率是v[i]

状态转移方程为:

dp[i][j] = dp[i-1][j-1] * v[] + dp[i-1][j] * (1-v[])

要么是第i个区间选择到了幸运数字,要么是前i-1个区间已经选够了j个幸运数字。

由此我们需要初始化dp[i][0]。因为1个幸运数字都没有选到,所以一路乘即可。

得到至少k%个幸运数字其实就是n * k%想上取整而已。设至少p个。之后把dp[][p] + dp[][p+1] + ... + dp[][n]即可。

现在剩下的问题是预处理区间内幸运数字的个数。这样就可以尽快查询到。(10^18很大啊)。

0~9:1个

10~99:10个

100~999:100个

......

取出左右端点将它们“剪裁”成上述片段即可。这里有用到前缀和。


AC代码:

#include<iostream>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<cstdio>
#include<algorithm>
using namespace std;
long long int T,n,k,num;
double dp[1005][1005];
double rate[1005];
long long int l,r;
long long int ll[20];//ll[i]代表最长i位数字所拥有的以1开头的数字的个数

long long int power(long long int x,long long int y)//x^y
{
    long long int j=1;
    while(y--)
        {
            j*=x;
        }
    return j;
}

void Init()
{
    ll[0]=0;
    for(long long int i=1;i<=18;i++)
    {
        ll[i]=power(10,i-1);
        ll[i]+=ll[i-1];
    }
}

long long int getnum(long long int x)
{
    char s[25];
    sprintf(s,"%lld",x);
    int len=strlen(s);
    if(len==1)
        {
            return x?1:0;
        }
    if(s[0]=='1')
        return x-power(10,len-1)+1+ll[len-1];
    else
        return ll[len];
}

//用这个来求解1开头的数字的个数会超时
/*long long int zg(long long int k)
{
    long long int p=k;
    long long int w=0,q=1;
    while(p)
    {
        w++;
        p/=10;
    }
    for(long long int i=1;i<=w-1;i++)
        q*=10;
    if(k/q==1)
        return 1;
    else
        return 0;
}*/

int main()
{
    Init();
    scanf("%lld",&T);
    while(T--)
    {
        //memset(rate,0,sizeof(rate));
        scanf("%lld",&n);
        for(long long int i=1;i<=n;i++)
        {
            scanf("%lld%lld",&l,&r);
            num=getnum(r)-getnum(l-1);
            rate[i]=1.0*num/(r-l+1);
        }
        scanf("%lld",&k);
        memset(dp,0,sizeof(dp));
        dp[0][0]=1;
        for(long long int i=1;i<=n;i++)
            dp[i][0]=dp[i-1][0]*(1.0-rate[i]);
        for(long long int i=1;i<=n;i++)
        {
            for(long long int j=1;j<=i;j++)
            {
                dp[i][j]+=dp[i-1][j-1]*rate[i]+dp[i-1][j]*(1.0-rate[i]);
            }
        }
        double ans=0.0;
        k=ceil(1.0*k/100.0*n);//ceil函数是求大于等于本身的最小整数
        for(long long int i=k;i<=n;i++)
            ans+=dp[n][i];
        printf("%.6f\n",ans);
    }
    return 0;
}


相关文章推荐

【概率DP】AHU-600 数字游戏

给出n个区间。设'1'开头的数字为幸运数字。从n个区间中取出n个数字。问取出的数字至少有k%是幸运数字的概率。

AOJ-AHU-OJ-191 删除数字

删除数字 Time Limit: 2000 ms   Case Time Limit: 2000 ms   Memory Limit: 64 MB Total Submission: 675   ...

AOJ-AHU-OJ-64 数字三角形

数字三角形问题 Time Limit: 1000 ms   Case Time Limit: 1000 ms   Memory Limit: 64 MB Description ...

AOJ 593分组玩游戏:利用STL 快速排序

n个小朋友在老师的带领下玩游戏。 其中需要分组,老师打算根据大家的衣服颜色分组。 已知衣服颜色用0-9的数字表示。 老师决定每k个人分为一组。每一组的小朋友必须身穿同一种颜色的衣服。 请你告诉...

C600音视频数字监控系统

  • 2012-06-05 13:33
  • 6.68MB
  • 下载

数字游戏

网上看到这样一道面试题,于是没事干

MFC 小游戏之猜数字

这个游戏简单而有趣,在做这个游戏中,又对MFC的一些知识更为了解了,又学到了不少东西。 感兴趣的筒靴可以下载源代码看看! 游戏规则请看下面的说明界面截图! 因水平有限,不好之处请见谅! ...

【cocos2dx-3.0beta-制作flappybird】对不起,我的数学老师是体育老师兼职的—游戏中所用到的数字类封装

一、引言 在游戏进行的过程当中,我们需要实时的显示得分,此外在游戏结束的时候,也需要显示最高分和本局得分。因此本节主要是讲解如何显示这些分数。 二、数字类的设计 这里我们写了一个Numb...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)