Description

题库链接

给出集合 \(S\) ,元素都是小于 \(M\) 的非负整数。问能够生成出多少个长度为 \(N\) 的数列 \(A\) ,数列中的每个数都属于集合 \(S\) ,并且

\[\prod_{i=1}^N A_i\equiv x \pmod{M}\]

答案对 \(1004535809\) 取模。

\(1\leq N\leq 10^9,3\leq M\leq 8000, M 为质数,0\leq x\leq M-1\)

Solution

显然能够得到 \(DP\) 的解法:令 \(f_{i,j}\) 为生成序列长度为 \(i\) 时,乘积在模 \(M\) 意义下为 \(j\) 的方案数。

显然 \(f_{i,j}\rightarrow f_{i+1,(j\times w)\mod M},w\in S\) 。

但 \(n\leq 10^9\) 显然不能递推。考虑优化。

由于乘法不太好搞,我们试着换种思路,我们不妨将集合内数取 \(\log\) 。那么 \(f_{i,\log j}\rightarrow f_{i+1,\log j+\log w},w\in S\) 。

但实数域上确实不好做,考虑取离散对数。由费马小定理,它是以 \(M-1\) 为周期的,那么只要 \(\text{NTT}\) 优化,加上快速幂。对模意义外的数讨论即可。

Code

#include <bits/stdc++.h>
using namespace std;
const int yzh = 1004535809;
const int N = 8000*4; int n, m, x, s, G, lg[N+5], a, len, L, R[N+5];
int A[N+5]; int quick_pow(int a, int b, int yzh) {
int ans = 1;
while (b) {
if (b&1) ans = 1ll*a*ans%yzh;
b >>= 1, a = 1ll*a*a%yzh;
}
return ans;
}
void get_G() {
int prime[N+5], tot = 0, x = m-1;
for (int i = 2, lim = sqrt(x)+1; i <= lim; i++)
if (x%i == 0) {
prime[++tot] = i;
while (x%i == 0) x /= i;
}
if (x != 1) prime[++tot] = x;
for (int i = 2; true; i++) {
int flag = 1;
for (int j = 1; j <= tot; j++)
if (quick_pow(i, (m-1)/prime[j], m) == 1) {
flag = 0; break;
}
if (flag == 1) {G = i; break; }
}
for (int i = 1, g = G; i < m; i++, g = 1ll*g*G%m) lg[g] = i;
}
void NTT(int *A, int o) {
for (int i = 0; i < len; i++) if (i < R[i]) swap(A[i], A[R[i]]);
for (int i = 1; i < len; i <<= 1) {
int gn = quick_pow(3, (yzh-1)/(i<<1), yzh), x, y;
if (o == -1) gn = quick_pow(gn, yzh-2, yzh);
for (int j = 0; j < len; j += (i<<1)) {
int g = 1;
for (int k = 0; k < i; k++, g = 1ll*g*gn%yzh) {
x = A[j+k], y = 1ll*g*A[j+k+i]%yzh;
A[j+k] = (x+y)%yzh, A[j+k+i] = (x-y+yzh)%yzh;
}
}
}
if (o == 1) return;
for (int i = 0, inv = quick_pow(len, yzh-2, yzh); i < len; i++)
A[i] = 1ll*A[i]*inv%yzh;
for (int i = m; i < len; i++) (A[i%(m-1) ? i%(m-1) : m-1] += A[i]) %= yzh, A[i] = 0;
}
void NTTpow(int *A, int b) {
int ans[N+5] = {0}; ans[0] = 1;
while (b) {
NTT(A, 1);
if (b&1) {
NTT(ans, 1);
for (int i = 0; i < len; i++) ans[i] = 1ll*ans[i]*A[i]%yzh;
NTT(ans, -1);
}
for (int i = 0; i < len; i++) A[i] = 1ll*A[i]*A[i]%yzh;
NTT(A, -1); b >>= 1;
}
for (int i = 0; i < len; i++) A[i] = ans[i];
}
void work() {
scanf("%d%d%d%d", &n, &m, &x, &s); get_G();
for (int i = 1; i <= s; i++) {scanf("%d", &a); ++A[lg[a]]; }
A[0] = 0;
for (len = 1; len <= (m<<1); len <<= 1) ++L;
for (int i = 0; i < len; i++) R[i] = (R[i>>1]>>1)|((i&1)<<(L-1));
NTTpow(A, n); printf("%d\n", A[lg[x]]);
}
int main() {work(); return 0; }

[SDOI 2015]序列统计的更多相关文章

  1. [BZOJ 3992] [SDOI 2015] 序列统计(DP+原根+NTT)

    [BZOJ 3992] [SDOI 2015] 序列统计(DP+原根+NTT) 题面 小C有一个集合S,里面的元素都是小于质数M的非负整数.他用程序编写了一个数列生成器,可以生成一个长度为N的数列,数 ...

  2. BZOJ 3992 [SDOI 2015] 序列统计 解题报告

    这个题最暴力的搞法就是这样的: 设 $Dp[i][j]$ 为前 $i$ 个数乘积为 $j$ 的方案数. 转移的话就不多说了哈... 当前复杂度 $O(nm^2)$ 注意到,$M$ 是个质数,就说明 $ ...

  3. [BZOJ 3992] [SDOI 2015] 序列统计

    Description 传送门 Solution [一] 设 \(f[i][j]\) 表示前 \(i\) 个数的乘积在模 \(p\) 意义下等于 \(j\) 的方案数,有 \[ f[i][j]=\su ...

  4. [BZOJ 3992][SDOI2015]序列统计

    3992: [SDOI2015]序列统计 Time Limit: 30 Sec  Memory Limit: 128 MBSubmit: 2275  Solved: 1090[Submit][Stat ...

  5. Mobius反演与积性函数前缀和演学习笔记 BZOJ 4176 Lucas的数论 SDOI 2015 约数个数和

    下文中所有讨论都在数论函数范围内开展. 数论函数指的是定义域为正整数域, 且值域为复数域的函数. 数论意义下的和式处理技巧 因子 \[ \sum_{d | n} a_d = \sum_{d | n} ...

  6. Bzoj 4403: 序列统计 Lucas定理,组合数学,数论

    4403: 序列统计 Time Limit: 3 Sec  Memory Limit: 128 MBSubmit: 328  Solved: 162[Submit][Status][Discuss] ...

  7. BZOJ4403 序列统计—Lucas你好

    绝对是全网写的最详细的一篇题解  题目:序列统计 代码难度:简单 思维难度:提高+-省选 讲下题面:给定三个正整数N.L和R,统计长度在1到N之间,元素大小都在L到R之间的单调不降序列的数量.输出答案 ...

  8. BZOJ 3992: [SDOI2015]序列统计 [快速数论变换 生成函数 离散对数]

    3992: [SDOI2015]序列统计 Time Limit: 30 Sec  Memory Limit: 128 MBSubmit: 1017  Solved: 466[Submit][Statu ...

  9. [SDOI2015]序列统计

    [SDOI2015]序列统计 标签: NTT 快速幂 Description 给你一个模m意义下的数集,需要用这个数集生成一个数列,使得这个数列在的乘积为x. 问方案数模\(1004535809\). ...

随机推荐

  1. 用Python满足满足自己的“小虚荣”

    首先声明,学习这个只是为了好玩,只是为了好玩,并不是想用这个弄虚作假,做一些不好的事情!一心想做技术人,自制自治! 我们有时候发布一篇日志,或者是一篇博文,总希望自己的浏览量能高点,这样看起来也倍有面 ...

  2. centos7 下通过nginx+uwsgi部署django应用

    1. 安装python3.6 1. 获取 wget https://www.python.org/ftp/python/3.6.2/Python-3.6.2.tgz tar -xzvf Python- ...

  3. Beta集合

    Beta冲刺day1 Beta冲刺day2 Beta冲刺day3 Beta冲刺day4 Beta冲刺day5 Beta冲刺day6 Beta冲刺day7 测试总结 总结合集 Beta预备

  4. 201421123042 《Java程序设计》第9周学习总结

    1. 本周学习总结 1.1 以你喜欢的方式(思维导图或其他)归纳总结集合与泛型相关内容. 2. 书面作业 本次作业题集集合 1. List中指定元素的删除(题集题目) 1.1 实验总结.并回答:列举至 ...

  5. selenium的Python使用(一)浏览器驱动的安装及使用

    一.selenium的安装 直接使用pip进行安装 pip install selenium    #(安装最新版本) pip install selenium==3.6.0   #(安装指定版本) ...

  6. python生成单词壁纸

    1.首先上结果: 其实就是一段简单的代码.加上英语单词表加上几张背景图生成许多类似的图片再设置成桌面背景,十分钟一换.有心的人闲的时候随手就能换换桌面背背单词.最不济也能混个脸熟. 3.上代码 #-* ...

  7. python 字符串的方法

    capitalize() 把字符串的第一个字符改为大写 casefold() 把整个字符串的所有字符改为小写 center(width) 将字符串居中,并使用空格填充至长度 width 的新字符串 c ...

  8. Linux知识积累(3)$()和${}和$(())和(())

    $()和${}和$(())和(()) $()和${}的用法:在 bash shell 中,$( ) 与 ` ` (反引号) 都是用来做命令替换用(command substitution)的.而 $( ...

  9. python 异常 反射

    异常 反射 一.异常处理: AttributeError 试图访问一个对象没有的树形,比如foo.x,但是foo没有属性x IOError 输入/输出异常:基本上是无法打开文件 ImportError ...

  10. Github学习笔记-不定时更新

    最近在搞一个外包项目,有个非常厉害的师兄带我,他写后端,我负责前端部分,项目合作需要用到github,这也是早晚要接触,记录下项目过程中对github/git使用的学习笔记. 1.在网上看了一些教程之 ...