2018.09.27 hdu5564Clarke and digits(数位dp+矩阵快速幂)
传送门
好题啊。
我只会写l,rl,rl,r都很小的情况(然而题上并没有这种数据范围)。
但这个dp转移式子可以借鉴。
我们用f[i][j][k]f[i][j][k]f[i][j][k]表示当前在第i位,模7余j,当前位是k。
显然有f[i+1][([j∗10+l)f[i+1][([j*10+l)f[i+1][([j∗10+l)%7][l]+=f[i][j][k]7][l]+=f[i][j][k]7][l]+=f[i][j][k]。
但是i上限1e91e91e9,直接做会凉。
于是我们构造矩阵来优化。
我们的矩阵只跟第二,三维有关。
准确的说,令t=10∗j+kt=10*j+kt=10∗j+k,t′=((j∗10+l)t'=((j*10+l)t′=((j∗10+l)%7)∗10+l7)*10+l7)∗10+l。
那么这个t->t’是合法的,就对矩阵对应的val[t][t′]val[t][t']val[t][t′]有1的贡献。
然后会发现需要记录矩阵的前缀和。
这个我们在矩阵最后添一列1来算就行了。
然后就快速幂弄一弄就行了。
代码:
#include<bits/stdc++.h>
#define ll long long
#define N 72
#define mod 1000000007
using namespace std;
int T,K;
ll L,R;
struct Matrix{ll val[N][N];Matrix(){memset(val,0,sizeof(val));}};
inline Matrix operator*(Matrix a,Matrix b){
Matrix c;
for(int i=0;i<N;++i)for(int j=0;j<N;++j)for(int k=0;k<N;++k)
(c.val[i][j]+=a.val[i][k]*b.val[k][j]%mod)%=mod;
return c;
}
inline Matrix operator^(Matrix A,int p){
Matrix ret;
for(int i=0;i<N;++i)ret.val[i][i]=1;
while(p){
if(p&1)ret=ret*A;
p>>=1,A=A*A;
}
return ret;
}
inline int calc(int i,int j){return i*10+j;}
int main(){
scanf("%d",&T);
while(T--){
scanf("%lld%lld%d",&L,&R,&K);
Matrix A,B;
for(int i=0;i<7;++i)for(int j=0;j<10;++j)for(int k=0;k<10;++k)
if(j+k!=K)++B.val[calc(i,j)][calc((i*10+k)%7,k)];
for(int i=0;i<10;++i)++B.val[i][N-1];
++B.val[N-1][N-1];
for(int i=1;i<10;++i)++A.val[0][calc(i%7,i)];
Matrix tmp1=A*(B^R),tmp2=A*(B^(L-1));
ll ans=((tmp1.val[0][N-1]-tmp2.val[0][N-1])%mod+mod)%mod;
printf("%lld\n",ans);
}
return 0;
}
2018.09.27 hdu5564Clarke and digits(数位dp+矩阵快速幂)的更多相关文章
- hdu5564--Clarke and digits(数位dp+矩阵快速幂)
Clarke and digits 问题描述 克拉克是一名人格分裂患者.某一天,克拉克变成了一个研究人员,在研究数字. 他想知道在所有长度在[l,r]之间的能被7整除且相邻数位之和不为k的正整数有多少 ...
- HUST 1569(Burnside定理+容斥+数位dp+矩阵快速幂)
传送门:Gift 题意:由n(n<=1e9)个珍珠构成的项链,珍珠包含幸运数字(有且仅由4或7组成),取区间[L,R]内的数字,相邻的数字不能相同,且旋转得到的相同的数列为一种,为最终能构成多少 ...
- BZOJ3329 Xorequ(数位dp+矩阵快速幂)
显然当x中没有相邻的1时该式成立,看起来这也是必要的. 于是对于第一问,数位dp即可.第二问写出dp式子后发现就是斐波拉契数列,矩阵快速幂即可. #include<iostream> #i ...
- BZOJ 3329 Xorequ:数位dp + 矩阵快速幂
传送门 题意 现有如下方程:$ x \oplus 3x = 2x $ 其中 $ \oplus $ 表示按位异或. 共 $ T $ 组数据,每组数据给定正整数 $ n $,任务如下: 求出小于等于 $ ...
- BZOJ3329: Xorequ(二进制数位dp 矩阵快速幂)
题意 题目链接 Sol 挺套路的一道题 首先把式子移一下项 \(x \oplus 2x = 3x\) 有一件显然的事情:\(a \oplus b \leqslant c\) 又因为\(a \oplus ...
- bnuoj 34985 Elegant String DP+矩阵快速幂
题目链接:http://acm.bnu.edu.cn/bnuoj/problem_show.php?pid=34985 We define a kind of strings as elegant s ...
- HDU 5434 Peace small elephant 状压dp+矩阵快速幂
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5434 Peace small elephant Accepts: 38 Submissions: ...
- 【BZOJ】2004: [Hnoi2010]Bus 公交线路 状压DP+矩阵快速幂
[题意]n个点等距排列在长度为n-1的直线上,初始点1~k都有一辆公车,每辆公车都需要一些停靠点,每个点至多只能被一辆公车停靠,且每辆公车相邻两个停靠点的距离至多为p,所有公车最后会停在n-k+1~n ...
- 【BZOJ】4861: [Beijing2017]魔法咒语 AC自动机+DP+矩阵快速幂
[题意]给定n个原串和m个禁忌串,要求用原串集合能拼出的不含禁忌串且长度为L的串的数量.(60%)n,m<=50,L<=100.(40%)原串长度为1或2,L<=10^18. [算法 ...
随机推荐
- IOUtils总结
参考:https://www.cnblogs.com/xing901022/p/5978989.html 常用的静态变量 在IOUtils中还是有很多常用的一些变量的,比如换行符等等 public s ...
- ADODB 调用存储过程
追加参数法调用存储过程 追加参数通过CreateParameter方法,用来指定属性创建新的Parameter对象.具体语法如下: Set parameter = command.CreatePara ...
- nginx的配置文件解析
worker_processes ;#工作进程的个数,一般与计算机的cpu核数一致 events { worker_connections ;#单个进程最大连接数(最大连接数=连接数*进程数) } h ...
- delphi修改QQ快捷方式的目标地址达到在启动QQ的同时也能运行自己想要启动的EXE可执行文件
delphi修改QQ快捷方式的目标地址达到在启动QQ的同时也能运行自己想要启动的EXE可执行文件. 直接上代码,自已体会 !! Unit1.pas代码如下: unit Unit1; interface ...
- servlet类第二篇
1servlet的生命周期是什么? 服务器启动时(web.xml中配置load-on-startup=1,默认为0)或者第一次请求该servlet时,就会初始化一个Servlet对象,也就是会执行初始 ...
- (转)游戏引擎中三大及时光照渲染方法介绍(以unity3d为例)
重要:在目前市面上常见的游戏引擎中,主要采用以下三种灯光实现方式: 顶点照明渲染路径细节 Vertex Lit Rendering Path Details 正向渲染路径细节 Forward Rend ...
- GBDT 将子树结果当成lr输出
http://scikit-learn.org/stable/auto_examples/ensemble/plot_feature_transformation.html#example-ensem ...
- Python运维开发基础10-函数基础
一,函数的非固定参数 1.1 默认参数 在定义形参的时候,提前给形参赋一个固定的值. #代码演示: def test(x,y=2): #形参里有一个默认参数 print (x) print (y) t ...
- Spring ThreadPoolTaskExecutor
1. ThreadPoolTaskExecutor配置 1 <!-- spring thread pool executor --> 2 <bean id="taskExe ...
- padding-top和margin-top的区别
学习前端知识,对我们查找页面元素很有帮助,而且自己在原公司时,有参与一个QA系统,自己去设计了这个产品,出了原型图,同时设计了几个页面,希望通过做这个产品提高自己的技术,可是因为离职,所以计划搁浅了, ...