結果

問題 No.171 スワップ文字列(Med)
ユーザー ayaoniayaoni
提出日時 2020-12-27 21:16:23
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,824 bytes
コンパイル時間 164 ms
コンパイル使用メモリ 82,420 KB
実行使用メモリ 55,580 KB
最終ジャッジ日時 2024-04-09 20:16:24
合計ジャッジ時間 1,435 ms
ジャッジサーバーID
(参考情報)
judge1 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 41 ms
54,512 KB
testcase_01 AC 40 ms
54,112 KB
testcase_02 WA -
testcase_03 AC 41 ms
54,480 KB
testcase_04 AC 41 ms
54,676 KB
testcase_05 AC 41 ms
55,408 KB
testcase_06 WA -
testcase_07 AC 42 ms
54,316 KB
testcase_08 AC 42 ms
54,868 KB
testcase_09 AC 42 ms
54,932 KB
testcase_10 AC 41 ms
54,520 KB
testcase_11 AC 43 ms
54,360 KB
testcase_12 AC 41 ms
54,176 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
from collections import Counter
sys.setrecursionlimit(10**7)
def I(): return int(sys.stdin.readline().rstrip())
def MI(): return map(int,sys.stdin.readline().rstrip().split())
def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))
def LI2(): return list(map(int,sys.stdin.readline().rstrip()))
def S(): return sys.stdin.readline().rstrip()
def LS(): return list(sys.stdin.readline().rstrip().split())
def LS2(): return list(sys.stdin.readline().rstrip())


def factorial(N,p):  # 0!~N!(mod p) のリスト
    res = [0]*(N+1)
    res[0] = 1
    for i in range(1,N+1):
        res[i] = (res[i-1]*i) % p
    return res


def ext_gcd(a, b):  # a*x+b*y == gcd(a,b) たる gcd(a,b),x,y
    if b > 0:
        d,x,y = ext_gcd(b,a % b)
        return d,y,x-(a//b)*y
    return a,1,0


# V = [(Xi,Yi),…]

def remainder(V):  # Z == Xi (mod Yi) たる Z, lcm(Yi)
    x,lcm = 0,1
    for X,Y in V:
        g,a,b = ext_gcd(lcm,Y)
        x,lcm = (Y*b*x+lcm*a*X)//g,lcm*(Y//g)
        x %= lcm
    return x,lcm


S = S()
N = len(S)
count = Counter(S)
X = list(count.values())

fac191 = factorial(190,191)


def f(x):  # x!が何回3で割り切れるか
    res = 0
    while x >= 3:
        res += x//3
        x //= 3
    return res


r = f(N)
ans3 = pow(2,N//3,3)
if N % 3 == 2:
    ans3 *= 2
    ans3 %= 3

for x in X:
    r -= f(x)
    ans3 *= pow(2,x//3,3)
    ans3 %= 3
    if x % 3 == 2:
        ans3 %= 2
        ans3 %= 3
if r:
    ans3 = 0

s = N//191
ans191 = 1
while N >= 191:
    ans191 *= -1
    ans191 %= 191
    N -= 191
ans191 *= fac191[N]
ans191 %= 191

for x in X:
    s -= x//191
    while x >= 191:
        ans191 *= -1
        x -= 191
    ans191 *= pow(fac191[x],189,191)
    ans191 %= 191
if s:
    ans191 = 0

ans,l = remainder([(ans3,3),(ans191,191)])
ans -= 1
ans %= 573
print(ans)
0