LightOJ - 1038 Race to 1 Again 基础期望概率 dp

最新推荐文章于 2019-11-06 11:12:43 发布

swust_fang

最新推荐文章于 2019-11-06 11:12:43 发布

阅读量154

点赞数 1

分类专栏： ACM

本文链接：https://blog.csdn.net/swust5120171204/article/details/96449629

版权

ACM 专栏收录该内容

73 篇文章 0 订阅

订阅专栏

传送门

刚刚学习期望&概率

我们设数X的期望改变次数为P[X]

如果要求X的期望，很容易想到找x的因子;

可以得到下式

$P[x]=\frac{\sum (P[ai]+1)+P[x]}{cnt}$ ,cnt为X因子个数，ai为X的因子

可以这么理解，当因子ai为1时，因为除于1不改变其期望值,上式上边分子部分为1+P[1]+P[x]

不为1时，上部分为p[ai]+1，相当于把状态转换过去，并且期望次数加1

分母为cnt，即为包括1和n的因子个数，因为该点可以由所有因子转化而来
经过化简可得

$P[x]=\frac{cnt+\sum P[i]}{cnt-1}$

哈哈，这个时候就是可以用素筛的思想来进行求解。

预处理所有p[i]

第一次自己推出来，开心

#include<stdio.h>
#include<string.h>
#include<string>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<set>
#include<stack>
#include<vector>
#include<map>
#include<queue>
#include<functional>
#pragma comment (linker,"/STACK:102400000,102400000")
#define myself i,l,r
#define lson i<<1
#define rson i<<1|1
#define Lson i<<1,l,mid
#define Rson i<<1|1,mid+1,r
#define half (l+r)/2
#define lowbit(x) x&(-x)
#define min4(a,b,c,d) min(min(a,b),min(c,d))
#define min3(x,y,z) min(min(x,y),min(y,z))
#define max4(a,b,c,d) max(max(a,b),max(c,d))
#define max3(x,y,z) max(max(x,y),max(y,z))
#define pii make_pair
#define pr pair<int,int>
typedef long long ll;
typedef unsigned long long ull;
const int inff=0x3f3f3f3f;
const long long inFF=9223372036854775807;
const int dir[4][2]={0,1,0,-1,1,0,-1,0};
const double eqs=1e-9;
const double E=2.718281828459;
const double pi=acos(-1.0);
const int maxn=1e5+5;
const int mod=1e9+7;
using namespace std;
double p[maxn];
int cnt[maxn];
void init()
{
    p[1]=0;
    for(int i=2;i<maxn;i++)
    {
        p[i]=(p[i]+cnt[i]+2)/(cnt[i]+1);
        for(int j=2*i;j<maxn;j+=i)
            p[j]+=p[i],cnt[j]++;
    }
}
int main()
{
    int t,n,cas=0;
    init();
    cin>>t;
    while(t--)
    {
        cin>>n;
        printf("Case %d: %.7f\n",++cas,p[n]);
    }
    return 0;
}

AC代码：