結果
| 問題 | No.981 一般冪乗根 |
| ユーザー |
 37zigen
|
| 提出日時 | 2020-02-01 01:21:17 |
| 言語 | Java21 (openjdk 21) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,942 bytes |
| コンパイル時間 | 2,456 ms |
| コンパイル使用メモリ | 79,080 KB |
| 実行使用メモリ | 61,144 KB |
| 最終ジャッジ日時 | 2024-10-09 14:01:19 |
| 合計ジャッジ時間 | 18,049 ms |
|
ジャッジサーバーID （参考情報） |
judge5 / judge3 |
（要ログイン）
テストケース
テストケース表示| 入力 | 結果 | 実行時間 実行使用メモリ |
|---|---|---|
| testcase_00 | TLE | - |
| testcase_01 | -- | - |
| testcase_02 | -- | - |
| testcase_03 | -- | - |
| testcase_04 | -- | - |
| testcase_05 | -- | - |
| testcase_06 | -- | - |
| testcase_07 | -- | - |
| testcase_08 | -- | - |
| testcase_09 | -- | - |
| testcase_10 | -- | - |
| testcase_11 | -- | - |
| testcase_12 | -- | - |
| testcase_13 | -- | - |
| testcase_14 | -- | - |
| testcase_15 | -- | - |
| testcase_16 | -- | - |
| testcase_17 | -- | - |
| testcase_18 | -- | - |
| testcase_19 | -- | - |
| testcase_20 | -- | - |
| testcase_21 | -- | - |
| testcase_22 | -- | - |
| testcase_23 | -- | - |
| testcase_24 | -- | - |
| testcase_25 | -- | - |
| testcase_26 | -- | - |
| testcase_27 | -- | - |
| testcase_28 | -- | - |
| evil_60bit1.txt | -- | - |
| evil_60bit2.txt | -- | - |
| evil_60bit3.txt | -- | - |
| evil_hack | -- | - |
| evil_hard_random | -- | - |
| evil_hard_safeprime.txt | -- | - |
| evil_hard_tonelli0 | -- | - |
| evil_hard_tonelli1 | -- | - |
| evil_hard_tonelli2 | -- | - |
| evil_hard_tonelli3 | -- | - |
| evil_sefeprime1.txt | -- | - |
| evil_sefeprime2.txt | -- | - |
| evil_sefeprime3.txt | -- | - |
| evil_tonelli1.txt | -- | - |
| evil_tonelli2.txt | -- | - |
ソースコード
import java.io.BufferedReader;
import java.io.File;
import java.io.FileNotFoundException;
import java.io.FileReader;
import java.io.IOException;
import java.math.BigInteger;
import java.util.Arrays;
import java.util.Random;
import java.util.Scanner;
class Main {
public static void main(String[] args) throws Exception {
new Main().run();
}
void run() {
Scanner sc = new Scanner(System.in);
int T=sc.nextInt();
for(int i=0;i<T;++i){
long p=sc.nextLong();
long k=sc.nextLong();
long a=sc.nextLong();
long root = kth_root(a,k,p);
System.out.println(root);
}
}
long kth_root(long a,long k,long p){
long g=gcd(k,p-1);
if(pow(a,(p-1)/g,p)!=1)return -1;
a=pow(a,inv(k/g,(p-1)/g),p);
for(long div=2;div*div<=g;++div) {
int sz=0;
while(g%div==0) {
g/=div;
++sz;
}
if(sz>0) {
long b=peth_root(a, pow(div,sz,p), p);
a=b;
}
}
if(g>1)a=peth_root(a,g,p);
return a;
}
long peth_root(long a,long k,long p) {
long sz=p-1;
int s=0;
while(sz%k==0){sz/=k;++s;}
long q1=sz;
long q2=(p-1)/sz;
long ans=pow(a,((k-1)*inv(q1,k)%k*q1+1)/k,p);
long c=2;
while(pow(c,(p-1)/k,p)==1)++c;
c=pow(c,q1,p);
while(true){
long err=pow(ans,k,p)*inv(a,p)%p;
int bit=s-cnt(err,k,p);
if(bit==s)break;
ans=ans*pow(c,pow(k,bit-1,p),p)%p;
}
return ans;
}
int cnt(long a,long base,long p){
int ret=0;
while(a!=1){
a=pow(a,base,p);
++ret;
}
return ret;
}
long inv(long a,long p){
a%=p;
long u=1,v=0;
long b=p;
while(b>0) {
long q=a/b;
a%=b;
u-=v*q%p;
u%=p;
{
u^=v;v^=u;u^=v;
a^=b;b^=a;a^=b;
}
}
return (u%p+p)%p;
}
long pow(long a,long n,long p){
n%=p-1;
long r=1;
for(;n>0;n>>=1,a=a*a%p)if(n%2==1)r=r*a%p;
return r;
}
long gcd(long a,long b) {
return a==0?b:gcd(b%a,a);
}
static void tr(Object... objects) {
System.out.println(Arrays.deepToString(objects));
}
}