博客
关于我
[SDOI2017]序列计数
阅读量:351 次
发布时间:2019-03-04

本文共 1343 字,大约阅读时间需要 4 分钟。

题目

思路

感谢提供了思路。和代码。

基本思路,是从 d p \tt{dp} dp 中脱胎换骨:

如果用 f ( i , j ) f(i,j) f(i,j) 表示 i i i 个数字、余数为 j j j 的答案,那么对于任意 L L L 均有 f ( i , j ) = ∑ x = 0 p − 1 f ( L , x ) × f ( i − L , ( j + p − x )   m o d   p ) f(i,j)=\sum_{x=0}^{p-1}f(L,x)\times f(i-L,(j+p-x)\bmod p) f(i,j)=x=0p1f(L,x)×f(iL,(j+px)modp)

也就是枚举前 L L L 个数字的余数,与剩下的余数相加。

然后我们写一个函数 g ( L , x ) = ∑ i = 0 p − 1 f ( L , i ) x i g(L,x)=\sum_{i=0}^{p-1}f(L,i)x^i g(L,x)=i=0p1f(L,i)xi ,因为 d p \tt{dp} dp 转移非常像卷积。

于是我们就可以写成快速幂的形式。求的答案就是 g n ( 1 , x ) g^n(1,x) gn(1,x) 中的某一项。

嗯?你问我怎么处理 “有一个质数” ?正难则反,减去全部是合数的情况即可。

代码

#include
#include
#define mod 20170408using namespace std;int n,m,p,vis[20000010],cnt,pre[2000000];long long f[210],F[210],g[210],G[210];long long c[210],x=1;void work(long long *a,long long *b,long long *d) { int i,j; for(i=0; i
>n>>m>>p; F[0]=G[0]=f[1]=g[1]=1; for(i=2; i<=m; i++) { f[i%p]++; //此时的f为f[1][0--p]的值 if(!vis[i])pre[++cnt]=i; else g[i%p]++; //通过欧拉筛记录g[1][0--p]的值 for(j=1; pre[j]*i<=m; j++) { vis[pre[j]*i]=1; if(!(i%pre[j])) break; } } while(n) { if(n&1) work(F,f,F),work(G,g,G); work(f,f,f); work(g,g,g); n>>=1; } cout<<(F[0]-G[0]+mod)%mod; //特别注意这里,学长特别提醒我了的。F[0]可能比G[0]小。 return 0;}

转载地址:http://lovq.baihongyu.com/

你可能感兴趣的文章
mysql 常见问题
查看>>
MYSQL 幻读(Phantom Problem)不可重复读
查看>>
mysql 往字段后面加字符串
查看>>
mysql 快照读 幻读_innodb当前读 与 快照读 and rr级别是否真正避免了幻读
查看>>
MySQL 快速创建千万级测试数据
查看>>
mysql 快速自增假数据, 新增假数据,mysql自增假数据
查看>>
MySql 手动执行主从备份
查看>>
Mysql 批量修改四种方式效率对比(一)
查看>>
Mysql 报错 Field 'id' doesn't have a default value
查看>>
MySQL 报错:Duplicate entry 'xxx' for key 'UNIQ_XXXX'
查看>>
Mysql 拼接多个字段作为查询条件查询方法
查看>>
mysql 排序id_mysql如何按特定id排序
查看>>
Mysql 提示:Communication link failure
查看>>
mysql 插入是否成功_PDO mysql:如何知道插入是否成功
查看>>
Mysql 数据库InnoDB存储引擎中主要组件的刷新清理条件:脏页、RedoLog重做日志、Insert Buffer或ChangeBuffer、Undo Log
查看>>
mysql 数据库中 count(*),count(1),count(列名)区别和效率问题
查看>>
mysql 数据库备份及ibdata1的瘦身
查看>>
MySQL 数据库备份种类以及常用备份工具汇总
查看>>
mysql 数据库存储引擎怎么选择?快来看看性能测试吧
查看>>
MySQL 数据库操作指南:学习如何使用 Python 进行增删改查操作
查看>>