POJ 2773 Happy 2006 欧拉函数的应用

题意:输入两个数m,k。求与m互素的第k个数。
题解:我们知道,如果a和m互素,那么k*m+a与m互素。也就是说【1,m-1】,与【k*m+1,k*m+m-1】区间存在一个对应关系:【1,m-1】中与m互素的个数 = 【k*m+1,k*m+m-1】中与m互素的个数。

另外,值得注意的是,我们在求欧拉函数的过程中实际上求出了m所有的素因子,既然如此,就可以通过筛选法筛掉所有与m不互素的数。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<ctime>
#include<algorithm>
using namespace std;

bool pri[1000009];

int Euler ( int n )
{
    int i, j, m = n, ret = n;
    memset(pri,0,n+1);

    for ( i = 2; i * i <= n; i++ )
    {
        if ( n % i == 0 )
        {
            n /= i;
            for ( j = 1; i * j <= m; j++ ) //用筛选法标记所有与n不互素的数
                pri[i*j] = 1;

            ret = ret - ret / i;
            while ( n % i == 0 )
                n = n / i;
        }
    }
    if ( n > 1 )  // 当n>1的时候,n本身也是一个素因子
    {
        for ( j = 1; n * j <= m; j++ )
            pri[n*j] = 1;
        ret = ret - ret / n;
    }
    return ret;
}


int main()
{
    int m, k;
    while ( scanf("%d%d",&m, &k) != EOF )
    {
        int p = Euler(m);
        int tmp = k / p;
        if ( k % p == 0 ) tmp--;
        k = k - p * tmp;
        int i, cnt = 0;

        for ( i = 1; i <= m; i++ )
        {
            if ( ! pri[i] ) cnt++;
            if ( cnt == k ) break;
        }
        printf("%d\n",i+m*tmp);
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值