結果
問題 | No.2119 一般化百五減算 |
ユーザー |
![]() |
提出日時 | 2022-11-04 22:59:07 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,077 bytes |
コンパイル時間 | 196 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 23,996 KB |
最終ジャッジ日時 | 2024-07-18 20:54:37 |
合計ジャッジ時間 | 9,106 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 22 TLE * 2 -- * 1 |
ソースコード
def inv_gcd(a,b): a=a%b if a==0: return (b,0) s=b;t=a m0=0;m1=1 while(t): u=s//t s-=t*u m0-=m1*u s,t=t,s m0,m1=m1,m0 if m0<0: m0+=b//s return (s,m0) def inv_mod(x,m): assert 1<=m z=inv_gcd(x,m) assert z[0]==1 return z[1] def crt(r,m,maxr): assert len(r)==len(m) n=len(r) r0=0;m0=1 for i in range(n): assert 1<=m[i] r1=r[i]%m[i] m1=m[i] if m0<m1: r0,r1=r1,r0 m0,m1=m1,m0 if (m0%m1==0): if (r0%m1!=r1): return (0,0) continue g,im=inv_gcd(m0,m1) u1=m1//g if ((r1-r0)%g): return (0,0) x=(r1-r0)//g % u1*im%u1 r0+=x*m0 m0*=u1 if r0<0: r0+=m0 if r0>maxr: return (0,0) return (r0,m0) n = int(input()) b = []; c = [] for _ in range(int(input())): bi, ci = map(int, input().split()) b.append(bi); c.append(ci) r, m = crt(c, b, n) print(r if m else "NaN")