1875. [SDOI2009]HH去散步【矩阵乘法】
Description
Input
Output
一行,表示答案。
Sample Input
0 1
0 2
0 3
2 1
3 2
Sample Output
4
第一道矩阵乘法……orz调了半天发现自己犯了一个傻逼错误
多余的不解释了,学习笔记里有思路。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#define MOD (45989)
using namespace std;
struct Mar
{
int a[][];
} unit,G;
struct node1
{
int to,next;
} edge[];
int n,head[],num_edge; void add(int u,int v)
{
edge[num_edge].to=v;
edge[num_edge].next=head[u];
head[u]=num_edge++;
} Mar Mul(Mar x,Mar y)
{
Mar c;
memset(c.a,,sizeof(c.a));//初始化,不然会GG
for (int i=; i<=num_edge-; ++i)
for (int j=; j<=num_edge-; ++j)
for (int k=; k<=num_edge-; ++k)
c.a[i][j]=(c.a[i][j]+x.a[i][k]*y.a[k][j])%MOD;
return c;
} Mar Mar_pow(Mar a,int p)
{
Mar ans=unit;
while (p!=)
{
if (p&)
ans=Mul(ans,a);
a=Mul(a,a);
p>>=;
}
return ans;
} int main()
{
int m,t,st,ed,u,v;
scanf("%d%d%d%d%d",&n,&m,&t,&st,&ed);
++st;++ed;
for (int i=; i<=n; ++i) head[i]=-;
for (int i=; i<=m; ++i)
{
scanf("%d%d",&u,&v);
++u;++v;
add(u,v);add(v,u);
}
for (int i=; i<=num_edge-; ++i)
unit.a[i][i]=;//给单位矩阵赋值
for (int i=; i<=num_edge-; ++i)
for (int j=head[edge[i].to]; j!=-; j=edge[j].next)
if ((i!=(j^)))
G.a[i][j]++; G=Mar_pow(G,t-);
int Final=;
for (int i=head[st]; i!=-; i=edge[i].next)
for (int j=head[ed]; j!=-; j=edge[j].next)
Final=(Final+G.a[i][j^])%MOD;//异或就是将二进制最后一位取反,即入边变出边,出边变入边
printf("%d",Final);
}
1875. [SDOI2009]HH去散步【矩阵乘法】的更多相关文章
- bzoj 1875: [SDOI2009]HH去散步 -- 矩阵乘法
1875: [SDOI2009]HH去散步 Time Limit: 20 Sec Memory Limit: 64 MB Description HH有个一成不变的习惯,喜欢饭后百步走.所谓百步走, ...
- 【bzoj1875】[SDOI2009]HH去散步 矩阵乘法
题目描述 一张N个点M条边的无向图,从A走到B,要求:每一次不能立刻沿着上一次的边的反方向返回.求方案数. 输入 第一行:五个整数N,M,t,A,B. N表示学校里的路口的个数 M表示学校里的路的条数 ...
- [SDOI2009] HH去散步 (矩阵乘法)
link $solution:$ 将边化为点后重新建矩阵,跑$T-1$幂即可(因为跑的是新边). 最后直接找与$x,y$所相连的边即可. #include<iostream> #inclu ...
- 【BZOJ】1875: [SDOI2009]HH去散步 矩阵快速幂
[题意]给定n个点m边的无向图,求A到B恰好经过t条边的路径数,路径须满足每条边都和前一条边不同.n<=20,m<=60,t<=2^30. [算法]矩阵快速幂 [题解]将图的邻接矩阵 ...
- BZOJ 1875: [SDOI2009]HH去散步( dp + 矩阵快速幂 )
把双向边拆成2条单向边, 用边来转移...然后矩阵乘法+快速幂优化 ------------------------------------------------------------------ ...
- 洛谷P2151 [SDOI2009] HH去散步 [矩阵加速]
题目传送门 HH去散步 题目描述 HH有个一成不变的习惯,喜欢饭后百步走.所谓百步走,就是散步,就是在一定的时间 内,走过一定的距离. 但是同时HH又是个喜欢变化的人,所以他不会立刻沿着刚刚走来的路走 ...
- BZOJ 1875: [SDOI2009]HH去散步(矩阵乘法)
首先,题意就把我们引向了矩阵乘法,注意边长m<=60,那么就按边建图,变成一个120个点的图,然后乱搞就行了。 PS:WA了N久改了3次终于A了QAQ CODE: #include<cst ...
- bzoj 1875 [SDOI2009]HH去散步(矩乘)
Description HH有个一成不变的习惯,喜欢饭后百步走.所谓百步走,就是散步,就是在一定的时间 内,走过一定的距离. 但是同时HH又是个喜欢变化的人,所以他不会立刻沿着刚刚走来的路走回. 又因 ...
- bzoj1875 [SDOI2009]HH去散步 矩阵快速幂
题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=1875 题解 如果没有这个"不能立刻沿着刚刚走来的路走回",那么这个题就是一 ...
随机推荐
- MarkDownPad 专业汉化破解
解压Pa_ttrar 运行Pa_ttrar.exe 点击下边第一个按钮“patch”——>弹出窗选择“YES” 选择“YES”后会选择一个文件,找到“C:\Users\用户名\AppD ...
- vue2.0 element-ui中input的@keyup.native.enter='onQuery'回车查询刷新整个表单的解决办法
项目中用的element-ui是v1.4.3版本 实现的功能是在input中输入查询的名称,按下键盘回车键,可以查询表格中数据 问题是,我输入名称,按下回车,会整个表单刷新,搜索条件也被清空:代码如下 ...
- vue-cli脚手架的安装
https://github.com/vuejs/vue-cli 官网 使用官方推荐的webpack 条件:node在4.以上,npm在3以上,查看版本号打开cmd输入,node -v npm ...
- drupal7 formAPI给元素加css样式
比如,我定义了一个表单元素,名字叫做包库开始日期, $form['starttime_baoyue']=array( '#type'=>'textfield', '#title'=>t(' ...
- Json反序列化为动态类型(dynamic)
方法(依赖于Newtonsoft.Json): /// <summary> /// 反序列化json字符串 /// </summary> /// <typeparam n ...
- Eclipse 导入 Android studio Exception Ljava/lang/UnsatisfiedLinkEror
android studio compile fileTree(dir: 'libs', include: ['*.jar']) 没有加载so文件 main 下加入 jniLibs---so文件即可 ...
- Visualizing CNN Layer in Keras
CNN 权重可视化 How convolutional neural networks see the world An exploration of convnet filters with Ker ...
- Reducing and Profiling GPU Memory Usage in Keras with TensorFlow Backend
keras 自适应分配显存 & 清理不用的变量释放 GPU 显存 Intro Are you running out of GPU memory when using keras or ten ...
- RESULT_OK,RESULT_CANCELED,RESULT_FIRST_USER
RESULT_OK是执行结果,有RESULT_OK,RESULT_CANCELED,RESULT_FIRST_USER 在调用系统app时返回时RESULT_CANCELED如字面意思代表取消,RES ...
- Path2.0中绚丽的的旋转菜单
我们看一下实现的效果图: 在上图中,我将菜单弹出的效果设置成直线型,最终的弹出或汇总点在下面的红色按钮中. 它的实现原理是设置动画的同时并利用动画中的插入器(interpolator)来实现弹力.主要 ...