期望概率dp lightoj 1038 Race to 1 Again

原创 2015年11月17日 16:54:33

从考虑第一步开始

p(x):代表从x到1的期望

n的所有因子为n,a1,a2....,1(可以打表求出),总个数为size, x = size - 2;

p(n) = (p(a1)+1+p(a2)+1+.....+p(ax)+1)*(1/size) + 1/size + (p(n)+1)*(1/size);

化简式子得到:

p(a1)+p(a2)+,,,+p(ax) = s;

p(n) = (s + size)/(size-1);

利用打表得到的size,少了1所以+1

/********************************************
Author         :Crystal
Created Time   :
File Name      :
********************************************/
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <climits>
#include <string>
#include <vector>
#include <cmath>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <sstream>
#include <cctype>
using namespace std;
typedef long long ll;
typedef pair<int ,int> pii;
#define MEM(a,b) memset(a,b,sizeof a)
#define CLR(a) memset(a,0,sizeof a);
const int inf = 0x3f3f3f3f;
const int MOD = 1e9 + 7;
//#define LOCAL
double dp[200006];
std::vector<int> v[100005];
void init(){
	for(int i=2;i<=100001;i++){
		for(int j=i;j<=100001;j+=i){
			v[j].push_back(i);
		}
	}
	dp[1] = 0;
	for(int i=2;i<=100001;i++){
		double s = 0;
		for(int j=0;j<v[i].size();j++){
			s += dp[i/v[i][j]];
		}
		dp[i] = (s+1)/v[i].size()+1;
		//递推公式,
	}
}
int main()
{
#ifdef LOCAL
	freopen("in.txt", "r", stdin);
//	freopen("out.txt","w",stdout);
#endif
	init();
	int t;cin >>  t;
	int kase = 1;
	while(t--){
		int n;cin >>n;
		printf("Case %d: %.10lf\n",kase++,dp[n]);

	}
	return 0;
}









相关文章推荐

LightOJ 1038 Race to 1 again 期望+概率dp

/* 题目描述:给定正整数n,每次任意选择n的一个因子并除掉,问平均一共要除多少次能把n变成1 思路:设dp[x]表示把x变为1的期望,f[x]表示x的因子个数,根据全期望公...
  • jijijix
  • jijijix
  • 2016年08月27日 10:15
  • 105

lightoj1038 - Race to 1 Again(期望DP)

题意 给出一个1≤N≤1051 \leq N \leq 10^5,每次选其一个约数相除,知道得到结果为1为止,求期望次数; 思路 期望dp,求x平均除多少次得到1;假设x有c个因子(含1和本身...

LightOJ 1038 Race to 1 Again【概率DP】

Race to 1 Again Crawling in process... Crawling failed Time Limit:2000MS     Memory Limit:32768K...
  • Bcwan_
  • Bcwan_
  • 2016年10月03日 16:53
  • 147

LIghtOJ1038---Race to 1 Again(概率dp)

Rimi learned a new thing about integers, which is - any positive integer greater than 1 can be divid...

Light oj 1038 Race to 1 Again(概率dp)

1038 - Race to 1 Again PDF (English) Statistics Forum Time Limit: 2 se...

light oj 1038 - Race to 1 Again(概率DP)

题目链接:http://lightoj.com/volume_showproblem.php?problem=1038 1038 - Race to 1 Again  ...

LIghtOJ 1038 Race to 1 Again

如果只是找n的所有因子,而不需要判断因子是不是素数,应该会快一些。

UVa 11762 - Race to 1 (概率 期望 DP 马尔可夫过程)

UVA - 11762 Race to 1 Time Limit: 10000MS Memory Limit: Unknown 64bit IO Format: %lld & ...

UVA 11762 Race to 1 dp+概率

#include #include #include #include using namespace std; const int maxn=1e6+10; const int maxc=1...

UVA 11762 Race to 1 概率DP

题目大意: 现在对于一个数D, 每次随机从不大于D的素数中选择一个, 如果选择的素数能够整除D, 将D除以这个数得到一个新的D, 否则什么都不做, 问最后将D变为1的选择步数的期望值 大致...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:期望概率dp lightoj 1038 Race to 1 Again
举报原因:
原因补充:

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