bzoj3992【SDOI2015】序列统计
3992: [SDOI2015]序列统计
Time Limit: 30 Sec Memory Limit: 128 MB
Submit: 673 Solved: 327
[Submit][Status][Discuss]
Description
Input
一行。四个整数。N、M、x、|S|,当中|S|为集合S中元素个数。
第二行,|S|个整数,表示集合S中的全部元素。
Output
一行,一个整数,表示你求出的种类数mod
1004535809的值。
Sample Input
1 2
Sample Output
HINT
Source
NTT第一题
首先能够发现1004535809=479*2^21+1,并且是一个质数,所以能够用NTT解决。
用f[i][j]表示i个数模m等于g^j的方案数(i为2的整数次幂,g为m的原根)。则f[i][j]=∑f[i/2][k]*f[i/2][j-k]。
这样就成了卷积形式,NTT搞定。
但n并不一定是2的整数次幂,这里就要用到高速幂的思想(详见代码)。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
#define ll long long
#define maxn 40000
#define mod 1004535809
using namespace std;
int n,m,num,s,mg,g,bit,inv;
int a[maxn],c[maxn],A[maxn],B[maxn],ind[maxn],rev[maxn];
inline int read()
{
int x=0,f=1;char ch=getchar();
while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
inline ll power(ll x,int y,int p)
{
ll ret=1;
for(;y;y>>=1,x=x*x%p) if (y&1) ret=ret*x%p;
return ret;
}
inline bool get_order(int x,int m)
{
int lim=sqrt(m),phi=m-1;
F(i,1,lim) if (phi%i==0)
{
if (power(x,i,m)==1){if (i!=m-1) return false;}
if (power(x,phi/i,m)==1){if (phi/i!=m-1) return false;}
}
return true;
}
inline int get_primitive_root(int x)
{
F(i,2,x) if (get_order(i,x)) return i;
}
inline void ntt(int *a,int flag)
{
F(i,0,(1<<bit)-1) if (rev[i]>i) swap(a[i],a[rev[i]]);
F(i,1,bit)
{
int y=(1ll*flag*(mod-1)/(1<<i)+mod-1)%(mod-1);
ll wn=power(g,y,mod);
for(int j=0;(j<1<<bit);j+=(1<<i))
{
ll w=1;
F(k,j,j+(1<<(i-1))-1)
{
int u=a[k],v=w*a[k+(1<<(i-1))]%mod;
a[k]=(u+v)%mod;
a[k+(1<<(i-1))]=((u-v)%mod+mod)%mod;
w=w*wn%mod;
}
}
}
if (flag<0) F(i,0,(1<<bit)-1) a[i]=1ll*a[i]*inv%mod;
}
inline void convol(int *a,int *b)
{
int len=1<<bit;
F(i,0,len-1) c[i]=b[i];
ntt(a,1);ntt(c,1);
F(i,0,len-1) a[i]=(ll)a[i]*c[i]%mod;
ntt(a,-1);
for(int i=m-1,j=0;i<len;i++,j++) a[j]=(a[j]+a[i])%mod,a[i]=0;
}
int main()
{
n=read();m=read();num=read();s=read();
F(i,1,s) a[i]=read();
mg=get_primitive_root(m);
g=get_primitive_root(mod);
int tmp=1;
F(i,0,m-2){ind[tmp]=i;tmp=tmp*mg%m;}
for(bit=0;(1<<bit)<(m-1)*2;bit++);
inv=power(1<<bit,mod-2,mod);
F(i,0,(1<<bit)-1) rev[i]=rev[i>>1]>>1|((i&1)<<(bit-1));
A[0]=1;
F(i,1,s) if (a[i]) B[ind[a[i]]]=1;
for(;n;convol(B,B),n>>=1) if (n&1) convol(A,B);
printf("%d\n",A[ind[num]]);
return 0;
}
bzoj3992【SDOI2015】序列统计的更多相关文章
- [BZOJ3992][SDOI2015]序列统计(DP+原根+NTT)
3992: [SDOI2015]序列统计 Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 1888 Solved: 898[Submit][Statu ...
- BZOJ3992: [SDOI2015]序列统计(NTT 原根 生成函数)
题意 题目链接 给出大小为\(S\)的集合,从中选出\(N\)个数,满足他们的乘积\(\% M = X\)的方案数 Sol 神仙题Orz 首先不难列出最裸的dp方程,设\(f[i][j]\)表示选了\ ...
- BZOJ3992: [SDOI2015]序列统计
Description 小C有一个集合S,里面的元素都是小于M的非负整数.他用程序编写了一个数列生成器,可以生成一个长度为N的数列,数列中的每个数都属于集合S. 小C用这个生成器生成了许多这样的数列. ...
- 2018.12.31 bzoj3992: [SDOI2015]序列统计(生成函数+ntt+快速幂)
传送门 生成函数简单题. 题意:给出一个集合A={a1,a2,...as}A=\{a_1,a_2,...a_s\}A={a1,a2,...as},所有数都在[0,m−1][0,m-1][0,m− ...
- 【动态规划】bzoj3992 [Sdoi2015]序列统计 10分
#include<cstdio> using namespace std; #define MOD 1004535809 int a[8001],f[1001][101],n,m,x,S; ...
- 【NTT】bzoj3992: [SDOI2015]序列统计
板子题都差点不会了 Description 小C有一个集合S,里面的元素都是小于M的非负整数.他用程序编写了一个数列生成器,可以生成一个长度为N的数 列,数列中的每个数都属于集合S.小C用这个生成器生 ...
- BZOJ3992 [SDOI2015]序列统计 【生成函数 + 多项式快速幂】
题目 小C有一个集合S,里面的元素都是小于M的非负整数.他用程序编写了一个数列生成器,可以生成一个长度为N的数 列,数列中的每个数都属于集合S.小C用这个生成器生成了许多这样的数列.但是小C有一个问题 ...
- 【BZOJ3992】[SDOI2015]序列统计 NTT+多项式快速幂
[BZOJ3992][SDOI2015]序列统计 Description 小C有一个集合S,里面的元素都是小于M的非负整数.他用程序编写了一个数列生成器,可以生成一个长度为N的数列,数列中的每个数都属 ...
- BZOJ 3992: [SDOI2015]序列统计 [快速数论变换 生成函数 离散对数]
3992: [SDOI2015]序列统计 Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 1017 Solved: 466[Submit][Statu ...
- [SDOI2015]序列统计
[SDOI2015]序列统计 标签: NTT 快速幂 Description 给你一个模m意义下的数集,需要用这个数集生成一个数列,使得这个数列在的乘积为x. 问方案数模\(1004535809\). ...
随机推荐
- Ubuntu 16.04桌面版GUI网络配置工具NetworkManager的命令行工具nm-tool无法使用的问题
说明: 1.Ubuntu中分桌面版和服务器版,而这两个版本在网络管理方面使用的工具都不一样,尤其是在桌面版,使用了NetworkManager进行管理. 2.服务器版使用的是命令行配置,而桌面版包含了 ...
- 【gcc】warning: control reaches end of non-void function
用gcc编译一个C程序的时候出现这样的警告: warning: control reaches end of non-void function 它的意思是:控制到达非void函数的结尾.就是说你的一 ...
- 【Linux】Centos7安装之后,双系统的情况下,怎么能在CentOS7下访问Windows的磁盘
想要在CentOS7下访问Windows的NTFS格式的磁盘,需要在Linux下下载ntfs-3g步骤1: 进入root用户下,使用yum命令下载ntfs-3g.[前提是已经添加了常用源:http:/ ...
- Go -- 升级go版本
先卸载go的旧版本, 参考卸载go; 然后用brew安装, 如果之前用brew安装的go, 可直接brew update go, 否则, 安装go新版: 执行一下 export PATH=$PATH: ...
- Go -- 卸载 Go
直接删除usr下的go文件夹即可 命令行: cd / 1 然后: open usr 1 删除go文件夹
- PHP模板解析入门
学习:李炎恢PHP视频第二季 模板引擎的特点: 1.鼓励分离:让更个系统的可读性和维护性得到提高. 2.促进分工:使得程序员和美工去专心处理自己的设计. 3.比PHP更容易解析:编译文件和缓存文件加载 ...
- scala函数返回值
1.使用returndef functionName ([参数列表]) : [return type] = { function body return [expr] } 2.直接把返回值写在最后: ...
- 对tensorflow 中的attention encoder-decoder模型调试分析
#-*-coding:utf8-*- __author = "buyizhiyou" __date = "2017-11-21" import random, ...
- TCP/IP协议组随笔
原文:https://my.oschina.net/xianggao/blog/654677 IP层负责网络主机的定位,数据传输的路由,由IP地址可以唯一的确定Internet上的一台主机. TCP层 ...
- Java高级特性—JVM
1).java监控工具使用 jconsole是一种集成了上面所有命令功能的可视化工具,可以分析jvm的内存使用情况和线程等信息 visualvm 提供了和jconsole的功能类似,提供了一大堆的插件 ...