uva 11762

最新推荐文章于 2017-01-06 16:17:11 发布

阿龙的圈子

最新推荐文章于 2017-01-06 16:17:11 发布

阅读量848

点赞数

CC 4.0 BY-SA版权

分类专栏： DP 文章标签： ACM

本文链接：https://blog.youkuaiyun.com/DoJintian/article/details/44519301

DP 专栏收录该内容

46 篇文章

订阅专栏

本文介绍了一种结合素数筛法与动态规划算法的解决方案，用于处理与整数分解相关的数学问题。通过预先生成素数标记，并利用动态规划优化计算过程，实现了高效的求解方法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

<pre name="code" class="cpp">#include <bits/stdc++.h>

using namespace std;
#define maxn 1000000+10

bool isprime[maxn];
int  mark[maxn];
double dp[maxn];
int T, n;

void make_prime(int num)
{
    int tmp=sqrt(num+0.5);
    memset(isprime,true,sizeof(isprime));
    for(int i=2; i<=tmp; i++)
        if(isprime[i])
            for(int j=i*i; j<=num; j+=i)
                isprime[j]=false;
    isprime[0] = false;
    isprime[1] = false;

    int cnt = 0;
    for(int i=2; i<=num; i++)
        if(isprime[i]) mark[i] = ++cnt;
        else mark[i] = cnt;
}

double solve(int m)
{
    int k = 0;
    if(m == 1) return 0;

    if(dp[m]) return dp[m];
    for(int i=2; i<=m; i++)
        if(isprime[i] && m%i == 0)
            {
                dp[m] += solve(m/i)/mark[m];
                k++;
            }

    dp[m] += mark[m]/(double)k;
    return dp[m];
}


int main()
{
    freopen("input.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
    cin>>T;
    make_prime(maxn);
    while(T--)
    {
        memset(dp, 0, sizeof(dp));
        cin>>n;
        double ans = solve(n);
        printf("%.10lf\n", ans);
    }
    return 0;
}