結果

問題 No.187 中華風 (Hard)
ユーザー ctyl_0ctyl_0
提出日時 2015-12-04 14:11:54
言語 PyPy2
(7.3.15)
結果
AC  
実行時間 1,301 ms / 3,000 ms
コード長 2,843 bytes
コンパイル時間 1,377 ms
コンパイル使用メモリ 77,232 KB
実行使用メモリ 81,976 KB
最終ジャッジ日時 2024-09-14 12:55:30
合計ジャッジ時間 20,262 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge6
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 292 ms
78,876 KB
testcase_01 AC 290 ms
79,512 KB
testcase_02 AC 1,144 ms
81,256 KB
testcase_03 AC 1,117 ms
80,928 KB
testcase_04 AC 1,248 ms
80,444 KB
testcase_05 AC 1,251 ms
80,652 KB
testcase_06 AC 1,301 ms
81,428 KB
testcase_07 AC 1,274 ms
80,960 KB
testcase_08 AC 901 ms
80,948 KB
testcase_09 AC 990 ms
81,976 KB
testcase_10 AC 984 ms
81,464 KB
testcase_11 AC 1,261 ms
80,528 KB
testcase_12 AC 1,258 ms
80,556 KB
testcase_13 AC 319 ms
79,480 KB
testcase_14 AC 337 ms
80,016 KB
testcase_15 AC 291 ms
79,676 KB
testcase_16 AC 294 ms
79,520 KB
testcase_17 AC 87 ms
78,108 KB
testcase_18 AC 260 ms
78,884 KB
testcase_19 AC 86 ms
78,332 KB
testcase_20 AC 1,031 ms
80,920 KB
testcase_21 AC 86 ms
78,216 KB
testcase_22 AC 1,258 ms
80,648 KB
testcase_23 AC 88 ms
78,764 KB
testcase_24 AC 86 ms
78,356 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def pow_mod(a, b, m):
    ret = 1
    while b > 0:
        if b & 1:
            ret = ret * a % m
        a = a * a % m
        b >>= 1
    return ret

def gcd(p, q):
    if q == 0:
        return p
    return gcd(q, p % q)

def ext_gcd(p, q): # return gcd, (x, y) which satisfies px + qy = 1
    if q == 0:
        return (p, 1, 0)
    g, y, x = ext_gcd(q, p % q)
    return (g, x, y - (p // q) * x)
    
def inv(p, q): # return p^(-1) mod q when (p, q) is co-prime
    g, x, y = ext_gcd(p, q)
    x %= q
    if x < 0:
        x += q
    return x
    
def garner(X, M, op = 1): # C: (X, Y) M: modulo
    # list up primes lower than 40000
    prime_sup = 40000
    primes = []
    is_prime = [0] * prime_sup
    
    for i in xrange(2, prime_sup):
        if is_prime[i] == 0:
            primes.append([i, 0])
            j = i ** 2
            while j < prime_sup:
                is_prime[j] = i
                j += i
    
    for i in xrange(len(X)):
        y = X[i][1]
        for j in xrange(len(primes)):
            cnt = 0
            while (y % primes[j][0] == 0):
                y /= primes[j][0]
                cnt += 1
            primes[j][1] = max(primes[j][1], cnt)
        if y > 1:
            primes.append([y, 1])
    
    for prime in primes[:][:]:
        if prime[1] == 0:
            primes.remove(prime)
    
    # op = 0
    if op == 0:
        ret = int(1)
        for i in xrange(len(primes)):
            ret *= pow_mod(primes[i][0], primes[i][1], M)
            ret %= M
        return ret
        
    # Does the equation have solutoion?    
    
    for i in range(len(X)):
        for j in range(i + 1, len(X)):
            g = gcd(X[i][1], X[j][1])
            if X[i][0] % g != X[j][0] % g:
                return -1
        
    # reguralization
    
    prime_check = [0] * len(primes)
    
    for i in xrange(len(X)):
        for j in xrange(len(primes)):
            if X[i][1] % pow_mod(primes[j][0], primes[j][1], M) == 0 and prime_check[j] == 0:
                prime_check[j] = 1
            else:
                while X[i][1] % primes[j][0] == 0:
                    X[i][1] /= primes[j][0]
                    X[i][0] %= X[i][1]
    
    X.sort(key = lambda x: (x[1]))
    
    # garner's algorithm
    
    for i in xrange(len(X)):
        for j in xrange(i):
            k = X[i][0] - X[j][0]
            if k < 0:
                k += X[i][1]
            X[i][0] = k * inv(X[j][1], X[i][1]) % X[i][1]
    
    ret = int(0)
    
    for i in reversed(xrange(len(X))):
        ret = (ret * X[i][1] + X[i][0]) % M
        
        
    return ret
    

N = int(input())
Z = []
zero = True
mod = int(1e9 + 7)

for i in xrange(N):
    Z.append(map(int, raw_input().strip().split(" ")))
    zero &= (Z[i][0] == 0)

if zero == True:
    print garner(Z, mod, 0)
else:
    print garner(Z, mod, 1)
0