Basic Gcd Problem

题目链接

题目要求：得到

i=1->i=x-1,所以gcd(i,x)的值即为x的所有因数，也就是说f(x)=max(c*f(xi)),xi为x的所有因数，我们从x往下推，得到x的因数的函数值，再得到x的因数的因数的函数值…最终状态的函数值为1，也就是说，所谓的函数值其实就是c的次幂，而幂的大小为x到其多次转移得到的因数为质数的过程次数。

所以我们要得到上述过程的最大次数，我们再由下往上推，当t为一素数转移n次得到的，那么f(t)=n,如果t可以转移到n,那么f(n)=f(t)+1.

首先我们要得到小于x的所有素数，我们这里预处理通过欧拉筛来得到所需的素数，时间复杂度为O(n)。

欧拉筛

bool is_prime[N];
int prime[N];
int cnt;

void get_prime(int maxm) {
	is_prime[1] = is_prime[0] = true;
	for (int i = 2; i <= maxm; i++) {
		if (!is_prime[i]) {
			prime[cnt++] = i;
		}
		for (int j = 0; j < cnt && i * prime[j] < maxm; j++) {
			is_prime[i * prime[j]] = true;
			if (i % prime[j]==0) break;
		}
	}
}

转移过程推广获得转移次数

for (int i = 0; i < cnt; i++) {
        fac[prime[i]] = 1;
        int x = prime[i];
        for (int j = 2; j * x <= maxm; j++) {
            fac[j * x] = max(fac[j * x], max(fac[x] + 1,fac[j]+1));//x,j均为x*j的因数
        }
    }

得到x的最大转移次数后，我们所求的答案即为c^fac[x],采用快速幂处理

AC代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
 
using namespace std;
 
const int mod = 1e9 + 7;
 
const int N = 1e6 + 10;
 
typedef long long ll;
 
int n;
 
int prime[N+10];
bool is_prime[N+10];
int cnt;
int fac[N];
 
void get_phi(int N) {
    for (int i = 2; i <= N; i++) {
        if (!is_prime[i]) {
            prime[cnt++] = i;
        }
        for (int j = 0; j < cnt && i * prime[j] <N; j++) {
            is_prime[i * prime[j]] = true;
            if (i % prime[j] == 0) {
                break;
            }
        }
    }
}
 
ll quick_pow(ll a, ll b) {
    ll ans = 1;
    a %= mod;
    while (b) {
        if (b & 1) ans = (ans * a) % mod;
        a = (a * a) % mod;
        b >>= 1;
    }
    return ans;
}
 
int s[N];
int c[N];
 
int main() {
    int t;
    scanf("%d", &t);
    int maxm = -1;
    for (int i = 1; i <= t; i++) {
        scanf("%d %d", &s[i], &c[i]);
        if (s[i] > maxm) maxm = s[i];
    }
    get_phi(maxm+1);
    for (int i = 0; i < cnt; i++) {
        fac[prime[i]] = 1;
        int x = prime[i];
        for (int j = 2; j * x <= maxm; j++) {
            fac[j * x] = max(fac[j * x], max(fac[x] + 1,fac[j]+1));
        }
    }
    for (int i = 1; i <= t; i++) {
        printf("%lld\n", quick_pow((ll)c[i], (ll)fac[s[i]]) % mod);
    }
}