国内四大高端建站公司,网站建设内部下单流程图,网站建设兆金手指花总,北京 营销型网站P2151 [SDOI2009]HH去散步
题意#xff1a;
HH有个一成不变的习惯#xff0c;喜欢饭后百步走。所谓百步走#xff0c;就是散步#xff0c;就是在一定的时间 内#xff0c;走过一定的距离。 但是同时HH又是个喜欢变化的人#xff0c;所以他不会立刻沿着刚刚走来的路走回…P2151 [SDOI2009]HH去散步
题意
HH有个一成不变的习惯喜欢饭后百步走。所谓百步走就是散步就是在一定的时间 内走过一定的距离。 但是同时HH又是个喜欢变化的人所以他不会立刻沿着刚刚走来的路走回。 又因为HH是个喜欢变化的人所以他每天走过的路径都不完全一样他想知道他究竟有多 少种散步的方法。
现在给你学校的地图假设每条路的长度都是一样的都是1问长度为t从给定地 点A走到给定地点B共有多少条符合条件的路径
题解
这题跟裸的矩阵快速幂唯一区别是本题中不允许刚经过的边本题中给的是无向边顿时没有什么思路。 我们可以点边互换将边成是点然后再重新建边如果原图中有4个边分别编号1~4然后我们依次判断第i个边的尾和第j个边的头是否连接如果连接就从第i个边向第j个边连一条边 如何处理走回去的情况只要我们同一条边拆出来的两个点互不连边即可
代码
// Problem: P2151 [SDOI2009]HH去散步
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P2151
// Memory Limit: 125 MB
// Time Limit: 1000 ms
// Data:2021-08-13 16:47:00
// By Jozky
#include bits/stdc.h
#include unordered_map
#define debug(a, b) printf(%s %d\n, a, b);
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pairint, int PII;
clock_t startTime, endTime;
//Fe~Jozky
const ll INF_ll 1e18;
const int INF_int 0x3f3f3f3f;
template typename T inline void read(T x)
{T f 1;x 0;char ch getchar();while (0 isdigit(ch)) {if (ch -)f -1;ch getchar();}while (0 ! isdigit(ch))x (x 1) (x 3) ch - 0, ch getchar();x* f;
}
template typename T inline void write(T x)
{if (x 0) {x ~(x - 1);putchar(-);}if (x 9)write(x / 10);putchar(x % 10 0);
}const int maxn 151, mod 45989;
int n, m, s, t, x[maxn], y[maxn], cnt;
ll T;
struct Ma
{ll m[maxn][maxn];
};
Ma mul(Ma a, Ma b, int n)
{Ma temp;memset(temp, 0, sizeof(temp));for (int i 1; i n; i) {for (int j 1; j n; j) {for (int k 1; k n; k) {temp.m[i][j] (temp.m[i][j] ((a.m[i][k] % mod) * (b.m[k][j]) % mod) % mod) % mod;}}}return temp;
}
Ma poww(Ma res, ll N, int n)
{Ma ans;for (int i 1; i n; i) {ans.m[i][i] 1;}while (N) {if (N 1)ans mul(ans, res, n);res mul(res, res, n);N 1;}return ans;
}
Ma A;
void add(int u, int v)
{x[cnt] u;y[cnt] v;
}
int main()
{//rd_test();scanf(%d%d%lld%d%d, n, m, T, s, t);s;t;x[cnt] 0;y[cnt] s;for (int i 1, u, v; i m; i) {scanf(%d%d, u, v);u, v;x[cnt] u, y[cnt] v;x[cnt] v, y[cnt] u;}for (int i 1; i cnt; i)for (int j 1; j cnt; j)if (i ! j i ! (j ^ 1)) {if (y[i] x[j])A.m[i][j] 1;}Ma Ans poww(A, T, cnt);int ans 0;for (int i 1; i cnt; i)if (y[i] t) {ans (ans Ans.m[1][i]) % mod;}cout ans endl;return 0;//Time_test();
}
