码迷,mamicode.com
首页 > 其他好文 > 详细

hdu 3589(二次剩余+雅可比符号)

时间:2016-07-03 20:06:56      阅读:146      评论:0      收藏:0      [点我收藏+]

标签:

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3589

题意:就是一个裸的雅可比符号:

具体参考:百度百科  百度百科

代码如下:

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
using namespace std;
int a[10];
int fac[100500],sum[100500];
int cnt;
void fen(int x){
    for(int i=2;i*i<=x;i++){
        if(x%i==0){
            fac[cnt]=i;
            while(x%i==0){
                x/=i;
                sum[cnt]++;
            }
            cnt++;
        }
    }
    if(x>0)fac[cnt]=x,sum[cnt++]++;
}
long long quick_mod(int a,int b,int p){
    long long ans=1;
    while(b){
        if(b&1)ans=ans*a%p;
        b>>=1;
        a=a*a%p;
    }
    return ans;
}
int main(){
    int a,n;
    while(scanf("%d%d",&a,&n)!=EOF){
        memset(sum,0,sizeof(sum));
        cnt=0;
        fen(n);
        int ans=1;
        for(int i=0;i<cnt;i++){
            int tmp=(int)quick_mod(a,(fac[i]-1)/2,fac[i]);
            if(tmp>1)tmp=-1;
            if(tmp==-1&&sum[i]%2==0)
            tmp=1;
            ans*=tmp;
        }
        cout<<ans<<endl;
    }
    return 0;
}


hdu 3589(二次剩余+雅可比符号)

标签:

原文地址:http://blog.csdn.net/qq_27599517/article/details/51792790

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!