2119: 【BZOJ3992】【SDOI2015】序列统计

Description

小\(C\)有一个集合\(S\),里面的元素都是小于\(M\)的非负整数。

他用程序编写了一个数列生成器,可以生成一个长度为\(N\)的数列,数列中的每个数都属于集合\(S\)。

小\(C\)用这个生成器生成了许多这样的数列。但是小\(C\)有一个问题需要你的帮助:给定整数\(x\),求所有可以生成出的,且满足数列中所有数的乘积\(\bmod M\)的值等于\(x\)的不同的数列的有多少个。

小\(C\)认为,两个数列\(\{A_i\}\)和\(\{B_i\}\)不同,当且仅当至少存在一个整数\(i\),满足\(A_i≠B_i\)。另外,小\(C\)认为这个问题的答案可能很大,因此他只需要你帮助他求出答案\(\bmod 1004535809\)的值就可以了。

Input

一行,四个整数,\(N\)、\(M\)、\(x\)、\(|S|\),其中\(|S|\)为集合\(S\)中元素个数。

第二行,\(|S|\)个整数,表示集合\(S\)中的所有元素。

Output

一行,一个整数,表示你求出的种类数\(\bmod 1004535809\)的值。

HINT

对于\(10\%\)的数据,\(1≤N≤1000\);

对于\(30\%\)的数据,\(3≤M≤100\);

对于\(60\%\)的数据,\(3≤M≤800\);

对于全部的数据,\(1≤N≤10^9\),\(3≤M≤8000\),\(M\)为质数,\(1≤x≤M−1\),输入数据保证集合\(S\)中元素不重复


先吐槽一波部分分吧...

明显有一个矩阵快速幂的做法,然后通过不了60pts,请问这个点的分怎么拿呢?

\(dp_{i,j}\)前\(i\)个数乘积为\(j\)的方案

\[dp_{i,j}=\sum dp_{i-1,ab\% m=j}
\]


思路:利用原根转换问题,然后生成函数用多项式快速幂求一下就行了。

对原根,若\(g\)是\(m\)的原根,则\(g^0,g^1,\dots,g^{\varphi(m)-1}\)遍历\(\bmod m\)的最小剩余系,就是它们有个一一对应关系

关于求原根,把\(\varphi(m)-1\)分解质因数得到每个质因子\(p_i\),然后从小到大枚举原根\(g\),检验每个质因子是否有\(g^{\frac{\varphi(m)-1}{p_i}}\equiv 1 \pmod m\),如果所有质因子都没有,\(g\)就是原根,因为原根非常的小,所以这个复杂度是正确的。

然后把前面\(ab\%m=j\)换掉,就成了\((a’+b')\%(m-1)=j'\)

然后你构造生成函数\(f_i=\sum is_ix^{i-1}\),\(is_i\)代表最开始转换后这个值存不存在。

然后发现结果就是这个\(\tt f^n\),写一个\(NTT\)快速幂就成了,复杂度\(O(m\log m\log n)\)


Code:

#include <cstdio>
#include <algorithm>
const int N=(1<<14)+10;
const int mod=1004535809,Gi=334845270;
#define add(x,y,p) ((x+y)%p)
#define mul(x,y,p) (1ll*(x)*(y)%p)
int qp(int d,int k,int p){int f=1;while(k){if(k&1) f=mul(f,d,p);d=mul(d,d,p),k>>=1;}return f;}
int A[N],B[N],F[N],D[N],n,m,turn[N],to[N],len=1;
int GetG(int x)
{
int s[1010]={0},phi=x-1,t=x-1;
for(int i=2;i*i<=t;i++)
{
if(t%i==0)
{
s[++s[0]]=i;
while(t%i==0) t/=i;
}
}
if(t!=1) s[++s[0]]=t;
for(int i=2;i<x;i++)
{
int flag=1;
for(int j=1;j<=s[0];j++)
if(qp(i,phi/s[j],x)==1)
{
flag=0;
break;
}
if(flag) return i;
}
return -1;
}
void NTT(int *a,int typ)
{
for(int i=0;i<len;i++) if(i<turn[i]) std::swap(a[i],a[turn[i]]);
for(int le=1;le<len;le<<=1)
{
int wn=qp(typ?3:Gi,(mod-1)/(le<<1),mod);
for(int p=0;p<len;p+=le<<1)
{
int w=1;
for(int i=p;i<p+le;i++,w=mul(w,wn,mod))
{
int tx=a[i],ty=mul(w,a[i+le],mod);
a[i]=add(tx,ty,mod);
a[i+le]=add(tx,mod-ty,mod);
}
}
}
if(!typ)
{
int inv=qp(len,mod-2,mod);
for(int i=0;i<len;i++) a[i]=mul(a[i],inv,mod);
}
}
void Rev()
{
int L=-1;for(int i=1;i<len;i<<=1) ++L;
for(int i=1;i<len;i++) turn[i]=turn[i>>1]>>1|(i&1)<<L;
}
void polymul(int *a,int *b)
{
for(int i=0;i<len;i++) A[i]=a[i],B[i]=b[i],a[i]=0;
NTT(A,1),NTT(B,1);
for(int i=0;i<len;i++) A[i]=mul(A[i],B[i],mod);
NTT(A,0);
for(int i=0;i<len;i++) a[i%(m-1)]=add(a[i%(m-1)],A[i],mod);
}
int main()
{
int s,x;
scanf("%d%d%d%d",&n,&m,&x,&s);
int g=GetG(m);
while(len<=(m<<1))len<<=1;Rev();
for(int f=1,i=0;i<m-1;i++,f=mul(f,g,m))
to[f]=i;
for(int x,i=1;i<=s;i++)
{
scanf("%d",&x);
if(x) F[to[x]]=1,D[to[x]]=1;
}
--n;
while(n)
{
if(n&1) polymul(F,D);
polymul(D,D);
n>>=1;
}
printf("%d\n",F[to[x]]);
return 0;
}

2018.12.17

【SDOI2015】序列统计 解题报告的更多相关文章

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

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

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

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

  3. [SDOI2015]序列统计

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

  4. 3992: [SDOI2015]序列统计

    3992: [SDOI2015]序列统计 链接 分析: 给定一个集和s,求多少个长度为n的序列,满足序列中每个数都属于s,并且所有数的乘积模m等于x. 设$f=\sum\limits_{i=0}^{n ...

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

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

  6. 【LG3321】[SDOI2015]序列统计

    [LG3321][SDOI2015]序列统计 题面 洛谷 题解 前置芝士:原根 我们先看一下对于一个数\(p\),它的原根\(g\)有什么性质(好像就是定义): \(g^0\%p,g^1\%p,g^2 ...

  7. [BZOJ3992][SDOI2015]序列统计(DP+原根+NTT)

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

  8. 【BZOJ3992】[SDOI2015]序列统计 NTT+多项式快速幂

    [BZOJ3992][SDOI2015]序列统计 Description 小C有一个集合S,里面的元素都是小于M的非负整数.他用程序编写了一个数列生成器,可以生成一个长度为N的数列,数列中的每个数都属 ...

  9. BZOJ 3992: [SDOI2015]序列统计 NTT+快速幂

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

随机推荐

  1. 搭建 Digital Ocean 服务器

    作为程序员,搭建属于自己的服务器可以部署自己的网站,可以配置代理***,了解国际前沿科技:一些高性能的服务器还可以完成一些云计算.深度学习模型的训练任务. DigitalOcean 服务启动 注册 D ...

  2. 从一个简单的寻路问题深入Q-learning

    这第一篇随笔实际上在我的科学网博客上是首发,我重新拿到博客园再发一次是希望以此作为我学习Q-learning的一个新的开始.以后这边主技术,科学网博客主理论.我也会将科学网那边技术类的文章转过来的.希 ...

  3. java web 3.1-web.xml文件配置

    <?xml version="1.0" encoding="UTF-8" ?> <web-app xmlns:xsi="http:/ ...

  4. AutoResetEvent 方法名称设计缺陷

    这个类和方法,让人乍一读是读不明白的.不能通过方法名称明白其含义.所以它的方法名称设计是欠考虑. 应该类似于这样: public static class MyAutoResetEvent { pub ...

  5. php 常用英语小汇

    bstract抽象的 -挨伯丝拽克特 access存取.访问 -挨克色丝 account账户 -厄靠恩特 action动作 -爱克身 activate激活 -爱克特维特 active活动的 -爱克得 ...

  6. webpack入门指南-step02

    webpack 安装 1)安装前的准备:webpack是基于node环境的项目,所以使用前必须先安装node和npm. 在安装 Webpack 前,你本地环境需要支持 node.js.如果电脑没有装过 ...

  7. 团队冲刺——Three

    第三天计划: 季方:学习爬虫的操作,以便后续功能实现: 司宇航:对当天实现的功能进行总的测试: 王金萱:数据库内数据的增删改查以及查看团队博客界面的实现: 马佳慧:学习css初步,进行页面绘制: 第二 ...

  8. c# using的作用

    using 关键字有两个主要用途:   (一).作为指令,用于为命名空间创建别名或导入其他命名空间中定义的类型.   (二).作为语句,用于定义一个范围,在此范围的末尾将释放对象. using指令   ...

  9. 课堂练习 psp表

    项目计划总结表: 日期 编程 完善程序 测试程序 参考资料 日总结 3.20 18:00---19:30       1.5 3.21   9:30----10:00 10:00---10:30   ...

  10. 【Coursera】应用机器学习的建议

    偏差方差权衡 使用较小的神经网络,类似于参数较少的情况,容易导致高偏差和欠拟合,但计算代价较小使用较大的神经网络,类似于参数较多的情况,容易导致高方差和过拟合,虽然计算代价比较大,但是可以通过归一化手 ...