URAL 1226 Little Chu(最大原根)

原创 2016年05月31日 01:21:35

题目链接:
URAL 1226 Little Chu*
题意:
找到最大的k,使得k, k^2, k^3…k^x 模n的结果各不相同且x尽可能大。
分析:
求最大原根。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <climits>
#include <cmath>
#include <ctime>
#include <cassert>
#include <bitset>
#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
using namespace std;
typedef long long ll;
const int MAX_N = 70000;

int T, n, prime_cnt, factor_cnt;
int prime[MAX_N], factor[MAX_N];
bitset<MAX_N> bs;

void GetPrime()
{
    bs.set();
    prime_cnt = 0;
    for(int i = 2; i < MAX_N; i++){
        if(bs[i] == 1) { prime[prime_cnt++] = i; }
        for(int j = 0; j < prime_cnt && i * prime[j] < MAX_N; j++){
            bs[i * prime[j]] = 0;
            if(i % prime[j] == 0) break;
        }
    }
}

void GetFactor(int x)
{
    factor_cnt = 0;
    for(int i = 0; i < prime_cnt && prime[i] * prime[i] <= x; i++){
        if(x % prime[i] == 0){
            factor[factor_cnt++] = prime[i];
            while(x % prime[i] == 0){
                x /= prime[i];
            }
        }
    }
    if(x > 1){
        factor[factor_cnt++] = x;
    }
}

int quick_pow(int a, int b, int m)
{
    ll res = 1, tmp = (ll)a, mod = (ll)m;
    while(b){
        if(b & 1) res = res * tmp % mod;
        tmp = tmp * tmp % mod;
        b >>= 1;
    }
    return (int)res;
}

int main()
{
    GetPrime();
    scanf("%d", &T);
    while(T--){
        scanf("%d", &n);
        GetFactor(n - 1);
        int ans ;
        for(int i = n - 1; i >= 2; i--){
            int flag = true;
            for(int j = 0; j < factor_cnt; j++){
                if(quick_pow(i, (n - 1) / factor[j], n) == 1){
                    flag = false;
                    //printf("i = %d factor[j] = %d\n", i, factor[j]);
                    break;
                }
            }
            if(flag){
                ans = i;
                break;
            }
        }
        printf("%d\n", ans);
    }
    return 0;
}
版权声明:缥缈玉京人,想语然、京兆眉妩。

相关文章推荐

ural 1099. Work Scheduling 一般图最大匹配

There is certain amount of night guards that are available to protect the local junkyard from possib...

ural 1146 Maximum Sum 最大连续和

题意:给个n*n的矩阵,所有子矩阵中 ,和最大等于多少。 做法: 首先要理解一个O(n)的算法 。 给一个数组 求连续和的最大值。 可以用一个sum来从下标0开始计算和,不断取最大值。当加和小...

原根和指数表

  • 2015-07-21 10:16
  • 2.02MB
  • 下载

阶与原根学习笔记

阶的定义: 正整数(a,m)=1,则最小的正整数r使得a^r=1(mod m) 为a模m的阶。 阶的性质: r | φ(m) 。(可用反证法,假设不整除,则……) 求阶: 给定互质的a,m,求a模m的...

[LOJ #6156][找原根][FFT || NTT] A*B Problem

题目要求对于所有0≤k<m0\leq k<m的kk都输出解,可以考虑anskans_k的生成函数。 ansk=∑i×j=kai∗ajans_k=\sum_{i\times j=k} a_i*a_j ...

离散对数和原根

原根与阶阶定义设 (a,m)=1(a,m)=1, 满足 ax≡1(mod m)a^x \equiv 1(mod\ m) 的最小的 xx,称为a对m的阶,记为 ordm(a)ord_{m}(a) 当...

原根产生算法

  • 2015-05-03 10:34
  • 1003B
  • 下载

51nod 1135 原根

#include using namespace std; const int MAXN=100100; int prime[MAXN],npri[MAXN],dive[MAXN]; int ps...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)