a乘b后取mod p

对a和b拆解,

例如a=7

a=2^{0}+2^{1}+2^{2}

例如b=5

b=2^{0}+2^{2}

相乘时

a x b=(2^{0}+2^{1}+2^{2}) \times(2^{0}+2^{2})=(2^{0}\times 2^{0}+2^{0}\times 2^{2})+(2^{1}\times 2^{0}+2^{1}\times 2^{2})+(2^{2}\times 2^{0}+2^{2}\times 2^{2})

=2^{0}+2^{2}+2^{1}+2^{3}+2^{4}+2^{4}

只需要对各个项的上标进行计算即可;

a\times b mod(p)=2^{0}mod(p)+2^{2}mod(p)+2^{1}mod(p)+2^{3}mod(p)+2^{4}mod(p)+2^{4}mod(p)

计算2^{n}可能会很大,但是n不会特别大,所以可以转变成n个2连乘,每次乘2都mod p就不会溢出。

#include"stdio.h"
long long power(long long b,long long p);
int main()
{
    long long a=0L;
    long long b=0L;
    long long p=0L;
    long long ans=0L;
    long long t=1L;
    int i=0,j,k=0,k1=0;
    int m1[100];//存b的二进制的每一项2的上标 比如7的2进制(111=1*2^0+1*2^1+1*2^2)上标0 1 2 即2的0次方,2的1次方,2的2二次方
    int m2[100];//存a的二进制的每一项2的上标

    scanf("%lld",&a);
    scanf("%lld",&b);
    scanf("%lld",&p);


    while(b>0)
    {
        if(b&1)
        {
            m1[k++]=i;

        }
        i++;

        b=b>>1;
    }
    /**
    for(j=0;j<k;j++)
    {
        printf("%d ",m1[j]);

    }
    printf("\n");

*/
    t=1L;
    ans=0L;
    i=0;
    while(a>0)
    {
        if(a&1)
        {

            m2[k1++]=i;
        }
        i++;
        a=a>>1;
        //printf("%lld    %d\n",ans,m2[i-1]);
    }
    /*
    for(j=0;j<k1;j++)
    {
        printf("%d ",m2[j]);

    }
        printf("\n");
        */
    ans=0L;
    for(i=0;i<k;i++)
    {
        for(j=0;j<k1;j++)
        {
            ans+=power(m1[i]+m2[j],p);
            ans=ans%p;

        }

    }
    printf("%lld  ",ans);
    //printf("\n");
   // printf("%lld ",power(2,10));
}
long long power(long long b,long long p)
{

    int i=0;
    long long a=1;
    for(i;i<b;i++)
    {
        a=a*2%p;
    }
    return a;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值