結果

問題 No.3486 Draw a Rainbow
コンテスト
ユーザー mkawa2
提出日時 2026-03-27 23:13:27
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
AC  
実行時間 692 ms / 4,000 ms
コード長 1,722 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 176 ms
コンパイル使用メモリ 85,120 KB
実行使用メモリ 233,360 KB
最終ジャッジ日時 2026-03-27 23:13:35
合計ジャッジ時間 6,974 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge1_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import sys

# sys.setrecursionlimit(200005)
# sys.set_int_max_str_digits(200005)
int1 = lambda x: int(x)-1
pDB = lambda *x: print(*x, end="\n", file=sys.stderr)
p2D = lambda x: print(*x, sep="\n", end="\n\n", file=sys.stderr)
def II(): return int(sys.stdin.readline())
def LI(): return list(map(int, sys.stdin.readline().split()))
def LLI(rows_number): return [LI() for _ in range(rows_number)]
def LI1(): return list(map(int1, sys.stdin.readline().split()))
def LLI1(rows_number): return [LI1() for _ in range(rows_number)]
def SI(): return sys.stdin.readline().rstrip()

dij = [(0, 1), (-1, 0), (0, -1), (1, 0)]
# dij = [(0, 1), (-1, 0), (0, -1), (1, 0), (1, 1), (1, -1), (-1, 1), (-1, -1)]
# inf = -1-(-1 << 31)
inf = -1-(-1 << 62)

# md = 10**9+7
md = 998244353

def zeta(aa):
    n = len(aa).bit_length()-1
    res = aa[:]
    for i in range(n):
        for s in range(1 << n):
            if s >> i & 1: res[s] += res[s ^ 1 << i]
    return res

def mobius(aa):
    n = len(aa).bit_length()-1
    res = aa[:]
    for i in range(n):
        for s in range(1 << n):
            if s >> i & 1: res[s] -= res[s ^ 1 << i]
    return res

def dodp(cnt):
    dp=[0]*(1<<l)
    dp[0]=1
    for s in range(1,1<<l):
        b=s.bit_length()-1
        dp[s]=dp[s^1<<b]*(pw[cnt[b]]-1)%md
    dp[0]=0
    return dp

pw=[1]
for _ in range(200005):pw.append(pw[-1]*2%md)

n,m,l=LI()
a2b=[[0]*l for _ in range(m)]
for _ in range(n):
    a,b=LI1()
    b-=m
    a2b[a][b]+=1
dps=[]
for a in range(m):
    dps.append(dodp(a2b[a]))
# print(dps)

ans=[1]*(1<<l)
for a in range(m):
    dps[a]=zeta(dps[a])
    for s in range(1<<l):
        ans[s]*=dps[a][s]
        ans[s]%=md
# print(dps)
# print(ans)

ans=mobius(ans)
print(ans[-1]%md)
0