結果

問題 No.187 中華風 (Hard)
ユーザー timi
提出日時 2020-12-21 11:53:28
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
RE  
実行時間 -
コード長 665 bytes
コンパイル時間 576 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 11,008 KB
最終ジャッジ日時 2024-09-21 12:50:20
合計ジャッジ時間 2,395 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 5 WA * 10 RE * 10
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
import math

def ext_gcd(p,q):
	if q==0:
		return (p,1,0)
	g,y,x = ext_gcd(q,p%q);
	y -= p/q*x
	return (g,x,y)

def cmt(a1,mo1,a2,mo2):
	g,x,y=ext_gcd(mo1,mo2)
	a1%=mo1
	a2%=mo2
	if a1%g != a2%g:
		return (-1,0)
	lcm=mo1*(mo2/g)
	v=a1+((a2-a1)%lcm)*x%lcm*(mo1/g)
	return (((v%lcm)+lcm) % lcm,lcm)

N=int(input())
P=[]
for i in range(N):
  x,y=map(int, input().split())
  P.append((x,y))
T=[P[0][0],P[0][1]]
for i in range(1,N):
  T2=cmt(T[0],T[1],P[i][0],P[i][1])
  #print(T[0],T[1],P[i][0],P[i][1])
  #print(T2,'modと最小公倍数')
  if T2[0]==-1:
    print(-1)
    exit()
  T=[int(T2[0]),int(T2[1])]

if T[0]==0:
  T[0]+=T[1]

print(T[0]%(10**9+7))
0