#include<iostream>
#include<vector>
#include<math.h>
using namespace std;
#define max_v 100000

__int64 p[50000];
vector<int>pr;

__int64 ph(int m)
{
	int x=m;
if(m<50&&p[m]) return p[m];
	__int64 r=1;
	for(int i=0;i<pr.size();i++){
		if(m%pr[i]==0){
			r*=pr[i]-1;m/=pr[i];
			while(m%pr[i]==0){m/=pr[i];r*=pr[i];}
		}
		if(m<50&&p[m]){r*=p[m]; break;}
	}
if(x<50)p[x]=r;
return r;
}

int gcd(int a,int b)
{
 if(b==0)return a;
 return gcd(b,a%b);
}

int main()
{	
	int i;
	for(i=2;i<50000;i++)
	if(!p[i])
	{
		pr.push_back(i);
		int x=i+i;
		while(x<50000){p[x]=1;x+=i;}
	}
	int x;
	memset(p,0,sizeof(p));
	p[1]=1;
	__int64 f;
	 x=1;
	while(cin>>x) {
		f=0;
		for(int i=1;i<=sqrt(x);i++){
			if(x%i==0){
				f+=(x/i)*ph(i);
				if(i*i!=x)f+=i*ph(x/i);
			
			}
	
		}
		printf("%I64d\n",f);
		

	}


}

• Re:你贴一下程序 (381) frkstyc 2005-07-15 01:54:37
• zumba fitness dvd (823) Unpanvary 2012-09-20 02:11:19
  • http://franchbottes.webnode.fr uggs pas cher (1060) Unpanvary 2012-12-15 00:46:48
    • Re:http://franchbottes.webnode.fr uggs pas cher (1072) 18760372288 2019-10-27 22:47:09
  • ugg mini (950) Unpanvary 2012-12-18 08:00:01
    • Shop yon a little jot to fathom cheap jerseys from china what other choices are available. (0) GesSuinsnus 2013-01-30 09:10:11
    • Shop yon a little jot to fathom cheap jerseys from china what other choices are available. (2527) GesSuinsnus 2013-01-30 09:10:16
  • ugg boots ebay (1068) Unpanvary 2013-01-02 02:20:15
  • ugg kensington boots (1066) Unpanvary 2013-01-02 02:20:26
• read more (1061) Unpanvary 2012-09-28 20:52:12
• uggs kopen (493) Unpanvary 2012-11-03 05:21:05
• ugg kensington (1161) Unpanvary 2012-12-19 15:24:51
• ugg handbags (1064) Unpanvary 2013-01-04 05:05:36