【bzoj4818】[Sdoi2017]序列计数 矩阵乘法
原文地址:http://www.cnblogs.com/GXZlegend/p/6825132.html
题目描述
输入
输出
样例输入
3 5 3
样例输出
33
题解
矩阵乘法
至少有1个质数的方案数=总方案数-不含质数的方案数
可以先在O(m)的时间内把1~m的质数筛出来。
然后考虑:从mod p=0,到mod p=0,可以由mod p=a和mod p=p-a两个阶段组成。
可以设f[i][j]表示从mod p=i到mod p=j的方案数,不难看出这是一个矩阵,而且自乘m次就能得到答案。
所以只要处理出2种f即可。
对于每个数可以更新所有的f[i](0≤i<n),不过这样会TLE
其实这样做没有必要,因为f都是循环出现的,只需要求f[0]即可,再据此推其余的f。
代码中,我将mod p=0当作了p来处理,其实看作0也是一样的。
#include <cstdio>
#include <cstring>
#include <algorithm>
#define mod 20170408
using namespace std;
typedef long long ll;
bool notprime[20000010];
int pri[5000010] , top;
struct matrix
{
int n , m;
ll num[110][110];
matrix()
{
n = m = 0;
memset(num , 0 , sizeof(num));
}
matrix operator*(const matrix a)const
{
matrix ans;
ans.n = n , ans.m = a.m;
int i , j , k;
for(i = 1 ; i <= ans.n ; i ++ )
for(j = 1 ; j <= ans.m ; j ++ )
for(k = 1 ; k <= m ; k ++ )
ans.num[i][j] = (ans.num[i][j] + num[i][k] * a.num[k][j]) % mod;
return ans;
}
}A , B;
matrix pow(matrix x , int y)
{
matrix ans;
int i;
ans.n = x.n , ans.m = x.m;
for(i = 1 ; i <= ans.n ; i ++ )
ans.num[i][i] = 1;
while(y)
{
if(y & 1) ans = ans * x;
x = x * x , y >>= 1;
}
return ans;
}
int main()
{
int n , m , p , i , j;
scanf("%d%d%d" , &n , &m , &p);
notprime[1] = 1;
for(i = 2 ; i <= m ; i ++ )
{
if(!notprime[i]) pri[++top] = i;
for(j = 1 ; j <= top && i * pri[j] <= m ; j ++ )
{
notprime[i * pri[j]] = 1;
if(i % pri[j] == 0) break;
}
}
A.n = A.m = B.n = B.m = p;
for(i = 1 ; i <= m ; i ++ )
{
A.num[p][(i - 1) % p + 1] ++ ;
if(notprime[i]) B.num[p][(i - 1) % p + 1] ++ ;
}
for(i = p - 1 ; i >= 1 ; i -- )
for(j = 1 ; j <= p ; j ++ )
A.num[i][j] = A.num[i + 1][j % p + 1] , B.num[i][j] = B.num[i + 1][j % p + 1];
printf("%lld\n" , (pow(A , n).num[p][p] - pow(B , n).num[p][p] + mod) % mod);
return 0;
}
【bzoj4818】[Sdoi2017]序列计数 矩阵乘法的更多相关文章
- BZOJ 4818 [Sdoi2017]序列计数 ——矩阵乘法
发现转移矩阵是一个循环矩阵. 然后循环矩阵乘以循环矩阵还是循环矩阵. 据说还有FFT并且更优的做法. 之后再看吧 #include <map> #include <cmath> ...
- luogu 3702 [SDOI2017]序列计数 矩阵乘法+容斥
现在看来这道题真的不难啊~ 正着求不好求,那就反着求:答案=总-全不是质数 这里有一个细节要特判:1不是质数,所以在算全不是质数的时候要特判1 code: #include <bits/stdc ...
- [bzoj4818][Sdoi2017]序列计数_矩阵乘法_欧拉筛
[Sdoi2017]序列计数 题目大意:https://www.lydsy.com/JudgeOnline/problem.php?id=4818. 题解: 首先列出来一个递推式子 $f[i][0]$ ...
- [Sdoi2017]序列计数 [矩阵快速幂]
[Sdoi2017]序列计数 题意:长为\(n \le 10^9\)由不超过\(m \le 2 \cdot 10^7\)的正整数构成的和为\(t\le 100\)的倍数且至少有一个质数的序列个数 总- ...
- [BZOJ 4818/LuoguP3702][SDOI2017] 序列计数 (矩阵加速DP)
题面: 传送门:https://www.lydsy.com/JudgeOnline/problem.php?id=4818 Solution 看到这道题,我们不妨先考虑一下20分怎么搞 想到暴力,本蒟 ...
- [BZOJ4818][SDOI2017]序列计数(动规+快速幂)
4818: [Sdoi2017]序列计数 Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 972 Solved: 581[Submit][Status ...
- 【BZOJ4818】【SDOI2017】序列计数 [矩阵乘法][DP]
序列计数 Time Limit: 30 Sec Memory Limit: 128 MB[Submit][Status][Discuss] Description Alice想要得到一个长度为n的序 ...
- 2019.02.11 bzoj4818: [Sdoi2017]序列计数(矩阵快速幂优化dp)
传送门 题意简述:问有多少长度为n的序列,序列中的数都是不超过m的正整数,而且这n个数的和是p的倍数,且其中至少有一个数是质数,答案对201704082017040820170408取模(n≤1e9, ...
- BZOJ4818 [SDOI2017] 序列计数 【矩阵快速幂】
题目分析: 一个很显然的同类项合并.注意到p的大小最大为100,考虑把模p意义下相同的求出来最后所有的减去没有质数的做矩阵快速幂即可. 代码: #include<bits/stdc++.h> ...
随机推荐
- Core Location :⽤用于地理定位
Core Location :⽤用于地理定位 在移动互联⽹网时代,移动app能解决⽤用户的很多⽣生活琐事,⽐比如 导航:去任意陌⽣生的地⽅方 周边:找餐馆.找酒店.找银⾏行.找电影院 在上述应⽤用中, ...
- Spring Cloud 入门Eureka -Consumer服务消费(声明式Feign)(三)
Spring Cloud Feign是一套基于Netflix Feign实现的声明式服务调用客户端.它使得编写Web服务客户端变得更加简单.我们只需要通过创建接口并用注解来配置它既可完成对Web服务接 ...
- django中的分页管理
有时,展示的对象太多,需要对他们进行分页展示,不能一页把所有的结果都展示出来吧,那样的话,哈哈,挺逗 使用Django分页器功能 从Django中导入Paginator模块(没有的话,自行下载,我是w ...
- uniqueidentifier数据类型转换
cast(id as varchar(36))
- 工作中遇到的比较奇怪的一些sql(一些子查询)
在列中进行子查询 1.在一个表中有多个员工ID,比如一个下单员工,一个修改订单的员工,可以使用在列中进行子查询,具体如下: ( SELECT staff_name FROM sp_staff_basi ...
- PHP.23-ThinkPHP框架的三种模型实例化-(D()方法与M()方法的区别)
三种模型实例化 原则上:每个数据表应对应一个模型类(Home/Model/GoodsModel.class.php --> 表tp_goods) 1.直接实例化 和实例化其他类库一样实例化模型类 ...
- 15 Django组件-中间件
中间件 中间件的概念 中间件顾名思义,是介于request与response处理之间的一道处理过程,相对比较轻量级,并且在全局上改变django的输入与输出.因为改变的是全局,所以需要谨慎实用,用不好 ...
- 15.3,redis持久化RDB与AOF
redis持久化 Redis是一种内存型数据库,一旦服务器进程退出,数据库的数据就会丢失,为了解决这个问题,Redis提供了两种持久化的方案,将内存中的数据保存到磁盘中,避免数据的丢失. RDB持久化 ...
- Dapper.Extension的基本使用
前言 上一篇随笔写了Dapper的简单的使用,这次写一下Dapper.Extension的使用,它是Dapper的简单的封装扩展,可以通过实例化的对象赋值后进行增删改的操作以及分页,但是却不能进 ...
- Android 打印方法调用堆栈
RuntimeException here = new RuntimeException("here"); here.fillInStackTrace(); Log.w(" ...