首页 > 代码库 > bzoj2004
bzoj2004
http://www.lydsy.com/JudgeOnline/problem.php?id=2004
矩阵快速幂。。。并没有想出来。。。
理解地不是很好 先挖个坑
#include<bits/stdc++.h> using namespace std; typedef long long ll; const int N = 210, mod = 30031; struct mat { int a[N][N]; } x, A; int n, k, p, cnt; int st[N]; mat operator * (mat A, mat B) { mat ret; memset(ret.a, 0, sizeof(ret.a)); for(int i = 1; i <= cnt; ++i) for(int j = 1; j <= cnt; ++j) for(int k = 1; k <= cnt; ++k) ret.a[i][j] = (ret.a[i][j] + A.a[i][k] * B.a[k][j]) % mod; return ret; } mat power(mat x, int t) { mat ret; memset(ret.a, 0, sizeof(ret.a)); for(int i = 1; i <= cnt; ++i) ret.a[i][i] = 1; for(; t; t >>= 1, x = x * x) if(t & 1) ret = ret * x; return ret; } void dfs(int d, int pos, int num) { if(d == k) { st[++cnt] = num; return; } for(int i = pos; i; --i) dfs(d + 1, i - 1, num | (1 << (i - 1))); } int main() { scanf("%d%d%d", &n, &k, &p); dfs(1, p - 1, 1 << (p - 1)); for(int i = 1; i <= cnt; ++i) for(int j = 1; j <= cnt; ++j) { int s = st[i] ^ (st[j] << 1) ^ (1 << p); if(s == (s & -s)) x.a[i][j] = 1; } x = power(x, n - k); A.a[1][1] = 1; A = x * A; printf("%d\n", A.a[1][1]); return 0; }
bzoj2004
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。