2020-2021第二届全国大学生算法设计与编程挑战赛 （秋季正式赛）-题解-10个

#include <cstdio>

const int N = 1e6 + 5;
int tot, p[N];
bool flg[N];

void sieve(int n){
   
	for(int i=2; i<=n; i++){
   
		if(!flg[i]){
   
			p[++tot] = i;
		}
		for(int j=1; j<=tot&&i*p[j]<=n; ++j){
   
			flg[i*p[j]] = 1;
			if(i%p[j]==0){
   
				break;
			}
		}
	}
} 

long long phi(long long x){
   
	long long ans = x;
	for(int i=1; i<=tot&& 1ll*p[i]*p[i]; ++i){
   
		if(x%p[i]){
   
			continue;
		}
		ans = ans/p[i]*(p[i]-1);
		while(x%p[i]==0){
   
			x /= p[i];
		}
	}
	if(x>1){
   
		ans = ans / (x * (x - 1));
	}
	return ans;
}

int main() {
   
	sieve(N - 5);
	long long n, k;
	scanf("%llld%lld", &n, &k);
	k = (k + 1) / 2;
	for(long long i=1; i<=k&&n>1; ++i){
   
		n = phi(n);
	}
	printf("%lld\n", n%1000000007);
	
	return 0;
}

#include <iostream>
#include <cmath>
using namespace std;
int work(int x){
   
	int L = 0, R = -1, lst = -1; //L最高位 R最低位 lst上一位数 
	while(x