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\). ...
随机推荐
- Bluetooth篇 开发实例之十 官网的Bluetooth Chat sample app.
运行的时候,会报错: java.lang.NullPointerException: Attempt to invoke virtual method 'void android.app.Action ...
- sql server 存储过程中使用变量表,临时表的分析(续)
最近,我有一朋友,对我说他的数据库中的很多存储过程,执行都是超时.让我替他看看是什么原因.我一看,原来他的存储过程中用了很多的临时表与变量表.于是我跟他说过犹不及. 在存储过程中使用临时表或变量表,使 ...
- delphi 浮点数转换成十六进制字符串的方法
我们在研究封包技术时,经常会碰到将浮点数转换成十六进制形式.比如在游戏中人物的座标,经常就用浮点数来表示.怎么将浮点数转换成十六进制字符串形式呢?下面我将写出其在DELPHI中的方法. 先 ...
- java内存缓存,节省内存
缓存的对象 这个问题就是我们上面提到的极端情况,在Java中,会对-128到127的Integer对象进行缓存,当创建新的Integer对象时,如果符合这个这个范围,并且已有存在的相同值的对象,则返回 ...
- django开发环境部署之pip、virtualenv、virtualenvwrapper
step1:安装pip 在python中可以使用easy_install和pip安装python拓展但推荐使用pip Don't use easy_install, unless you like s ...
- 2017.6.29 移除再导入maven module到IDEA中时提示: Unable to proceed. Nothing found to import.
解决办法来自:https://stackoverflow.com/questions/18278016/re-importing-modules-into-intellij 场景: 将其中一个modu ...
- AlamoFireDemo
// // ViewController.swift // AlamFireDemo // // import UIKit import Alamofire class ViewController: ...
- Redis与Reactor模式
Redis与Reactor模式 Jan 9, 2016 近期看了Redis的设计与实现,这本书写的还不错,看完后对Redis的理解有非常大的帮助. 另外,作者整理了一份Redis源代码凝视,大家能够c ...
- 【转】Intellij IDEA常用配置详解
1. IDEA内存优化 先看看你机器本身的配置而配置. \IntelliJ IDEA 8\bin\idea.exe.vmoptions -------------------------------- ...
- 【SpringMVC学习08】SpringMVC中实现文件上传
之前有写过一篇struts2实现的文件上传,这一篇博文主要来总结下springmvc实现文件上传的步骤.首先来看一下单个文件的上传,然后再来总结下多个文件上传. 1. 环境准备 springmvc上传 ...