经验首页 前端设计 程序设计 Java相关 移动开发 数据库/运维 软件/图像 大数据/云计算 其他经验
当前位置:技术经验 » 程序设计 » C++ » 查看文章
51nod 1135 原根(原根)
来源:cnblogs  作者:自为风月马前卒  时间:2018/11/28 9:54:29  对本文有异议

题意

题目链接

Sol

可以证明素数的原根不会超过他的\(\frac{1}{4}\)

那么预处理出\(P - 1\)的所有的质因数\(p_1, p_2 \dots p_k\),暴力判断一下,如果$\exists i, a^{\frac{P - 1}{p_i}} \equiv 1 \pmod {P - 1} $

那么说明\(a\)不是\(P\)的原根,因为根据原根的定义,需要保证\(P-1\)是第一个满足\(a^{P - 1} \equiv 1 \pmod {P - 1}\)的数

  1. #include<cstdio>
  2. using namespace std;
  3. const int MAXN = 1e6 + 10;
  4. int fp(int a, int p, int mod) {
  5. int base = 1;
  6. while(p) {
  7. if(p & 1) base = 1ll * base * a % mod;
  8. a = 1ll * a * a % mod; p >>= 1;
  9. }
  10. return base;
  11. }
  12. int GetG(int x) {
  13. static int q[MAXN]; int tot = 0, tp = x - 1;
  14. for(int i = 2; i * i <= tp; i++) {//这里是i * i
  15. if(!(tp % i)) {
  16. q[++tot] = i;
  17. while(!(tp % i)) tp /= i;
  18. }
  19. }
  20. if(tp > 1) q[++tot] = tp;
  21. for(int i = 2, j; i <= x - 1; i++) {
  22. for(j = 1; j <= tot; j++) if(fp(i, (x - 1) / q[j], x) == 1) break;
  23. if(j == tot + 1) return i;
  24. }
  25. }
  26. int main() {
  27. int P;
  28. scanf("%d", &P);
  29. printf("%d", GetG(P));
  30. return 0;
  31. }
  32. /*
  33. 1000000007
  34. */
 友情链接:直通硅谷  点职佳  北美留学生论坛

本站QQ群:前端 618073944 | Java 606181507 | Python 626812652 | C/C++ 612253063 | 微信 634508462 | 苹果 692586424 | C#/.net 182808419 | PHP 305140648 | 运维 608723728

W3xue 的所有内容仅供测试,对任何法律问题及风险不承担任何责任。通过使用本站内容随之而来的风险与本站无关。
关于我们  |  意见建议  |  捐助我们  |  报错有奖  |  广告合作、友情链接(目前9元/月)请联系QQ:27243702 沸活量
皖ICP备17017327号-2 皖公网安备34020702000426号