D.Easy problem

莫比乌斯反演
欧拉降幂

#include <bits/stdc++.h>
using namespace std;
typedef long long ll ;
const int mod = 59964251 ;
const int N = 1e5 + 10 ;
int prime[N] , vis[N] , mu[N] , tot , sum[N] ;
ll qmi(ll a , ll b)
{
	ll res = 1 ; 
	while(b)
	{
		if(b & 1) res = res * a % mod ;
		a = a * a % mod ;
		b >>= 1 ;
	}
	return res ;
}
void get()
{
	mu[1] = 1 ;
	for(int i = 2; i < N ;i ++)
	 {
	 	if(!vis[i]) prime[++ tot] = i , mu[i] = -1 ;
	 	for(int j = 1; j <= tot && i * prime[j] < N ;j ++)
	 	 {
	 	 	vis[i * prime[j]] = 1 ;
	 	 	if(i % prime[j] == 0) break ;
	 	 	mu[i * prime[j]] = -mu[i] ; 
		  }
	 }
}
int get_phi(int n)
{
	ll res = n ;
	for(int i = 2; i * i <= n ;i ++)
	 {
	 	if(n % i == 0)
	 	 {
	 	 	res = res / (i) * (i - 1) ;
	 	 	while(n % i == 0)
	 	 	 n /= i ;
		 }
		if(n == 1) break ;
	 }
	 if(n > 1) res = res / n * (n - 1) ;
	 return res ; 
}
int main()
{
	int T ;
	get() ;

	scanf("%d" , &T) ;
	while(T --)
	{
		char s[110000] ;
		int m , d , k ; 
		
		scanf("%s" , s) ;
	
		scanf("%d%d%d" , &m , &d , &k) ;
		int len = strlen(s) ;
		ll n = 0 ; 
		ll phi = get_phi(mod) ;
		
		for(int i = 0 ;i < len ;i ++)
		 {
		 	n = n * 10 + s[i] - '0' ;
		 	n %= phi ;
		 }
	    m /=d ;
	    for(int i = 1; i <= m ;i ++)
	     sum[i] = (sum[i - 1] + qmi(i , k)) % mod ;
	    ll ans = 0 ;
	    for(int i = 1; i <=  m ;i ++)
	     
	     	ans = (ans + mu[i] * qmi(i , k * n % phi + phi) % mod * qmi(sum[m / i] , n % phi + phi) % mod + mod) % mod ;
		 
		 ans = (ans * qmi(d , k * n % phi + phi)) % mod ;
		 printf("%d
" , ans) ;
		 
	}
	return 0 ; 
}