POJ 2480 Longge's problem 欧拉函数

原创 2012年03月25日 00:27:42
题意: Given an integer N(1 < N < 2^31),you are to calculate ∑gcd(i, N) 1<=i <=N.

题解:

公式:f(N)=∑x*φ(N/x),x | N (x是N的约数)
因为在1···N中,gcd(i,N) = x, 的个数的等于φ(N / x)

另外还可以利用函数的积性:

对于正整数n的一个函数 f(n),当中f(1)=1且当a,b互质,f(ab)=f(a)f(b),在数论上就称它为积性函数。若某函数f(n)符合f(1)=1,且就算a,b不互质,f(ab)=f(a)f(b),则称它为完全积性函数。


不妨令M, N互素
f(M) = ∑d1 * φ(M / d1), d1 | M
f(N) = ∑d2 * φ(N / d2), d2 | N

f(MN) = ∑d * φ(MN / d), d | MN
因为M, N互素,则每个d都可以唯一分解为M中的因子d1, 和N中的因子d2
即d = d1 * d2, d1 | M, d2 | N, d1与d2互素
则d * φ(MN / d) = d1 * d2 * φ(M / d1) * φ(N / d2)
f(MN)中的项与f(M) * f(N)中的项一一对应

解法一:47MS
#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);
    }
}


解法二:利用积性16ms
#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(欧拉函数,必须回顾的题)

欧拉函数的思想!!!! 刚开始写跪了,一直wa,太晚了,就只好先这样了 #include #include #include #include #define ll long long usin...
  • zhhx2001
  • zhhx2001
  • 2016年06月26日 22:42
  • 392

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

题意: 求f(n)=∑gcd(i, N) 1 分析: f(n)是积性的数论上有证明,且f(n)=sigma{1 代码: //poj 2480 //sep9 /* f(pi^ai) = Φ(pi^ai...
  • sepNINE
  • sepNINE
  • 2015年06月03日 00:11
  • 780

poj_2480

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

POJ2480 欧拉函数的应用

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

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

题意:求Σgcd(i,N)。 思路:若gcd(A,N)=i,那么gcd(A/i,N/i)=1,因此,对于N的每一个因子,sum+=φ(N/i)*i。 代码: #include #include ...
  • qian99
  • qian99
  • 2014年03月15日 16:18
  • 878

【poj】 2480 Longge's problem (欧拉函数)

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

POJ 2480 Longge's problem 欧拉函数

题意: Given an integer N(1 题解: 公式:f(N)=∑x*φ(N/x),x | N (x是N的约数) 因为在1···N中,gcd(i,N) = x, 的个数的等于φ(...
  • Tsaid
  • Tsaid
  • 2012年03月25日 00:27
  • 1146

POJ2480 积性函数

传送门 先前的大神已经写得很好了。转载下。 在数论中的积性函数:对于正整数n的一个函数 f(n),当中f(1)=1且当a,b互质,f(ab)=f(a)f(b),在数论上就称它为积性...
  • qq_36553623
  • 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
  • helloiamclh
  • 2016年07月22日 21:48
  • 291

POJ - 2480 Longge's problem 欧拉函数

题目:给出N,求 ∑gcd(i, N) 1 0 思路: phi(p^k)=(p-1)*p^(k-1) 考虑如果N=P^k的时候,那么F[N]=k*p^(k-1)*(p-1)+p^k...
  • ACVector
  • ACVector
  • 2017年09月05日 15:26
  • 81
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:POJ 2480 Longge's problem 欧拉函数
举报原因:
原因补充:

(最多只允许输入30个字)