結果
問題 | No.990 N×Mマス計算(Kの倍数) |
ユーザー | kohei2019 |
提出日時 | 2021-02-03 23:11:35 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 229 ms / 2,000 ms |
コード長 | 4,565 bytes |
コンパイル時間 | 144 ms |
コンパイル使用メモリ | 82,244 KB |
実行使用メモリ | 93,932 KB |
最終ジャッジ日時 | 2024-06-30 06:12:18 |
合計ジャッジ時間 | 3,116 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 43 ms
55,384 KB |
testcase_01 | AC | 42 ms
54,712 KB |
testcase_02 | AC | 44 ms
55,600 KB |
testcase_03 | AC | 43 ms
55,328 KB |
testcase_04 | AC | 44 ms
55,700 KB |
testcase_05 | AC | 43 ms
55,684 KB |
testcase_06 | AC | 43 ms
55,820 KB |
testcase_07 | AC | 45 ms
56,724 KB |
testcase_08 | AC | 43 ms
55,412 KB |
testcase_09 | AC | 43 ms
55,108 KB |
testcase_10 | AC | 113 ms
83,372 KB |
testcase_11 | AC | 139 ms
87,276 KB |
testcase_12 | AC | 229 ms
91,596 KB |
testcase_13 | AC | 107 ms
82,644 KB |
testcase_14 | AC | 124 ms
84,112 KB |
testcase_15 | AC | 110 ms
78,596 KB |
testcase_16 | AC | 123 ms
82,960 KB |
testcase_17 | AC | 100 ms
78,660 KB |
testcase_18 | AC | 226 ms
91,652 KB |
testcase_19 | AC | 139 ms
82,052 KB |
testcase_20 | AC | 155 ms
93,932 KB |
ソースコード
import math import sys sys.setrecursionlimit(10**7) #競技プログラミング対整数問題のライブラリーです class integerlib(): def __init__(self): pass def primeset(self,N): #N以下の素数をsetで求める.エラトステネスの篩O(√Nlog(N)) lsx = [1]*(N+1) for i in range(2,int(-(-N**0.5//1))+1): if lsx[i] == 1: for j in range(i,N//i+1): lsx[j*i] = 0 setprime = set() for i in range(2,N+1): if lsx[i] == 1: setprime.add(i) return setprime def defprime(self,N):#素数かどうかの判定、エラトステネスの篩O(√Nlog(N)) return N in self.primeset(N) def gcd(self,ls):#最大公約数 ls = list(ls) ans = 0 for i in ls: ans = math.gcd(ans,i) return ans def lmc(self,ls):#最小公倍数 ls = list(ls) ans = self.gcd(ls) for i in ls: ans = self.lmcsub(ans,i) return ans def lmcsub(self,a,b): gcd = math.gcd(a,b) lmc = (a*b)//gcd return lmc def factorization(self,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 factorizationset(self,N):#素因数分解√N,含まれている素因数の種類 if N == 1: return set() ls = self.factorization(N) setf = set() for j in ls: setf.add(j[0]) return setf def divisorsnum(self,N):#約数の個数 ls = [] for i in self.factorization(N): ls.append(i[1]) d = 1 for i in ls: d *= i+1 return d def Eulerfunc(self,N):#オイラー関数正の整数Nが与えられる。1,2,…,Nのうち、Nと互いに素であるものの個数を求めよ。 ls = list(self.factorizationset(N)) ls2 = [N] for i in ls: ls2.append(ls2[-1]-ls2[-1]//i) return ls2[-1] def make_divisors(self,N):#約数列挙O(√N) lower_divisors , upper_divisors = [], [] i = 1 while i*i <= N: if N % i == 0: lower_divisors.append(i) if i != N // i: upper_divisors.append(N//i) i += 1 return lower_divisors + upper_divisors[::-1] def invmod(self,a,mod):#mod逆元 if a == 0: return 0 if a == 1: return 1 return (-self.invmod(mod % a, mod) * (mod // a)) % mod def cmbmod(self,n, r, mod):#nCr % mod inv = [0,1] for i in range(2, n + 1): inv.append((-inv[mod % i] * (mod // i)) % mod) cmd = 1 for i in range(1,min(r,n-r)+1): cmd = (cmd*(n-i+1)*inv[i])%mod return cmd def permmod(self,n, r, mod):#nPr % mod perm = 1 for i in range(n,r-1,-1): perm = (perm*i)%mod return perm def modPow(self,a,n,mod):#繰り返し二乗法 a**n % mod if n==0: return 1 if n==1: return a%mod if n % 2 == 1: return (a*self.modPow(a,n-1,mod)) % mod t = self.modPow(a,n//2,mod) return (t*t)%mod import collections N,M,K = map(int,input().split()) lsB = input().split() f = lsB.pop(0) lsB = [int(i) for i in lsB] lsB.sort() lsA = [int(input()) for i in range(N)] IT = integerlib() if f == '*': lsgcdb = [] for i in range(M): lsgcdb.append(math.gcd(lsB[i],K)) Cb = collections.Counter(lsgcdb) lsrma = [] for i in range(N): a = lsA[i] gcda = math.gcd(a,K) rm = K//gcda lsrma.append(rm) Ca = collections.Counter(lsrma) cnt = 0 for keya,valuea in Ca.items(): for key,value in Cb.items(): if key%keya == 0: cnt += value*valuea print(cnt) else: for i in range(N): lsA[i] %= K for i in range(M): lsB[i] %= K Cb = collections.Counter(lsB) cnt = 0 for i in range(N): if lsA[i] == 0: cnt += Cb[0] else: cnt += Cb[K-lsA[i]] print(cnt)