結果

問題 No.2504 NOT Path Painting
ユーザー chineristACchineristAC
提出日時 2023-10-13 22:39:28
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 782 ms / 2,000 ms
コード長 2,137 bytes
コンパイル時間 340 ms
コンパイル使用メモリ 82,452 KB
実行使用メモリ 107,156 KB
最終ジャッジ日時 2024-09-15 18:24:47
合計ジャッジ時間 14,359 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 51 ms
56,448 KB
testcase_01 AC 583 ms
79,304 KB
testcase_02 AC 581 ms
79,484 KB
testcase_03 AC 589 ms
79,264 KB
testcase_04 AC 581 ms
79,076 KB
testcase_05 AC 590 ms
79,056 KB
testcase_06 AC 574 ms
79,096 KB
testcase_07 AC 575 ms
78,592 KB
testcase_08 AC 574 ms
78,984 KB
testcase_09 AC 566 ms
78,632 KB
testcase_10 AC 580 ms
78,976 KB
testcase_11 AC 585 ms
78,976 KB
testcase_12 AC 618 ms
79,716 KB
testcase_13 AC 580 ms
78,540 KB
testcase_14 AC 603 ms
79,020 KB
testcase_15 AC 749 ms
101,932 KB
testcase_16 AC 782 ms
106,816 KB
testcase_17 AC 744 ms
105,116 KB
testcase_18 AC 694 ms
106,908 KB
testcase_19 AC 746 ms
104,692 KB
testcase_20 AC 685 ms
107,156 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
from itertools import permutations
import heapq
from collections import deque
import random

input = lambda :sys.stdin.readline().rstrip()
mi = lambda :map(int,input().split())
li = lambda :list(mi())

mod = 998244353

def solve_brute(n,edge):

    
    res = 1
    

    

    deq = deque([0])
    topo = []
    parent = [-1] * n
    while deq:
        v = deq.popleft()
        topo.append(v)
        for nv in edge[v]:
            if nv == parent[v]:
                continue
            parent[nv] = v
            deq.append(nv)
    
    sz = [1] * n
    for v in topo[::-1]:
        for nv in edge[v]:
            if nv == parent[v]:
                continue
            sz[v] += sz[nv]
        
        ok = n * (n+1)//2 - (n-sz[v]) * (n-sz[v]+1) //2
        for nv in edge[v]:
            if nv == parent[v]:
                continue
            ok -= sz[nv] * (sz[nv]+1) // 2 
        
        
        
        res += ok * pow(n*(n+1)//2-ok,mod-2,mod)
        res %= mod
    
    for v in range(1,n):
        s = sz[v]
        t = N - sz[v]

        ok = s * t

        res -= ok * pow(n*(n+1)//2-ok,mod-2,mod)
        res %= mod


    
    return res % mod


def brute(N):
    def inv(r):
        return pow(r,mod-2,mod)
    
    dp = [[0]*N for i in range(N)]
    for i in range(N):
        k = i * (i+1) //2 + (N-i-1) * (N-i)//2
        dp[i][i] = N*(N+1)//2 * pow(k,mod-2,mod) % mod
    
    
    
    for l in range(N)[::-1]:
        for r in range(l+1,N):
            a,b = 0,0
            for x in range(N):
                for y in range(x,N):
                    nl,nr = max(l,x),min(r,y)
                    if (nl,nr)!=(l,r):
                        b += dp[nl][nr] * inv(N*(N+1)//2)
                        b %= mod
                    else:
                        a += 1
            
            dp[l][r] = (b+1) * inv(1-(a * inv(N*(N+1)//2))) % mod
        
    return dp[0][N-1]

for _ in range(int(input())):
    N = int(input())
    edge = [[] for v in range(N)]
    for _ in range(N-1):
        u,v = mi()
        edge[u-1].append(v-1)
        edge[v-1].append(u-1)

    print(solve_brute(N,edge))
0