UVA11426 GCD - Extreme (II) 欧拉函数`扩展欧几里德应用

这题目太难了 对我这个大菜狗来说，想不到跟欧拉函数有关系，后来看了书，马上懂了，还好懂的快 给了我一点心理安慰，思路很活跃，要求思维也要很活跃，好题目 留下来以后看看

解析：

设f(n)=gcd(1,n)+gcd(2,n)+....+gcd(n-1,n),则所求的答案为 s(n)=f(1)+f(2)+...+f(n),因为这个式子可以变成递推式 s(n)=s(n-1)+f(n),同时又可以找到边界条件 f[2]==s[2],所以只需要求出f(n);

用g(n,i)表示满足 gcd(x,n)=i的 x的个数 （x小于n），则 f(n)=sum{i*g(n,i)};

同时又利用 扩展欧几里德的性质  gcd(x,n)=i  的充要条件是 gcd(x/i,n/i)==1，所以 满足 x/i的解有 phi(n/i)个，说明  g(n,i)=phi(n/i),

接下来就简单咯  直接见代码

#include<iostream>
#include<cstdio>
#include<list>
#include<algorithm>
#include<cstring>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<cmath>
#include<memory.h>
#include<set>

#define ll long long
#define LL __int64
#define eps 1e-8

//const ll INF=9999999999999;

#define inf 0xfffffff

using namespace std;

//vector<pair<int,int> > G;
//typedef pair<int,int> P;
//vector<pair<int,int>> ::iterator iter;
//
//map<ll,int>mp;
//map<ll,int>::iterator p;
//
//vector<int>G[30012];

const int maxp=4000001;

ll phi[maxp],f[maxp],s[maxp];

void init()//注意不能用线性 的模版 因为会爆内存，题目给了10s所以不用担心使用这个模版会超时
{
	memset(phi,0,sizeof(phi));
	phi[1]=1;
	for(int i=2;i<maxp;i++)
		if(!phi[i])
		{
			for(int j=i;j<maxp;j+=i)
			{
				if(!phi[j])
					phi[j]=j;
				phi[j]=phi[j]/i*(i-1);
			}
		}
		return ;
}

int main(void)
{
	init();
	memset(f,0,sizeof(f));
	for(int i=1;i<maxp;i++)
		for(int x=i*2;x<maxp;x+=i)
			f[x]+=i*phi[x/i];
	s[2]=f[2];
	for(int x=3;x<maxp;x++)
		s[x]=s[x-1]+f[x];
	int n;
	while(scanf("%d",&n),n)
	{
		printf("%lld\n",s[n]);
	}
}