UVA11762 Race to 1

关注

原创最新推荐文章于 2026-03-05 00:20:27 发布 455 阅读 AI 写同款· GEO 优化 ›

从0开始使用ArcGIS实现克里金（Kriging）插值（适用数学建模竞赛，详解）适用：数学建模学习中遇到需要进行克里金（Kriging）插值的任务。希望尽快上手，选择出恰当的半函数模型，并获得插值后的图片。详细步骤：一、准备数据在excel中准备需要进行克里金插值的内容，应该包括至少三列：X坐标，Y坐标，数据。在每一列的列首做好标记，方便导入ArcMap后对应。注意存储时存为xls文件。二、将数据导如导ArcMap 点击下面这个标记，把数据文件添加到ArcMap。如果是初次使用的化，可能会不知道进入添加数据对话框之后如何操作。这时应该先点击红框标记的图标，连接到文件夹。阅读详情

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int N=1000005;
int p[N],tmp[N],maxn;
double dp[N];
void prime(){
    for(int i=2;i<N;i++)
        for(int j=i+i;j<N;j+=i)
            p[j]=1;
    for(int i=2;i<N;i++){
        if(!p[i]){
            tmp[maxn++]=i;
        }
    }
}
double dfs(int n){
    if(dp[n]!=-1) return dp[n];
    dp[n]=0;
    if(n==1) return dp[n];
    int tot=0,cnt=0;
    for(int i=0;i<maxn&&tmp[i]<=n;i++){
        tot++;
        if(n%tmp[i]==0){
            cnt++;
            dp[n]+=dfs(n/tmp[i]);
        }
    }
    dp[n]=(dp[n]+tot)/cnt;
    return dp[n];
}
//令 dp[ i ] 表示 D = i 的时候的期望, 即从i 转移到1 的次数期望
//dp(n)=(1−cnt/tot)∗dp(n)+∑dp(n/prime)∗(1/tot)
//cnt素数质因子数,tot总素数\
y为x的素数因子，sum[x]表示1到x中有多少个素数，p[x]表示x\
的素数因子个数
//dp[x] = 1+dp[x]*( (sum[x]-p[x])/sum[x] )+∑ (dp[x/y] / sum[x] )\
========>>>>dp[x]=(tot+∑ (dp[x/y])/cnt

int main()
{
    int T,n,cas=0;
    //memset(dp,-1,sizeof(dp));
    for(int i=0;i<N;i++) dp[i]=-1;//不可用上面那句初始化
    //cout<<dp[9]<<endl;
    cin>>T;
    prime();
    //for(int i=0;i<50;i++)cout<<i<<' '<<p[i]<<endl;
    while(T--){
        scanf("%d",&n);
        printf("Case %d: %.10lf\n",++cas,dfs(n));
    }
    return 0;
}

收录于