LightOJ 1038 Race to 1 Again

每次一个数,问不断除以约数,到1的时候,除的个数的期望。

怕超时,打表,递推公式:

(n的约数的个数(不含本身))*a[n]=a[n的约数(除了本身)](求和)+约数的个数(包含本身);

用a[8]举个例子;

a[8]=1/4*(a[1]+1)+1/4*(a[2]+1)+1/4*(a[4]+1)+1/4*(a[8]+1);

每次约去就让次数加一~

完美~

#include <iostream>
#include <functional>
#include <algorithm>
#include <complex>
#include <cstdlib>
#include <cstring>
#include <fstream>
#include <iomanip>
#include <sstream>
#include <utility>
#include <bitset>
#include <cctype>
#include <cstdio>
#include <limits>
#include <memory>
#include <string>
#include <vector>
#include <cmath>
#include <ctime>
#include <queue>
#include <stack>
#include <list>
#include <map>
#include <set>
using namespace std;

int main()
{
    double a[100005];
    a[1]=0;
    for (int i=2;i<=100000;i++)
    {
        double temp=0;
        double cnt=0;
        for (int j=1;j*j<=i;j++)
        {
            if (i%j==0)
            {
                temp+=a[j];
                cnt++;
                if (i/j!=j&&j!=1)
                {
                    temp+=a[i/j];
                    cnt++;
                }
            }
        }
        a[i]=(temp+cnt+1.0)/cnt;
    }
//    for (int i=2;i<=100000;i++)
//    {
//        if (i%10==0) system("pause");
//        printf ("%d=%f\n",i,a[i]);
//    }
    int T,ncas=1;
    scanf ("%d",&T);
    while (T--)
    {
        int n;
        scanf ("%d",&n);
        printf ("Case %d: %.10f\n",ncas++,a[n]);
    }
    return 0;
}

时间: 2024-10-13 04:43:58

LightOJ 1038 Race to 1 Again的相关文章

LightOJ 1038 Race to 1 Again 期望 记忆化dp

题目链接:点击打开链接 1038 - Race to 1 Again PDF (English) Statistics Forum Time Limit: 2 second(s) Memory Limit: 32 MB Rimi learned a new thing about integers, which is - any positive integer greater than 1 can be divided by its divisors. So, he is now playin

Lightoj 1038 - Race to 1 Again (概率DP)

题目链接: Lightoj  1038 - Race to 1 Again 题目描述: 给出一个数D,每次可以选择数D的一个因子,用数D除上这个因子得到一个新的数D,为数D变为1的操作次数的期望为多少? 解题思路: 概率DP咯,对于只知道期望是:E(X) = X1*p(X1) + X2*p(X2) + …… + Xn*p(Xn)的窝,拿这个题目没有一点办法.然后看了讨论版,发现总会有一些神人存在. 求操作次数的期望时,先设定第i个因子给期望的贡献为Ti,那么有:E = (T1 + T2 + T3

LightOJ 1038 - Race to 1 Again 【DP】

题目链接:http://www.lightoj.com/volume_showproblem.php?problem=1038 题意:题目很短,不叙述了. 解法:dp 代码: #include <stdio.h> #include <ctime> #include <math.h> #include <limits.h> #include <complex> #include <string> #include <functio

LightOJ 1038 - Race to 1 Again(期望+DP)

题目链接:http://lightoj.com/volume_showproblem.php?problem=1038 题意是:给你一个N (1 ≤ N ≤ 105) 每次N都随机选一个因子d,然后让N=N/d, 求N变成1的次数的期望: 当 N = 2 时 2有两个因子:1,2 E[2] = E[1]/2 + E[2]/2 + 1;因此可以求出E[2]; 当N = 8 时 8有4个因子1 2 4 8; E[8] = E[1]/4 + E[2]/4 + E[4]/4 + E[8]/4+ 1;因此

LightOJ 1038 - Race to 1 Again (给一个数,用这个数的因数除以这个数,直到为1时,求除的次数的期望。)(概率)

题意:http://www.lightoj.com/volume_showproblem.php?problem=1038 题意:给一个数,用这个数的因数除以这个数,直到为1时,求除的次数的期望. 设一个数的约数有M个,E[n] = (E[a[1]]+1)/M+(E[a[2]]+1)/M+...+(E[a[M]]+1)/M 一个数最大的约数是它自己. 则有,E[n] = (E[a[1]]+1)/M+(E[a[2]]+1)/M+...+(E[n]+1)/M (M-1)*E[n]=E[a[1]]+E

LightOJ - 1038 Race to 1 Again 递推+期望

题目大意:给出一个数,要求你按一定的规则将这个数变成1 规则如下,假设该数为D,要求你在[1,D]之间选出D的因子,用D除上这个因子,然后继续按该规则运算,直到该数变成1 问变成1的期望步数是多少 解题思路:递推,设该数为D,有N个因子,分别是1,n1,n2,n3-nn-2,D, 那么选到每个因子的概率都是1/N,除非选到D,不然选到其他因子的话都要多1步,然后再计算D除以该因子的期望 这就能得到公式了,设dp[D]为数D按规则变成1的期望步数 那么dp[D] = 1/N * (dp[D/1]

LightOJ 1326 - Race(第二类斯特林数啊 )

题目链接:http://lightoj.com/volume_showproblem.php?problem=1326 百度百科:斯特林数 ACdreamer:第一类Stirling数和第二类Stirling Disky and Sooma, two of the biggest mega minds of Bangladesh went to a far country. They ate, coded and wandered around, even in their holidays.

[LOJ 1038] Race to 1 Again

C - Race to 1 Again Time Limit:2000MS     Memory Limit:32768KB     64bit IO Format:%lld & %llu Description Rimi learned a new thing about integers, which is - any positive integer greater than 1 can be divided by its divisors. So, he is now playing w

lightoj 1326 - Race dp+预处理

1326 - Race PDF (English) Statistics Forum Time Limit: 1 second(s) Memory Limit: 32 MB Disky and Sooma, two of the biggest mega minds of Bangladesh went to a far country. They ate, coded and wandered around, even in their holidays. They passed severa