結果
問題 |
No.2125 Inverse Sum
|
ユーザー |
![]() |
提出日時 | 2021-05-13 16:08:41 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 111 ms / 2,000 ms |
コード長 | 848 bytes |
コンパイル時間 | 151 ms |
コンパイル使用メモリ | 82,984 KB |
実行使用メモリ | 93,036 KB |
最終ジャッジ日時 | 2024-09-20 01:29:35 |
合計ジャッジ時間 | 2,827 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
def fact(n): #nの素因数分解 arr=[] temp=n for i in range(2,int(-(-n**0.5//1))+1): if temp%i==0: cnt=0 while temp%i==0: cnt+=1 temp//=i arr.append([i, cnt]) if temp!=1: arr.append([temp, 1]) if arr==[]: arr.append([n, 1]) return arr def div(n): #n^2の約数列挙 arr=[1] for i,j in fact(n): tmp=arr.copy() arr=[] for a in range(2*j+1): for k in tmp: arr.append(k*pow(i,a)) return sorted(arr) def solve(p,q): #(pa-q)*(pb-q)=q^2 を解く ans=[] check=set() for i in div(q): j=q**2//i if (i+q)%p==0: if (i+q)//p not in check: ans.append([(i+q)//p,(j+q)//p]) check.add((i+q)//p) print(len(ans)) for i in ans: print(*i) import math p,q=map(int,input().split()) solve(p//math.gcd(p,q),q//math.gcd(p,q))