ACdreamer

欢迎关注知乎:https://www.zhihu.com/people/acdreamers/activities

HDU2879(积性函数)

题目:http://acm.hdu.edu.cn/showproblem.php?pid=2879



#include <iostream>
#include <string.h>
#include <stdio.h>

using namespace std;
typedef long long LL;
const int N = 10000005;
const int M = 1000005;

bool prime[N];
int p[N];
int k;

void isprime()
{
    k = 0;
    memset(prime,true,sizeof(prime));
    for(int i=2;i<N;i++)
    {
        if(prime[i])
        {
            p[k++] = i;
            for(int j=i+i;j<N;j+=i)
                prime[j] = false;
        }
    }
}

LL quick_mod(LL a,LL b,LL m)
{
    LL ans = 1;
    a %= m;
    while(b)
    {
        if(b&1)
        {
            ans = ans * a % m;
            b--;
        }
        b >>= 1;
        a = a * a % m;
    }
    return ans;
}

int main()
{
    int n,m,t;
    isprime();
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        LL cnt = 0;
        for(int i=0;i < k;i++)
        {
            if(p[i] > n) break;
            cnt += n / p[i];
        }
        LL ans = quick_mod(2,cnt,m);
        printf("%I64d\n",ans);
    }
    return 0;
}



阅读更多
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/ACdreamers/article/details/17101861
个人分类: 数论
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

不良信息举报

HDU2879(积性函数)

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭