HDU 3589 Jacobi symbol
彻底对数学绝望了
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 | #include <cstdio> #include <cmath> int flag[1005],p[500],a; int d[100]; int init( int s) { int len=0,tmp,h= sqrt (s+0.5); for ( int i=0;p[i]<=h;i++) if (s%p[i]==0) { if (a%p[i]==0) return -1; while (s%p[i]==0) { d[len++]=p[i]; s/=p[i]; } if (s==1) return len; } if (s>1) { if (a%s==0) return -1; d[len++]=s; } return len; } int getans( int x, int s) { int tmp=s/2; __int64 ans=1,b=x; while (tmp>0) { if (tmp&1)ans=ans*b%s; b=b*b%s; tmp/=2; } if (ans!=1) return -1; return 1; } int main() { int i,j,k=0,n; for (i=2;i<1000;i++) { if (!flag[i]) { p[k++]=i; for (j=i*i;j<1000;j+=i) flag[j]=1; } } while ( scanf ( "%d%d" ,&a,&n)!=-1) { int len=init(n); if (len==-1) { printf ( "0\n" ); continue ; } int ans=1; for (i=0;i<len;i++) ans*=getans(a,d[i]); printf ( "%d\n" ,ans); } return 0; } |
愿你出走半生,归来仍是少年
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步