BZOJ 1875 [SDOI2009]HH去散步 ——动态规划 矩阵乘法
发现t非常大,所以大概就是快速幂一类的问题了,
然后根据k^3logn算了算,发现k大约是边数的时候复杂度比较合适。
发现比较麻烦的就是前驱的记录,所以直接把边看做点,不能走反向边,但是可以走重边,然后t--
之后弄出状态转移矩阵递推即可。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 | #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; #define F(i,j,k) for (int i=j;i<=k;++i) #define D(i,j,k) for (int i=j;i>=k;--i) #define ll long long const ll md=45989; int n,m,t,A,B; int h[150],to[150],ne[150],fr[150],en=0,ans; void add( int a, int b) {to[en]=b;ne[en]=h[a];fr[en]=a;h[a]=en++;} struct Matrix{ ll x[125][125]; void init(){ memset (x,0, sizeof x);} Matrix operator * (Matrix y) { Matrix ret; ret.init(); F(i,0,en-1) F(k,0,en-1) { F(j,0,en-1) ret.x[i][k]+=x[i][j]*y.x[j][k]; ret.x[i][k]%=md; } return ret; } }a,b; int main() { memset (h,-1, sizeof h); scanf ( "%d%d%d%d%d" ,&n,&m,&t,&A,&B);t--; a.init(); b.init(); F(i,1,m) { int a,b; scanf ( "%d%d" ,&a,&b); add(a,b);add(b,a); } F(i,0,en-1) { int st=to[i]; for ( int j=h[st];j>=0;j=ne[j]) if ((i^1)!=j)b.x[i][j]++; } for ( int i=h[A];i>=0;i=ne[i]) a.x[0][i]=1; while (t) { if (t&1) a=a*b; b=b*b; t>>=1; } F(i,0,en-1) if (to[i]==B) (ans+=a.x[0][i])%=md; printf ( "%d\n" ,ans); } |
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
· 基于Microsoft.Extensions.AI核心库实现RAG应用
· 10年+ .NET Coder 心语 ── 封装的思维:从隐藏、稳定开始理解其本质意义
· 地球OL攻略 —— 某应届生求职总结
· 提示词工程——AI应用必不可少的技术
· Open-Sora 2.0 重磅开源!
· 周边上新:园子的第一款马克杯温暖上架