結果

問題 No.1655 123 Swaps
ユーザー kozykozy
提出日時 2021-08-21 01:26:56
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,038 ms / 2,000 ms
コード長 2,972 bytes
コンパイル時間 158 ms
コンパイル使用メモリ 81,948 KB
実行使用メモリ 208,788 KB
最終ジャッジ日時 2024-10-14 10:15:33
合計ジャッジ時間 25,110 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 155 ms
175,304 KB
testcase_01 AC 156 ms
174,940 KB
testcase_02 AC 198 ms
188,196 KB
testcase_03 AC 159 ms
176,228 KB
testcase_04 AC 158 ms
174,836 KB
testcase_05 AC 156 ms
174,848 KB
testcase_06 AC 156 ms
175,916 KB
testcase_07 AC 156 ms
176,428 KB
testcase_08 AC 156 ms
175,736 KB
testcase_09 AC 159 ms
177,448 KB
testcase_10 AC 184 ms
184,788 KB
testcase_11 AC 172 ms
179,848 KB
testcase_12 AC 185 ms
186,372 KB
testcase_13 AC 160 ms
174,524 KB
testcase_14 AC 159 ms
175,600 KB
testcase_15 AC 1,033 ms
207,968 KB
testcase_16 AC 1,031 ms
207,988 KB
testcase_17 AC 1,033 ms
208,228 KB
testcase_18 AC 1,028 ms
208,144 KB
testcase_19 AC 1,035 ms
208,080 KB
testcase_20 AC 1,036 ms
208,100 KB
testcase_21 AC 1,038 ms
208,008 KB
testcase_22 AC 1,038 ms
208,360 KB
testcase_23 AC 1,024 ms
208,788 KB
testcase_24 AC 1,027 ms
208,660 KB
testcase_25 AC 1,036 ms
208,416 KB
testcase_26 AC 1,029 ms
207,200 KB
testcase_27 AC 627 ms
199,628 KB
testcase_28 AC 607 ms
197,928 KB
testcase_29 AC 1,022 ms
205,100 KB
testcase_30 AC 1,011 ms
204,204 KB
testcase_31 AC 621 ms
197,768 KB
testcase_32 AC 161 ms
174,872 KB
testcase_33 AC 165 ms
174,932 KB
testcase_34 AC 163 ms
175,152 KB
testcase_35 AC 164 ms
176,644 KB
testcase_36 AC 598 ms
197,780 KB
testcase_37 AC 604 ms
197,904 KB
testcase_38 AC 607 ms
198,004 KB
testcase_39 AC 607 ms
197,684 KB
testcase_40 AC 1,027 ms
207,876 KB
testcase_41 AC 164 ms
175,628 KB
testcase_42 AC 166 ms
175,544 KB
testcase_43 AC 165 ms
175,296 KB
testcase_44 AC 162 ms
175,580 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

mod=924844033
root=pow(5,441,mod)
rootinv=pow(root,mod-2,mod)
bbeki=[0]*22
bbeki2=[0]*22
beki_dic=dict()
this=1
for i in range(22):
  beki_dic[this]=i
  this*=2
def beki_mae():
  x=root
  beki=1
  for i in range(21,-1,-1):
    bbeki[i]=x
    x=x**2
    x%=mod
  x=rootinv
  beki=1
  for i in range(21,-1,-1):
    bbeki2[i]=x
    x=x**2
    x%=mod
beki_mae()
def rev(i,k):
  ans=0
  s=bin(i)[2:]
  for i in range(len(s)):
    if s[len(s)-1-i]=="1":
      ans+=2**(k-1-i)
  return ans
def NTT_len(a):#a<=2^Nを満たす最小のNを返す
  k=1
  b=0
  for i in range(100):
    if k>=a:
      break
    else:
      k*=2
      b+=1
  return b
def NTT_change0(A,k):#長さ2^kにする
  N=len(A)
  A=A+[0]*(2**k-N)
  NTT_change(A)
  return A
def NTT_change(A):#Aをstart<=k<finまでNTT変換,f(x^i)(i=0~N-1)を求める
  N=len(A)
  le=N
  while le>1:
    x=bbeki[beki_dic[le]]
    for st in range(N//le):
      this=1
      for i in range(le//2):
        l=A[st*le+i]
        r=A[st*le+i+le//2]
        A[st*le+i]=(l+r)%mod
        A[st*le+i+le//2]=((l-r)*this)%mod
        this*=x
        this%=mod
    le//=2

def NTT_invchange0(A,k):#長さ2^kにする
  N=len(A)
  A=A+[0]*(2**k-N)
  NTT_invchange(A)
  return A
def NTT_invchange(A):#Aをstart<=k<finまでNTT変換,f(x^i)(i=0~N-1)を求める
  N=len(A)
  le=2
  while le<=N:
    x=bbeki2[beki_dic[le]]
    for st in range(N//le):
      this=1
      for i in range(le//2):
        l=A[st*le+i]
        r=A[st*le+i+le//2]*this
        r%=mod
        A[st*le+i]=(l+r)%mod
        A[st*le+i+(le//2)]=(l-r)%mod
        this*=x
        this%=mod
    le*=2
  invN=pow(N,mod-2,mod)
  for i in range(N):
    A[i]*=invN
    A[i]%=mod
def NTT_time(A,B):#A,Bの畳み込み
  n=len(A)
  m=len(B)
  k=NTT_len(n+m-1)
  A=NTT_change0(A,k)
  B=NTT_change0(B,k)
  c=list()
  for i in range(len(A)):
    c.append(A[i]*B[i]%mod)
  c=NTT_invchange0(c,k)
  return c[:n+m-1]

def NTTinv(f):
  """
  1/f=gをreturn
  """
  le=len(f)
  this_roop=NTT_len(len(f))
  c=pow(int(f[0]),mod-2,mod)
  g=[c]
  a=1
  for i in range(this_roop):
      a*=2
      S=NTT_time(g,f)[:a]
      S=[-i%mod for i in S]
      S[0]+=2
      g=NTT_time(g,S)[:a]
  return g[:le]
g1 = [1, 1]
g2 = [1, 1]
inverse = [0, 1]
N=5*(10**5)
for i in range( 2, N + 1 ):
    g1.append( ( g1[-1] * i ) % mod )
    inverse.append( ( -inverse[mod % i] * (mod//i) ) % mod )
    g2.append( (g2[-1] * inverse[-1]) % mod )
a,b,c=map(int,input().split())
A=a
B=b
C=c
N=a+b+c
if N%2==1:
  print(0)
  exit()
L=[[0]*(A+1) for i in range(3)]
R=[[0]*(B+1) for i in range(3)]
for x in range(A+1):
  L[x%3][x]=g2[x]*g2[A-x]%mod
for y in range(B+1):
  R[y%3][y]=g2[y]*g2[B-y]%mod
C=[0]*(A+B+1)
for x in range(3):
  S=NTT_time(L[x],R[((B+2*x-A)*2)%3])
  for i in range(len(S)):
    C[i]+=S[i]
    C[i]%=mod
for i in range(A+B+1):
  C[i]*=(g1[(N//2)]**2)*g2[(N//2)-i]*g2[(N//2)-A-B+i]%mod
  C[i]%=mod
ans=0
for i in range((N//2)-c,(N//2)+1):
  if 0<=i<=A+B:
    ans+=C[i]
    ans%=mod
print(ans)
0