ACM-ICPC 2018 南京赛区网络预赛--J Sum(数论)

题意:

6 = 1*6 = 2*3 = 3*2 = 6*1,12 = 2*6 = 6*2(12可以被2^2整除),这样 f(6) = 4,f(12) = 2,求 f(1) ... f(n) 的和。

题解:

其实我们把数字拆成质数相乘就会发现只要这个质数的个数不超过2个就有值,比如 6 的有 2,3,即 6 = 1 * 2^1 * 3^1, 24 = 1 * 2^3 * 3^1。

而且我们可以发现数字如果是质数则 f(n) = 2。

我们假设 n = p^i * x,如果 i == 2 则 f(n) = f(x);如果 i == 1 则 f(n) = f(p)*f(x)。

通过这个公式打个表就可以输出答案了。

#include 
#include  
#include   
#include   
#include    
#include    
#include    
#include    
#include     
#include     
#include     
#include     
#include      
#include       
#include       
#define mem(a, b) memset(a, b, sizeof(a))
#define inf 0x7ffffff
#define pi acos(-1)
using namespace std;
typedef long long ll;

const int maxn = 2e7+10;
int vis[maxn], p[maxn], f[maxn], sum[maxn];

void init(){
	f[1] = 1;
	mem(vis, 0);
	int k = 0;
	for(int i = 2; i < maxn; i++){
		if(vis[i] == 0){
			p[k++] = i;
			f[i] = 2;
		}
		for(int j = 0; p[j]*i < maxn && j < k; j++){
			int now = p[j] * i;
			vis[now] = 1;
			if(i % p[j]){
				f[now] = f[i]*f[p[j]];
			}
			else if(i % (p[j]*p[j]) == 0){
				f[now] = 0;
			}
			else{
				f[now] = f[now/(p[j]*p[j])];
				break;
			}
		}
	}
	for(int i = 1; i < maxn; i++){
		sum[i] = sum[i-1]+f[i];
	}
}

int main(){
	int t;
	init();
	scanf("%d", &t);
	while(t--){
		int n;
		scanf("%d", &n);
		printf("%d\n", sum[n]);
	}
}

 

你可能感兴趣的:(ICPC题解)