hdu 2879 HeHe 欧筛&&积性函数

In the equation X^2≡X(mod N) where x∈[0,N-1], we define He[N] as the number of solutions. 
And furthermore, define HeHe[N]=He[1]*……*He[N] 
Now here is the problem, write a program, output HeHe[N] modulo M for a given pair N, M.

Input

First line: an integer t, representing t test cases. 
Each test case contains two numbers N (1<=N<=10^7) and M (0<M<=10^9) separated by a space.

Output

For each test case, output one line, including one integer: HeHe[N] mod m.

Sample Input

1
2 3

Sample Output

2

积性函数定义：
对于正整数n的一个算术函数 f(n)，若f(1)=1，且当a,b互质时f(ab)=f(a)f(b)，在数论上就称它为积性函数。若对于某积性函数 f(n) ，就算a, b不互质，也有f(ab)=f(a)f(b)，则称它为完全积性的
常见积性函数： 
φ(n) －欧拉函数，计算与n互质的正整数之数目　 
μ(n) －莫比乌斯函数，关于非平方数的质因子数目 
gcd(n,k) －最大公因子，当k固定的情况　 
d(n) －n的正因子数目 
σ(n) －n的所有正因子之和　 
σk(n) － 因子函数，n的所有正因子的k次幂之和，当中k可为任何复数。　 
Idk(n) －幂函数，对于任何复数、实数k，定义为Idk(n) = n^k （完全积性）　 
λ(n) －刘维尔函数，关于能整除n的质因子的数目　 
γ(n)，定义为γ(n)=(-1)^ω(n)，在此加性函数ω(n)是不同能整除n的质数的数目
 

题解：通过证明可得，he是个积性函数，然后he[p]=2，这样求hehe就可以找一下1-n质因子的总个数k

然后结果是2^k

#include <bits/stdc++.h>
#define maxn 10000000+5
using namespace std;
typedef long long int ll;
int tot,prime[maxn],is_not_pr[maxn];
void get_list(){
    tot=0;
    memset(is_not_pr,0,sizeof(is_not_pr));
    for(int i=2;i<maxn;i++){
        if(!is_not_pr[i]) prime[++tot]=i;
        for(int j=1;j<=tot&&i*prime[j]<=maxn;j++){
            is_not_pr[i*prime[j]]=1;
                //合数标为1，同时，prime[j]是合数i*prime[j]的最小素因子
            if(i%prime[j]==0) break;
                    //即比一个合数大的质数和该合数的乘积可用一个更大的合数和比其小的质数相乘得到
        }
    }
}
ll q_pow(ll a,ll b,ll mod){
    ll ans = 1;
    while(b){
        if(b & 1)
            ans = ans * a % mod;
        b >>= 1;
        a = a * a % mod;
    }
    return ans;
}
int main(){
    get_list();
    /*for(int i=1;i<=2;i++){
        cout<<prime[i]<<endl;
    }*/
    ll n,m;
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%lld%lld",&n,&m);
        ll coun=0;
        for(int i=1;i<=tot&&prime[i]<=n;i++){
            coun+=n/prime[i];
        }
        ll ans=q_pow(2,coun,m);
        printf("%lld\n",ans);
    }
    return 0;
}