# POJ 2480 Longge's problem 欧拉函数

f(M) = ∑d1 * φ(M / d1), d1 | M
f(N) = ∑d2 * φ(N / d2), d2 | N

f(MN) = ∑d * φ(MN / d), d | MN

f(MN)中的项与f(M) * f(N)中的项一一对应

#include<cstdio>
#include<cstring>
using namespace std;

#define MAXN 200000
#define lint __int64
struct Factor { lint b, e; };
Factor f[MAXN]; lint fnum;
lint a[MAXN], p[MAXN], pn;
lint n, ret;

void Prime()
{
lint i, j; pn = 0;
memset(a,0,sizeof(a));
for ( i = 2; i < MAXN; i++ )
{
if ( a[i] == 0 ) p[pn++] = i;
for ( j = 0; j < pn && i*p[j] < MAXN && (p[j]<=a[i] || !a[i]); j++ )
a[i*p[j]] = p[j];
}
}

lint Euler ( lint n )
{
lint ret = n;
for ( int i = 0; p[i] * p[i] <= n; i++ )
{
if ( n % p[i] == 0 )
{
ret = ret - ret / p[i];
while ( n % p[i] == 0 ) n /= p[i];
}
}
if ( n > 1 )
ret = ret - ret / n;
return ret;
}

void split ( lint n )
{
fnum = 0;
for ( int i = 0; p[i] * p[i] <= n; i++ )
{
if ( n % p[i] ) continue;
f[fnum].b = p[i]; f[fnum].e = 0;
while ( n % p[i] == 0 )
{
f[fnum].e++;
n /= p[i];
}
fnum++;
}
if ( n > 1 )
f[fnum].b = n, f[fnum++].e = 1;

}

void DFS ( lint val, int index ) //求n的每一个约数，然后利用欧拉函数
{
if ( index == fnum )
{
ret += Euler(n/val) * val;   //Euler(n/val)的值表示1-n中gcd(n,i）= val的个数
return;
}
for ( lint i = 0, tmp = 1; i <= f[index].e; i++, tmp *= f[index].b )
DFS ( val*tmp, index+1 );
}

int main()
{
Prime();
while ( scanf("%I64d",&n) != EOF )
{
split ( n );
ret = 0;
DFS ( 1, 0 );
printf("%I64d\n",ret);
}
}


#include<cstdio>
#include<cstring>
using namespace std;

#define MAXN 200000
#define lint __int64
struct Factor { lint b, e, mult; };
Factor f[MAXN]; lint fnum;
lint a[MAXN], p[MAXN], pn;

void Prime()
{
lint i, j; pn = 0;
memset(a,0,sizeof(a));
for ( i = 2; i < MAXN; i++ )
{
if ( a[i] == 0 ) p[pn++] = i;
for ( j = 0; j < pn && i*p[j] < MAXN && (p[j]<=a[i] || !a[i]); j++ )
a[i*p[j]] = p[j];
}
}

lint Euler ( lint n )
{
lint ret = n;
for ( int i = 0; p[i] * p[i] <= n; i++ )
{
if ( n % p[i] == 0 )
{
ret = ret - ret / p[i];
while ( n % p[i] == 0 ) n /= p[i];
}
}
if ( n > 1 )
ret = ret - ret / n;
return ret;
}

void split ( lint n )
{
fnum = 0;
for ( int i = 0; p[i] * p[i] <= n; i++ )
{
if ( n % p[i] ) continue;
f[fnum].b = p[i]; f[fnum].e = 0;
f[fnum].mult = 1;
while ( n % p[i] == 0 )
{
f[fnum].e++;
f[fnum].mult *= p[i];
n /= p[i];
}
fnum++;
}
if ( n > 1 )
f[fnum].b = f[fnum].mult = n, f[fnum++].e = 1;

}

int main()
{
Prime(); lint n;
while ( scanf("%I64d",&n) != EOF )
{
split ( n );
lint ret = 1, tmp, sum;
for ( int i = 0; i < fnum; i++ )
{
tmp = 1, sum = Euler(f[i].mult); //所有与f[i].mult互素的数先加起来
for ( int j = 1; j <= f[i].e; j++ )
{
tmp *= f[i].b;
sum += Euler(f[i].mult/tmp) * tmp;
}
ret *= sum;
}
printf("%I64d\n",ret);
}
}


• 本文已收录于以下专栏：

## poj2480（欧拉函数，必须回顾的题）

• zhhx2001
• 2016年06月26日 22:42
• 392

## poj 2480 Longge's problem 积性函数性质+欧拉函数

• sepNINE
• 2015年06月03日 00:11
• 780

## poj_2480

#include #include #include #include #include #include #include using namespace std; #define M...
• bryant03
• 2015年08月01日 00:03
• 486

## POJ2480 欧拉函数的应用

POJ 2480 欧拉函数的应用题目链接 乍一看和欧拉函数没什么关系，但数论就是这样，处处联系。我们可以枚举i(1...
• MasterAn
• 2017年11月11日 20:09
• 70

## POJ 2480 Longge's problem (欧拉函数)

• qian99
• 2014年03月15日 16:18
• 878

## 【poj】 2480 Longge's problem （欧拉函数）

【poj】 2480 Longge’s problem （欧拉函数） 先打表找规律，发现对于一个n，最大公因子相同的项很多，于是考虑从因子角度来求，最大公因数为1的项的数量显然为phi(n),对于某个...
• xuziling_
• 2017年08月21日 11:19
• 52

## POJ 2480 Longge's problem 欧拉函数

• Tsaid
• 2012年03月25日 00:27
• 1146

## POJ2480 积性函数

• qq_36553623
• 2017年08月05日 14:51
• 119

## POJ 2480 Longge's problem （欧拉函数）

Longge’s problem Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 8107 Accepted...
• helloiamclh
• 2016年07月22日 21:48
• 291

## POJ - 2480 Longge's problem 欧拉函数

• ACVector
• 2017年09月05日 15:26
• 81

举报原因： 您举报文章：POJ 2480 Longge's problem 欧拉函数 色情 政治 抄袭 广告 招聘 骂人 其他 (最多只允许输入30个字)