uva11426(数学题)

  题目的大概意思是给你一个n，要求gcd(1,2) + gcd(1,3) + ... +  gcd(i, j) + ... + gcd(n-1, n)， 1<= i < j <= n。

  先按照递推的思路，设f(n) = gcd(1,n) + gcd(2,n) + ... + gcd(n-1,n)， 那么ans(n) = ans(n-1) + f(n), 所以题目的关键就是求f(n)，求f(n)就是要找到n的约数，一个一个的找，那样的话是n^2的复杂度，n的范围是4000000，所以很明显是不行的。

  我们可以设g(n,i) 为gcd(n,x) = i的个数，那么g(n,i) = phi(n/i)，所以f(n) = sum(g(n/i) * i), ( n % i == 0)，这样求的话也是需要得到n所有的约数，然后想到约数可以在筛选法的时候可以得到，所有f就可以在筛选的时候做出来。

#include "stdio.h"
#include "stdlib.h"
#include "string.h"
#include "math.h"
#include "time.h"
#include <string>
#include <queue>
#include <stack>
#include <vector>
#include <map>
#include <algorithm>
#include <iostream>
using namespace std;

#define MAXM 1
#define MAXN 1
#define max(a,b) a > b ? a : b
#define min(a,b) a < b ? a : b
#define Mem(a,b) memset(a,b,sizeof(a))
int Mod = 1000000007;
double pi = acos(-1.0);
double eps = 1e-6;

typedef struct{
	int f,t,w,next;
}Edge;
Edge edge[MAXM];
int head[MAXN];
int kNum;
typedef long long LL;

void addEdge(int f, int t, int w)
{
	edge[kNum].f = f;
	edge[kNum].t = t;
	edge[kNum].w = w;
	edge[kNum].next = head[f];
	head[f] = kNum ++;
}

int phi[4000005];
LL f[4000005];
LL s[4000005];

void getPhi(){
	Mem(phi, 0);
	for(int i = 2; i <= 4000000; i ++){
		if( phi[i] ) continue;
		for(int j = i; j <= 4000000; j += i){
			if( !phi[j] ) phi[j] = j;
			phi[j] = phi[j] / i * ( i - 1 );
		}
	}
}

void getf(){
	Mem(f, 0);

	phi[1] = 0;
	for(int i = 1; i <= 4000000; i ++){
		for(int j = i; j <= 4000000; j += i){
			f[j] += i * phi[j/i];
		}
	}
}

void getS(){
	getPhi();
	getf();
	for(int i = 2; i <= 4000000; i ++){
		s[i] = s[i-1] + f[i];
	}
}

int N;




int main()
{
//	freopen("d:\\test.txt", "r", stdin);
	getS();
	while(cin>>N){
		if( N <= 0 )
			break;
		cout<<s[N]<<endl;
	}

	return 0;
}