G = []
# 2, 3, ..., 50 からなる有向サイクル
G += [(n, n+1) for n in range(2, 50)] + [(50,2)]
# 52, 53, ..., 100 からなる完全グラフ
G += [(a, b) for a in range(52, 101) for b in range(52, 101)]
# 余分な弦、1 → サイクル、サイクル → v、v → 完全
G += [(10,12), (1,20), (30,51), (51,52)]
N = 100
M = len(G)
print(f'{N=}, {M=}')

N = 100
dp = [0] * (N+1)
dp[1] = 1
max_T = 10000
ans = [0] * (max_T + 1)
for t in range(max_T + 1):
    ans[t] = sum(dp)
    newdp = [0] * (N+1)
    for a, b in G:
        newdp[b] |= dp[a]
    dp = newdp

# 各時刻ごとの答。
# 解が収束して周期にはいるまでに 2300 秒ほどかかることが確かめられる。
# そこまでの dp 遷移の回数は、M 倍して 7_000_000 回ほど。
for t, x in enumerate(ans):
    print(t, x)