大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。
Jetbrains全家桶1年46,售后保障稳定
今天学了数论。。。求原根真的好暴力
设模数为p,
我们把 p − 1 p-1 p−1分解质因数, p 1 , p 2 … p m p_1,p_2 \ldots p_m p1,p2…pm对于每一个 2 ≤ i ≤ p − 1 2 \leq i \leq p-1 2≤i≤p−1 ,判断 i p − 1 p j % p i^{p-1 \over p_j} \%p ipjp−1%p是否为1,如果是,那么这个数就不是原根,否则就是
AC Code
#include<cstdio>
#include<iostream>
#include<cmath>
using namespace std;
int zys[50];
int p;
int o;
int n;
int ksm(int a,int c,int b)
{
long long re=1;
long long t=a;
while(c)
{
if(c&1)re=re%b*t%b;
t=t%b*t%b;
c>>=1;
}
return re;
}
int pd(int x)
{
for(int i=1;i<=o;i++)
{
if(ksm(x,(p-1)/zys[i],p)%p==1)return 0;
}
return 1;
}
int main()
{
scanf("%d",&p);
int pp=p;
p--;
for(int i=2;i<=sqrt(p);i++)
{
if(p==1)break;
if(p%i==0)
{
o++;
zys[o]=i;
while(p%i==0)
{
p/=i;
}
}
}
if(p!=1)
{
o++;
zys[o]=p;
}
p=pp;
for(int i=2;i<p;i++)
{
if(pd(i))
{
printf("%d",i);
break;
}
}
return 0;
}
发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/215435.html原文链接:https://javaforall.cn
【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛
【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...