ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Race to 1 Again(期望dp&除因数到1)

2021-06-08 17:05:09  阅读:222  来源: 互联网

标签:约数 Again cnt const int Race include dp define


题目:https://vjudge.z180.cn/problem/LightOJ-1038#author=0

题意:一个数n除以它的因数k,n=n/k,直到n=1,问除以回合的期望是多少。

题解:因为题目给的T很大,n也很大,所以要在T外面把结果都预处理好,从小往大推就行,f[4]=(f[1]+f[2]+f[4)/cnt+1,cnt是4的约数有几个,这里就是3,左右两边同时乘以cnt,然后对f[4]移项得,f[4]=(f[1]+f[2]+cnt)/(cnt-1),按照这个公式直接暴力递推下去就行。

#include <algorithm>
#include <bitset>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <deque>
#include <functional>
#include <iostream>
#include <map>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <vector>
//#include <unordered_map>
//#include <unordered_set>
//#include <bits/stdc++.h>
//#define int long long
#define pb push_back
#define PII pair<int, int>
#define mpr make_pair
#define ms(a, b) memset((a), (b), sizeof(a))
#define x first
#define y second
typedef long long ll;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
const int N = 1e6 + 7;
using namespace std;


double f[N];
// signed main(){
int main(int argc, char const *argv[]) {
    for (int i = 2; i <= 100000; i++) {
        int cnt = 0;//记录约数个数
        double s = 0;//记录约数期望的和
        for (int j = 1; j <= sqrt(i); j++) {
            if (i % j == 0) {
                cnt++;
                s += f[j]; 
                if (i / j != j) {//判断这俩约数是否相同
                    cnt++;
                    s += f[i / j];
                }
            }
        }
        f[i] = (cnt + s) / (cnt - 1); //公式
    }
    int T;
    int Case=0;
    scanf("%d",&T);
    while(T--){
        int n;
        scanf("%d",&n);
        printf("Case %d: ",++Case);
        printf("%lf\n",f[n]);
    }
    return 0;
}

标签:约数,Again,cnt,const,int,Race,include,dp,define
来源: https://blog.csdn.net/Are_you_ready/article/details/117707646

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有