結果

問題 No.990 N×Mマス計算（Kの倍数）
コンテスト
ユーザー neterukun
提出日時 2020-02-14 21:54:03
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 1,564 bytes
コンパイル時間 210 ms
コンパイル使用メモリ 82,160 KB
実行使用メモリ 99,176 KB
最終ジャッジ日時 2024-11-16 00:36:44
合計ジャッジ時間 3,190 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 9 WA * 3 RE * 7
権限があれば一括ダウンロードができます

ソースコード

diff # 

def gcd(a: int, b: int) -> int:
    """a, bの最大公約数(greatest common divisor: GCD)を求める
    計算量: O(log(min(a, b)))
    """
    if b == 0:
        return a
    return gcd(b, a%b)


def make_divisors(n):
    """自然数nの約数を列挙したリストを出力する
    計算量: O(sqrt(N))
    入出力例: 12 -> [1, 2, 3, 4, 6, 12]
    """
    divisors = []
    for k in range(1, int(n**0.5) + 1):
        if n % k == 0:
            divisors.append(k)
            if k != n // k:
                divisors.append(n // k)
    divisors = sorted(divisors)
    return divisors

n, m, k = list(map(int, input().split()))
op = list(input().split())
b = [int(op[i]) for i in range(1, m + 1)]
a = [int(input()) for i in range(n)]

ans = 0
if op[0] == "+":
    a = [a[i] % k for i in range(n)]
    b = [b[i] % k for i in range(m)]
    cnt_b = {}
    for num in b:
        if num not in cnt_b:
            cnt_b[num] = 1
        else:
            cnt_b[num] += 1
    for num in range(n):
        if num == 0:
            ans += cnt_b[num]
        if k - num in cnt_b:
            ans += cnt_b[k - num]
else:
    li = make_divisors(k)
    to_ind = {v: i for i, v in enumerate(li)}
    a_li = [0] * len(li)
    b_li = [0] * len(li)
    for num in a:
        tmp = gcd(num, k)
        a_li[to_ind[tmp]] += 1  
    for num in b:
        tmp = gcd(num, k)
        b_li[to_ind[tmp]] += 1 

    for i, num1 in enumerate(li):
        for j, num2 in enumerate(li):
            if (num1 * num2) % k == 0:
                ans += a_li[i] * b_li[j]

print(ans)
0