QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#589512 | #5139. DFS Order 2 | Soestx | ML | 0ms | 210848kb | C++17 | 3.0kb | 2024-09-25 18:13:58 | 2024-09-25 18:14:00 |
Judging History
answer
#include<bits/stdc++.h>
using namespace std;
#define int long long
#define fi first
#define se second
#define lowbit(x) (x&(-x))
typedef long long ll;
//typedef unsigned long long ull;
const int N = 1e6 + 10, M = 1e5,mod=998244353;
int n, m, k;
int res;
int dp[510][510][510],dap[510],tp[510][510];
int siz[510],son[510];
int h[510],e[1010],ne[1010],idx;
void add(int a,int b)
{
e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
int fac[510];
void ini()
{
idx=0;
fac[0]=1;
for(int i=1;i<=n;i++)
{
h[i]=-1;
fac[i]=fac[i-1]*i%mod;
}
}
ll qpow(ll a,ll b)
{
ll res=1;
while(b)
{
if(b&1) res=res*a%mod;
a=a*a%mod;
b>>=1;
}
return res;
}
int tap[510];
void dfs(int id,int f)
{
siz[id]=1;
dap[id]=1;
int cnt=0;
dp[id][0][0]=1;
for(int i=h[id];i!=-1;i=ne[i])
{
int j=e[i];
if(j==f) continue;
dfs(j,id);
for(int o=siz[id];o>=0;o--)
{
for(int p=cnt;p>=0;p--)
{
int to=o+siz[j];
dp[id][to][p+1]=(dp[id][to][p+1]+1ll*dp[id][o][p]*(p+1)%mod)%mod;
}
}
cnt++;
siz[id]+=siz[j];
dap[id]=1ll*dap[id]*dap[j]%mod;
}
son[id]=cnt;
tap[id]=dap[id];
dap[id]=1ll*dap[id]*fac[cnt]%mod;
}
int tmp[510][510],arr[510];
void DFS(int id,int f)
{
for(int i=h[id];i!=-1;i=ne[i])
{
int j=e[i];
if(j==f) continue;
//cout<<j<<"-----------------------"<<endl;
for(int a=0;a<=siz[id];a++) for(int b=0;b<=son[id];b++) tmp[a][b]=0,arr[a]=0;
for(int o=0;o<=siz[id];o++)
{
for(int p=0;p<=son[id];p++)
{
if(!dp[id][o][p]) continue;
int t=dp[id][o][p];;
if(o>=siz[j]&&p)
t=(1ll*dp[id][o][p]-p*tmp[o-siz[j]][p-1]%mod);
//cout<<o<<" "<<p<<" "<<dp[id][o][p]<<" "<<dp[id][o-siz[j]][p-1]<<" "<<t<<endl;
t=(t%mod+mod)%mod;
tmp[o][p]=(tmp[o][p]+t)%mod;
}
}
int ty=1ll*tap[id]*qpow(dap[j],mod-2)%mod;
for(int o=0;o<=siz[id];o++)
{
int sum=0;
for(int p=0;p<=son[id];p++)
{
//cout<<o<<" "<<p<<" "<<tmp[o][p]<<endl;
sum=(1ll*sum+tmp[o][p]*fac[son[id]-1-p]%mod)%mod;
}
arr[o]=1ll*sum*ty%mod;
}
for(int o=0;o<=n;o++)
{
for(int p=0;p<=siz[id];p++)
{
int to=o+p+1;
tp[j][to]=(tp[j][to]+1ll*tp[id][o]*arr[p]%mod)%mod;
}
}
}
for(int i=h[id];i!=-1;i=ne[i])
{
int j=e[i];
if(j==f) continue;
DFS(j,id);
}
}
void solve()
{
cin>>n;
ini();
for(int i=1;i<n;i++)
{
int a,b;
cin>>a>>b;
add(a,b);
add(b,a);
}
tp[1][1]=1;
dfs(1,0);
DFS(1,0);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
//cout<<i<<" "<<j<<" "<<tp[i][j]<<endl;
//cout<<tp[i][j]<<" "<<dap[i]<<endl;
cout<<1ll*tp[i][j]*dap[i]%mod<<" ";
//printf("%d ",1ll*tp[i][j]*dap[i]%mod);
}
cout<<endl;
}
}
signed main() {
ios::sync_with_stdio(0);
cin.tie(0);
cout.tie(0);
int T;
T = 1;
//cin>>T;
while (T--) solve();
return 0;
}
/*
1
2 0
10 15
*/
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 0ms
memory: 18036kb
input:
5 1 2 1 3 3 4 3 5
output:
4 0 0 0 0 0 2 0 0 2 0 2 2 0 0 0 0 1 2 1 0 0 1 2 1
result:
ok 25 numbers
Test #2:
score: 0
Accepted
time: 0ms
memory: 26164kb
input:
10 9 2 9 6 10 5 1 5 1 6 9 3 5 8 4 3 7 9
output:
24 0 0 0 0 0 0 0 0 0 0 0 0 4 2 2 8 2 2 4 0 0 0 4 4 4 4 4 4 0 0 0 0 0 4 4 4 4 4 4 0 12 0 0 0 0 0 12 0 0 0 12 0 0 12 0 0 0 0 0 0 0 0 4 2 2 8 2 2 4 0 0 6 6 0 0 0 0 6 6 0 0 12 0 0 12 0 0 0 0 0 0 6 6 0 0 0 0 6 6
result:
ok 100 numbers
Test #3:
score: 0
Accepted
time: 0ms
memory: 210848kb
input:
100 18 100 91 87 28 83 11 98 51 52 24 91 72 53 18 19 89 16 77 35 26 25 73 16 96 70 56 44 69 10 63 30 54 95 39 66 58 98 8 71 58 65 74 73 2 64 12 19 32 81 31 54 43 41 84 59 55 75 72 81 59 37 10 94 93 2 64 47 13 32 36 84 28 22 30 28 25 77 47 6 80 52 54 17 23 40 47 88 49 53 65 27 99 59 25 70 91 9 74 1 7...
output:
8388559 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 62914557 0 62914557 62914557 0 62914557 0 0 0 62914557 0 62914557...
result:
ok 10000 numbers
Test #4:
score: -100
Memory Limit Exceeded
input:
500 382 156 418 376 91 15 142 274 449 174 375 82 118 175 421 458 361 222 14 474 11 324 368 341 227 424 231 249 81 435 250 271 118 38 147 61 124 408 135 1 244 316 301 80 39 313 90 118 290 465 465 250 277 341 8 105 319 373 305 379 309 200 180 398 47 489 463 259 173 492 494 343 251 193 111 32 401 270 4...
output:
219078761 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ...