【ACM-ICPC 2018 南京赛区网络预赛 J】Sum

【链接】 我是链接,点我呀:) 【题意】

在这里输入题意

【题解】

线性筛求出每个数的最小质因子x for 从1-n 对于i,它的最小质因子为x 考虑i=a*b 如果i能被x^3整除那么这x怎么分配给它的两个因子a,b都不行,(都有一个因子会分配至少2个以上,那么我们提出来一个x^2的话,a或b就能被x^2整除了)。所以ans[i] = 0 如果i能被x^2整除(先判断x^3再判断x^2这样能保证x就只能被x^2整除),那么ans[i] = ans[i/(x^2)] 因为这个x显然只能分到两边去(a,b各分配一个,不能都分配到a或都分配到b) 所以答案和i/(x^2)的一样.i/(x^2)有多少个答按a'*b',那么把a'和b'都各乘上x就ok 如果i只能被x整除 ans[i] = ans[i/x]*2 因为这个x有两种可能,分到a或者b。所以i/x的每种方案a',b'。加一个x乘在a'或者b'都能凑成i的方案.且乘1个x显然不会让它有平方因子。

ans[1]=1 然后按照上面的规则递推就好 (再求个前缀和

【代码】

/*
	N在const里面,可以改动;
	求出2..N之间的所有质数;
	放在zsb数组里面;
	时间复杂度是O(N)的;
	但是需要空间也是O(N)的;
*/

#include <bits/stdc++.h>
using namespace std;

const int N = 2e7;
bool iszs[N + 10];
int mizs[N+10];
int ans[N+10];

vector <int> zsb;

int main() {
    #ifdef LOCAL_DEFINE
	    freopen("rush_in.txt", "r", stdin);
	#endif
	memset(iszs, true, sizeof iszs);
	for (int i = 2; i <= N; i++) {
		if (iszs[i]) zsb.push_back(i);
		int len = zsb.size();
		for (int j = 0; j <= len - 1; j++) {
			int t = zsb[j];
			if (i*t > N) break;
			iszs[i*t] = false;
			mizs[i*t] = t;
			if (i%t == 0) break;
		}
	}
    ans[1] = 1;
    for (int i = 2;i <= N;i++){
        int t = mizs[i];
        if (t<=1000 && (i%(t*t*t)==0)){
            ans[i] = 0;
        }else if (t<=10000&&(i%(t*t)==0)){
            ans[i] = ans[i/(t*t)];
        }else {
            ans[i] = ans[i/t]*2;
        }
    }

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

来源：oschina

链接：https://my.oschina.net/u/4387995/blog/3837942

标签

black