結果
問題 |
No.3178 free sort
|
ユーザー |
|
提出日時 | 2025-06-13 21:36:42 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,059 bytes |
コンパイル時間 | 583 ms |
コンパイル使用メモリ | 82,820 KB |
実行使用メモリ | 248,588 KB |
最終ジャッジ日時 | 2025-06-13 21:37:02 |
合計ジャッジ時間 | 17,090 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | WA * 40 |
ソースコード
def cmb(n, r, p):#nCr mod p if (r < 0) or (n < r): return 0 r = min(r, n - r) return fact[n] * factinv[r] * factinv[n-r] % p p = 998244353 binomial_max = 10 ** 6 # N は必要分だけ用意する fact = [1, 1] factinv = [1, 1] inv = [0, 1] for i in range(2, binomial_max + 1): fact.append((fact[-1] * i) % p) inv.append((-inv[p % i] * (p // i)) % p) factinv.append((factinv[-1] * inv[-1]) % p) N = list(input()) zero = N.count("0") cnt_l = [N.count(str(i)) for i in range(10)] pt = fact[len(N)] for i in range(10): pt *= factinv[cnt_l[i]] pt %= p pt2 = 0 for i in range(1,zero+1): rest_zero = zero-i for j in range(1,10): if cnt_l[j] > 0: tmp = fact[len(N)-1-i] tmp *= factinv[rest_zero] for k in range(1,10): if j == k: tmp *= factinv[cnt_l[k]-1] else: tmp *= factinv[cnt_l[k]] pt2 += tmp #print(j,tmp) pt -= pt2 pt %= p print(pt)