hihoCoder#1743:K-偏差排列(矩阵快速幂+状压dp)
题意
如果一个 \(1\to N\) 的排列 \(P=[P_1, P_2, ... P_N]\) 中的任意元素 \(P_i\) 都满足 \(|P_i-i| ≤ K\) ,我们就称 \(P\) 是 \(K\)-偏差排列。
给定 \(N\) 和 \(K\) ,请你计算一共有少个不同的排列是 \(K\)-偏差排列。
例如对于 \(N=3\) ,有 \(3\) 个 \(1\)-偏差排列:\([1, 2, 3], [1, 3, 2], [2, 1, 3]\)。
由于答案可能非常大,你只需要输出答案模 \(1000000007\) 的余数。
对于 \(70\%\) 的数据,\(1 ≤ N ≤ 1000\)
对于 \(100\%\) 的数据,\(1 ≤ N ≤ 1000000000, 1 ≤ K ≤ 3\)
题解
一道好题~
这是它的最初版本 #1732 : 1-偏差排列 .
那个找规律就是 斐波那契数列 了, dp 的话也是一样的结果 .
对于这个题我们可以沿用那题思路, 考虑一个位置 \(i\) 能放哪些数, 根据定义能放 \([i-k, i+k]\) 中共 \(2k+1\) 个数.
考虑状压到 \(i\) 这个点, 这些数中的哪些被放了, 每次转移的时候考虑放入一个数, 这个数之前不能出现, 这样就是合法转移了.
最后到 \(n\) 的时候, 不能放比 \(n\) 大的数, 且小于等于 \(n\) 的数都要放进去, 只会有那个位置存在正确答案, 这个状态 \(sta=2 ^ {k + 1} - 1\) (也就是意味着 \([n - k, n]\)都得选) .
当 \(n < k\) 的时候要特判掉一些诡异的特殊情况 .
然后这样直接写就有 \(70pts\) 了.
有一些不合法状态不能转移, 也就是要放的数不存在于 \([1, n]\) 之间.
这样的话, 就是矩阵快速幂套路优化了, 考虑对这个转移系数建立矩阵, 然后它的 \(n\) 次幂中的 \((sta,sta)\) 这个位置就会存在最后的答案咯...
代码
\(70pts:\)
#include <bits/stdc++.h>
#define For(i, l, r) for(register int i = (l), i##end = (int)(r); i <= i##end; ++i)
using namespace std;
typedef long long ll;
const ll Mod = 1e9 + 7;
int n, k, all;
ll ans = 0, dp[2][1500] = {0};
int main () {
cin >> n >> k;
if (n <= 2) return printf ("%d\n", n), 0;
all = (1 << (2 * k + 1)) - 1;
dp[0][0] = 1;
int cur = 0;
For (i, 1, n) {
For (j, 0, all) if(dp[cur][j]) {
int sta = (j >> 1);
For (s, 0, 2 * k) if (!(sta & (1 << s))) {
int tmp = i - k + s;
if (tmp < 1 || tmp > n) continue ;
(dp[cur ^ 1][sta | (1 << s)] += dp[cur][j]) %= Mod;
}
dp[cur][j] = 0;
}
cur ^= 1;
}
int Sta = (1 << (k + 1)) - 1;
printf ("%lld\n", dp[cur][Sta]);
return 0;
}
\(100pts:\)
#include <bits/stdc++.h>
#define For(i, l, r) for(register int i = (l), i##end = (int)(r); i <= i##end; ++i)
#define Set(a, v) memset(a, v, sizeof(a))
using namespace std;
void File() {
#ifdef zjp_shadow
freopen ("P1743.in", "r", stdin);
freopen ("P1743.out", "w", stdout);
#endif
}
const int Mod = 1e9 + 7, Maxn = 130;
int n, k, all;
struct Matrix {
int a[Maxn][Maxn]; Matrix() { Set(a, 0); }
void Unit() { For (i, 0, all) a[i][i] = 1; }
};
inline Matrix operator * (Matrix a, Matrix b) {
Matrix res;
For (i, 0, all) For (k, 0, all) if (a.a[i][k])
For (j, 0, all) (res.a[i][j] += 1ll * a.a[i][k] * b.a[k][j] % Mod) %= Mod;
return res;
}
inline Matrix fpm(Matrix x, int power) {
Matrix res; res.Unit();
for (; power; power >>= 1, x = x * x)
if (power & 1) res = res * x;
return res;
}
Matrix Bas, Ans;
int ans = 0;
int main () {
File();
cin >> n >> k;
if (n <= 2) return printf ("%d\n", n), 0;
all = (1 << (2 * k + 1)) - 1;
For (i, 0, all) {
int j = (i >> 1);
For (s, 0, 2 * k) if (!(j & (1 << s)))
++ Bas.a[i][j | (1 << s)];
}
Ans = fpm(Bas, n);
int Sta = (1 << (k + 1)) - 1;
printf ("%d\n", Ans.a[Sta][Sta]);
return 0;
}
hihoCoder#1743:K-偏差排列(矩阵快速幂+状压dp)的更多相关文章
- Codeforces Round #554 (Div. 2) F2. Neko Rules the Catniverse (Large Version) (矩阵快速幂 状压DP)
题意 有nnn个点,每个点只能走到编号在[1,min(n+m,1)][1,min(n+m,1)][1,min(n+m,1)]范围内的点.求路径长度恰好为kkk的简单路径(一个点最多走一次)数. 1≤n ...
- bzoj2004 [Hnoi2010]Bus 公交线路 矩阵快速幂+状压DP
题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=2004 题解 如果 \(N\) 没有那么大,考虑把每一位分配给每一辆车. 假设已经分配到了第 \ ...
- POJ 3744 【矩阵快速幂优化 概率DP】
搞懂了什么是矩阵快速幂优化.... 这道题的重点不是DP. /* 题意: 小明要走某条路,按照个人兴致,向前走一步的概率是p,向前跳两步的概率是1-p,但是地上有地雷,给了地雷的x坐标,(一维),求小 ...
- 【BZOJ5010】【FJOI2017】矩阵填数 [状压DP]
矩阵填数 Time Limit: 10 Sec Memory Limit: 128 MB[Submit][Status][Discuss] Description 给定一个 h*w 的矩阵,矩阵的行 ...
- poj 3744 Scout YYF I (矩阵快速幂 优化 概率dp)
题目链接 分析&&题意来自 : http://www.cnblogs.com/kuangbin/archive/2012/10/02/2710586.html 题意: 在一条不满地雷的 ...
- hdu 5564 Clarke and digits 矩阵快速幂优化数位dp
Clarke and digits Time Limit: 5000/3000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others ...
- POJ 3744 Scout YYF I(矩阵快速幂优化+概率dp)
http://poj.org/problem?id=3744 题意: 现在有个屌丝要穿越一个雷区,雷分布在一条直线上,但是分布的范围很大,现在这个屌丝从1出发,p的概率往前走1步,1-p的概率往前走2 ...
- 食物(矩阵快速幂)(DP)
这个题..我们可以想到用递推写!!qwq(好吧,其实我的DP水平不高啊qwq) 就是我们以两个为单位(一共九种组合情况),然后往后面推下一位的情况. 通过手动模拟,我们可以找到它们之间的递推关系(详见 ...
- [BZOJ5010][FJOI2017]矩阵填数(状压DP)
5010: [Fjoi2017]矩阵填数 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 90 Solved: 45[Submit][Status][ ...
随机推荐
- Luogu4156 WC2016 论战捆竹竿 KMP、同余类最短路、背包、单调队列
传送门 豪华升级版同余类最短路-- 官方题解 主要写几个小trick: \(1.O(nm)\)实现同余类最短路: 设某一条边长度为\(x\),那么我们选择一个点,在同余类上不断跳\(x\),可以形成一 ...
- 浅谈CDQ分治与偏序问题
初识CDQ分治 CDQ分治是一个好东西,一直听着dalao们说所以就去学了下. CDQ分治是我们处理各类问题的重要武器.它的优势在于可以顶替复杂的高级数据结构,而且常数比较小:缺点在于必须离线操作. ...
- 使用ajax方法实现form表单的提交
作者:13 GitHub:https://github.com/ZHENFENG13 版权声明:本文为原创文章,未经允许不得转载. 写在前面的话 在使用form表单的时候,一旦点击提交触发submit ...
- Facebook React 和 Web Components(Polymer)对比优势和劣势
目录结构 译者前言 Native vs. Compiled 原生语言对决预编译语言 Internal vs. External DSLs 内部与外部 DSLs 的对决 Types of DSLs - ...
- Iptables防火墙规则使用梳理
iptables是组成Linux平台下的包过滤防火墙,与大多数的Linux软件一样,这个包过滤防火墙是免费的,它可以代替昂贵的商业防火墙解决方案,完成封包过滤.封包重定向和网络地址转换(NAT)等功能 ...
- nginx的access.log文件详解
事实证明,日志文件真的是很重要很重要的.能够帮助我们快速的定位问题,并且知道用户访问的状态,浏览器,Ip,接口地址等,简直可怕.. 一.nginx的access.log(1)对博主而言,日志文件存放在 ...
- 这是一个数学题牛客训练赛E
题目描述 https://www.nowcoder.net/acm/contest/78/E 已知有一个n+1个数的数列,对于给定的A0和An ,当i满足当1<=i<=n-1时有 ...
- 《Linux内核分析》课程第八周学习总结
姓名:何伟钦 学号:20135223 ( *原创作品转载请注明出处*) ( 学习课程:<Linux内核分析>MOOC课程http://mooc.study.163.com/course/U ...
- R和python语言如何求平均值,中位数和众数
均值是通过取数值的总和并除以数据序列中的值的数量来计算. R语言平均值公式: mean(x, trim = 0, na.rm = FALSE, ...)#x - 是输入向量.trim - 用于从排序的 ...
- org.springframework.beans.factory.BeanDefinitionStoreException: IOException parsing XML document from class path resource [applicationContext.xml]; nested exception is java.io.FileNotFoundException: c
//这个是 配置文件放错了地方 org.springframework.beans.factory.BeanDefinitionStoreException: IOException parsing ...