题解:
我们相当于寻找走到每个点的长度的期望,对于每个点对总期望的贡献就是能走到这个点的次数乘以走到这个点的总长度。
#include <bits/stdc++.h>
#define int long long
using namespace std;
const int mod=998244353;
const int N=1e5+10;
int q_pow(int a,int b)
{
int res=1;
while(b){
if(b&1) res=res*a%mod;
a=a*a%mod;
b>>=1;
}
return res%mod;
}
vector<int> g[N];
int cnt[N],f[N];
void dfs(int u)
{
if(cnt[u]) return ;
cnt[u]=1;
for(auto it:g[u]){
dfs(it);
cnt[u]=(cnt[u]+cnt[it])%mod;
f[u]=(f[u]+f[it]+cnt[it])%mod;
}
}
signed main()
{
int n,m; scanf("%lld%lld",&n,&m);
for(int i=1;i<=m;i++){
int a,b; scanf("%lld%lld",&a,&b);
g[a].push_back(b);
}
// cout<<"-----"<<endl;
for(int i=1;i<=n;i++) if(!cnt[i]) dfs(i);
int sum1=0,sum2=0;
for(int i=1;i<=n;i++) sum1=(sum1+f[i])%mod,sum2=(sum2+cnt[i])%mod;
printf("%lld\n",sum1*q_pow(sum2,mod-2)%mod);
}