P3708 koishi的数学题

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

const int N=1000010;

bool flag[N];
long long f[N],ans[N];
int tot,n,prim[N],t[N];

void init() {
	f[1]=1;
	for(int i=2; i<=n; i++) {
		if (!flag[i]) {
			prim[++tot]=i;
			f[i]=i+1;
			t[i]=i;
		}
		for(int j=1; j<=tot; j++) {
			int k=prim[j]*i;
			if(k>n)
				break;
			flag[k]=1;
			if(i%prim[j]==0) {
				t[k]=t[i]*prim[j];
				if(k==t[k])
					for(int l=1; l<=k; l*=prim[j])
						f[k]+=l;
				else
					f[k]=f[t[k]]*f[k/t[k]];
				break;
			}
			t[k]=prim[j];
			f[k]=f[prim[j]]*f[i];
		}
	}
}

int main() {
	scanf("%d",&n);
	init();
	ans[1]=n-1;
	for(int i=2; i<=n; i++)
		ans[i]=ans[i-1]+n-f[i];
	for(int i=1; i<=n; i++)
		printf("%lld ",ans[i]);
	printf("\n");
	return 0;
}