結果
問題 | No.171 スワップ文字列(Med) |
ユーザー | rpy3cpp |
提出日時 | 2015-05-19 02:28:12 |
言語 | Python3 (3.11.6 + numpy 1.26.0 + scipy 1.11.3) |
結果 |
AC
|
実行時間 | 20 ms / 1,000 ms |
コード長 | 1,423 bytes |
コンパイル時間 | 92 ms |
コンパイル使用メモリ | 11,076 KB |
実行使用メモリ | 8,600 KB |
最終ジャッジ日時 | 2023-09-20 10:01:15 |
合計ジャッジ時間 | 1,097 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge12 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 18 ms
8,412 KB |
testcase_01 | AC | 18 ms
8,476 KB |
testcase_02 | AC | 18 ms
8,600 KB |
testcase_03 | AC | 17 ms
8,432 KB |
testcase_04 | AC | 19 ms
8,464 KB |
testcase_05 | AC | 19 ms
8,400 KB |
testcase_06 | AC | 19 ms
8,464 KB |
testcase_07 | AC | 19 ms
8,512 KB |
testcase_08 | AC | 19 ms
8,540 KB |
testcase_09 | AC | 19 ms
8,464 KB |
testcase_10 | AC | 18 ms
8,428 KB |
testcase_11 | AC | 17 ms
8,432 KB |
testcase_12 | AC | 20 ms
8,420 KB |
ソースコード
import collections def nCb_mod_p(n, b, mod): ''' nCb % mod を返す。 mod が素数のときのみ有効。 ''' if b > n - b: b = n - b num = 1 mod_counter = 0 for k in range(n, n-b, -1): if k % mod: num *= k else: mod_counter += 1 num *= (k // mod) num %= mod den = 1 for k in range(1, b+1): if k % mod: den *= k else: mod_counter -= 1 den *= (k // mod) den %= mod if mod_counter > 0: return 0 r = inv(den, mod) return (num * r) % mod def inv(a, m): ''' a * b = 1 (mod m) となる b を返す。 m は素数とする。 b = a**(m-2) (mod m) を使って計算。 ''' n = m - 2 b = bin(n)[2:][-1::-1] ret = 1 tmp = a if b[0] == '1': ret = a for bi in b[1:]: tmp *= tmp tmp %= m if bi == '1': ret *= tmp ret %= m return ret def solve(): S = input().rstrip() dic = collections.Counter(S) n = len(S) # 573 = 3 * 191 r3 = 1 r191 = 1 for b in dic.values(): r3 *= nCb_mod_p(n, b, 3) r3 %= 3 r191 *= nCb_mod_p(n, b, 191) r191 %= 191 n -= b for r in range(r191, 573, 191): if r % 3 == r3: return (r - 1) % 573 if __name__ == '__main__': print(solve())