记录编号 |
606669 |
评测结果 |
AAAAAAAAAA |
题目名称 |
3946.信使 |
最终得分 |
100 |
用户昵称 |
dream |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
12.225 s |
提交时间 |
2025-10-01 17:10:06 |
内存使用 |
11.87 MiB |
显示代码纯文本
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=205;
int n,m,q;
ll mod;
ll ans[55][N][N],g[N][N],f[N][N][55],h[N][N][55];//h i j:i-i /ij |f i j:i-j /ij
int u,v,d;
int main(){
freopen("messenger.in","r",stdin);
freopen("messenger.out","w",stdout);
ios::sync_with_stdio(0);
cin>>n>>m>>mod;
for(int i=1;i<=m;i++){
int a,b;
cin>>a>>b;
ans[1][a][b]=g[a][b]=1;
f[a][b][1]=1;
}
for(int i=1;i<=n;i++){
if(g[i][i]){
for(int j=1;j<=n;j++){
if(i==j) continue;
h[i][j][1]=1;
}
}
}
for(int id=1;id<=50;id++){
for(int k=1;k<=n;k++){
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
ans[id][i][j]+=ans[id-1][i][k]*g[k][j]%mod;
ans[id][i][j]%=mod;
}
}
}
}
for(int idx=2;idx<=50;idx++){
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
ll s1=0,s2=0;
for(int k=1;k<idx;k++){
s1+=f[i][j][k]*ans[idx-k][j][j]%mod+h[i][j][k]*ans[idx-k][i][j]%mod;
s2+=f[i][j][k]*ans[idx-k][j][i]%mod+h[i][j][k]*ans[idx-k][i][i]%mod;
s1%=mod;
s2%=mod;
}
f[i][j][idx]=ans[idx][i][j]+mod-s1;
f[i][j][idx]%=mod;
h[i][j][idx]=ans[idx][i][i]+mod-s2;
h[i][j][idx]%=mod;
}
}
}
cin>>q;
while(q--){
cin>>u>>v>>d;
cout<<f[u][v][d]%mod<<"\n";
}
return 0;
}