欧拉函数---数论

最新推荐文章于 2024-01-05 23:21:43 发布

谛听-

最新推荐文章于 2024-01-05 23:21:43 发布

阅读量336

点赞数

分类专栏：题

本文链接：https://blog.csdn.net/u012319493/article/details/51436130

版权

题专栏收录该内容

456 篇文章 1 订阅

订阅专栏

时间限制:10000ms
单点时限:1000ms
内存限制:256MB
描述

小Hi和小Ho有时候会用密码写信来互相联系，他们用了一个很大的数当做密钥。小Hi和小Ho约定了一个区间[L,R]，每次小Hi和小Ho会选择其中的一个数作为密钥。

小Hi：小Ho，这次我们选[L,R]中的一个数K。

小Ho：恩，小Hi，这个K是多少啊？

小Hi：这个K嘛，不如这一次小Ho你自己想办法算一算怎么样？我这次选择的K满足这样一个条件：

假设φ(n)表示1..n-1中与n互质的数的个数。对于[L,R]中的任意一个除K以外的整数y，满足φ(K)≤φ(y)且φ(K)=φ(y)时，K

#include "iostream"
#include "algorithm"
#include "vector"
#include "fstream"
#include "string.h"
#include "stdio.h"
using namespace std;

int oula[5000001];
vector<int> prime;
int isPrime[5000001];

int main()
{
    int l, r;
    cin >> l >> r;
    memset(isPrime, 1, sizeof(isPrime));
    isPrime[0] = isPrime[1] = 0;
    int count = 0;
    for(int i=2; i<=r; i++)
    {
        if(isPrime[i])
        {
            prime.push_back(i);
            oula[i] = i-1;
            count++;
        }
        for(int j=0; j<count; j++)
        {
            if(i * prime[j] > r)
                break;
            isPrime[i * prime[j]] = 0;
            if(i % prime[j] == 0)
                oula[i * prime[j]] = oula[i] * prime[j];
            else
                oula[i * prime[j]] = oula[i] * (prime[j] - 1);
        }
    }
    int ans = l;
    for(int k=l+1; k<=r; k++)
        if(oula[k] < oula[ans])
            ans = k;
    cout << ans;
    return 0;
}