結果
問題 | No.2125 Inverse Sum |
ユーザー |
👑 ![]() |
提出日時 | 2022-11-18 21:38:44 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 94 ms / 2,000 ms |
コード長 | 1,362 bytes |
コンパイル時間 | 141 ms |
コンパイル使用メモリ | 81,816 KB |
実行使用メモリ | 82,248 KB |
最終ジャッジ日時 | 2024-09-20 02:02:24 |
合計ジャッジ時間 | 2,595 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
def Prime_Factorization(N): if N==0: return [[0,1]] if N<0: R=[[-1,1]] else: R=[] N=abs(N) if N&1==0: C=0 while N&1==0: N>>=1 C+=1 R.append([2,C]) if N%3==0: C=0 while N%3==0: N//=3 C+=1 R.append([3,C]) k=5 Flag=0 while k*k<=N: if N%k==0: C=0 while N%k==0: C+=1 N//=k R.append([k,C]) k+=2+2*Flag Flag^=1 if N!=1: R.append([N,1]) return R #素因数分解の結果から, 約数を全て求める. def Divisors_from_Prime_Factor(P,sorting=False): X=[1] for p,e in P: q=1 n=len(X) for _ in range(e): q*=p for j in range(n): X.append(X[j]*q) if sorting: X.sort() return X def solve(): from math import gcd P,Q=map(int,input().split()) g=gcd(P,Q) P//=g; Q//=g E=[(p,2*e) for p,e in Prime_Factorization(Q)] Ans=[] for a in Divisors_from_Prime_Factor(E, True): if (Q+a)%P==0 and (Q+Q*Q//a)%P==0: Ans.append(((Q+a)//P, (Q+Q*Q//a)//P)) print(len(Ans)) for N,M in Ans: print(N,M) #================================================== solve()