显然存在方案的数一定是L的因数,考虑对其因子预处理答案,O(1)回答。

  考虑每个质因子,设其在g中有x个,l中有y个,则要求所有选中的数该质因子个数都在[x,y]中,且存在数的质因子个数为x、y。对于后一个限制,显然可以简单地容斥,即[x,y]-[x+1,y]-[x,y-1]+[x+1,y-1],枚举这个至多是48的,这个取最大值时因子个数是28。暴力枚举数数即可。复杂度总之O(能过)。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define P 1000000007
#define ll long long
int read()
{
int x=,f=;char c=getchar();
while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
return x*f;
}
int n,g,l,u,m,prime[],cnt[][],d[],ans[],p[][],q[][],a[],tot,t,sum;
int ksm(int a,int k)
{
int s=;
for (;k;k>>=,a=1ll*a*a%P) if (k&) s=1ll*s*a%P;
return s;
}
void get(int k,int s)
{
if (k>t) {sum++;return;}
ll x=ksm(prime[k],cnt[][k]);
for (int i=cnt[][k];i<=cnt[][k];i++)
{
if (s*x<=n) get(k+,s*x);else break;
x=1ll*x*prime[k];
}
}
void build(int k,int s)
{
if (k>t)
{
p[++tot][]=s;
for (int i=;i<=t;i++) p[tot][i]=a[i];
return;
}
ll x=ksm(prime[k],cnt[][k]);
for (int i=cnt[][k];i<=cnt[][k];i++)
{
a[k]=i;
if (s*x<=n) build(k+,s*x);else break;
x=1ll*x*prime[k];
}
}
void calc(int op)
{
//for (int i=1;i<=t;i++) cout<<prime[i]<<' '<<cnt[0][i]<<' '<<cnt[1][i]<<endl;cout<<endl;
sum=;get(,);
for (int i=;i<=tot;i++)
{
bool flag=;
for (int j=;j<=t;j++)
if (p[i][j]<cnt[][j]||p[i][j]>cnt[][j]) {flag=;break;}
if (flag)
{
ans[i]+=op*ksm(,sum-);
if (ans[i]<) ans[i]+=P;if (ans[i]>=P) ans[i]-=P;
}
}
}
void dfs(int k,int op)
{
if (k>t) {calc(op);return;}
dfs(k+,op);
cnt[][k]++;dfs(k+,-op);
cnt[][k]--;dfs(k+,op);
cnt[][k]--;dfs(k+,-op);
cnt[][k]++;
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("bzoj5019.in","r",stdin);
freopen("bzoj5019.out","w",stdout);
const char LL[]="%I64d\n";
#else
const char LL[]="%lld\n";
#endif
n=read(),g=read(),l=read(),m=read();
if (l%g) {for (int i=;i<=m;i++) printf("0\n");return ;}
u=l;
for (int i=;i*i<=u;i++)
if (u%i==)
{
prime[++t]=i,cnt[][t]=;u/=i;
while (u%i==) cnt[][t]++,u/=i;
}
if (u>) prime[++t]=u,cnt[][t]=;
u=g;
for (int i=;i<=t;i++)
while (u%prime[i]==) u/=prime[i],cnt[][i]++;
build(,);
for (int i=;i<=tot;i++) d[i]=p[i][];
sort(d+,d+tot+);
for (int i=;i<=tot;i++)
for (int j=;j<=t;j++)
q[i][j]=p[i][j];
for (int i=;i<=tot;i++)
{
int x=lower_bound(d+,d+tot+,q[i][])-d;
for (int j=;j<=t;j++) p[x][j]=q[i][j];
}
dfs(,);
while (m--)
{
int x=read(),y=lower_bound(d+,d+tot+,x)-d;
if (d[y]!=x) {printf("0\n");continue;}
else printf("%d\n",ans[y]);
}
return ;
}

BZOJ5019 SNOI2017遗失的答案(容斥原理)的更多相关文章

  1. BZOJ5019[Snoi2017]遗失的答案——FWT+状压DP

    题目描述 小皮球在计算出答案之后,买了一堆皮肤,他心里很开心,但是一不小心,就忘记自己买了哪些皮肤了.==|||万 幸的是,他还记得他把所有皮肤按照1-N来编号,他买来的那些皮肤的编号(他至少买了一款 ...

  2. bzoj5019: [Snoi2017]遗失的答案

    Description 小皮球在计算出答案之后,买了一堆皮肤,他心里很开心,但是一不小心,就忘记自己买了哪些皮肤了.==|||万 幸的是,他还记得他把所有皮肤按照1-N来编号,他买来的那些皮肤的编号( ...

  3. 【BZOJ5019】[SNOI2017]遗失的答案(FWT,动态规划)

    [BZOJ5019][SNOI2017]遗失的答案(FWT,动态规划) 题面 BZOJ 题解 发现\(10^8\)最多分解为不超过\(8\)个本质不同质数的乘积. 而\(gcd\)和\(lcm\)分别 ...

  4. bzoj 5019: [Snoi2017]遗失的答案【dp+FWT】

    满足GL的组合一定包含GL每个质因数最大次幂个最小次幂,并且能做限制这些数不会超过600个 然后质因数最多8个,所以可以状压f[s1][s2]为选s1集合满足最大限制选s2集合满足最小限制 dfs一下 ...

  5. bzoj 5019 [Snoi2017]遗失的答案

    题面 https://www.lydsy.com/JudgeOnline/problem.php?id=5019 题解 如果L不是G的倍数 答案为0 下面考虑G|L的情况 将G,L质因数分解 设$L= ...

  6. LOJ2257 SNOI2017 遗失的答案 容斥、高维前缀和

    传送门 数字最小公倍数为\(L\)的充分条件是所有数都是\(L\)的约数,而\(10^8\)内最多约数的数的约数也只有\(768\)个.所以我们先暴力找到所有满足是\(L\)的约数.\(G\)的倍数的 ...

  7. luogu P5366 [SNOI2017]遗失的答案

    luogu 首先gcd为\(G\),lcm为\(L\),有可能出现的数(指同时是\(G\)的因数以及是\(L\)的倍数)可以发现只有几百个.如果选出的数要能取到gcd,那么对于每种质因子,都要有一个数 ...

  8. 洛谷$P5366\ [SNOI2017]$遗失的答案 数论+$dp$

    正解:数论$dp$ 解题报告: 传送门$QwQ$ 考虑先质因数分解.所以$G$就相当于所有系数取$min$,$L$就相当于所有系数取$max$ 这时候考虑,因为数据范围是$1e8$,$1e8$内最多有 ...

  9. [SNOI2017]遗失的答案

    题目 首先\(G,L\)肯定会满足\(G|L\),否则直接全部输出\(0\) 之后我们考虑一下能用到的质因数最多只有\(8\)个 同时我们能选择的数\(x\)肯定是\(L\)的约数,还得是\(G\)的 ...

随机推荐

  1. day59

    轮播图作业 <!DOCTYPE html> <html lang="zh"> <head> <meta charset="UTF ...

  2. [Usaco2005 Open]Disease Manangement 疾病管理 BZOJ1688

    分析: 这个题的状压DP还是比较裸的,考虑将疾病状压,得到DP方程:F[S]为疾病状态为S时的最多奶牛数量,F[S]=max{f[s]+1}; 记得预处理出每个状态下疾病数是多少... 附上代码: # ...

  3. net 表格控件

    一个开源的表格控件,界面像Excel,看来很好,有机会在项目中使用:ReoGrid https://reogrid.net/

  4. python基础1之python介绍、安装、变量和字符编码、数据类型、输入输出、数据运算、循环

    开启python之路 内容概要: 一.python介绍 二.安装 三.第一个python程序 四.变量和字符编码 五.用户输入 六.数据类型 七.一切皆对象 八.数据运算 九.if else 流程判断 ...

  5. 双面间谍(spy)

    双面间谍 链接 分析: 戳这 代码: #include<cstdio> #include<algorithm> #include<cstdio> #include& ...

  6. Gitlab+Jenkins学习目录

    Gitlab+Jenkins基础篇 Gitlab+Jenkins学习之路(一)之Git基础 Gitlab+Jenkins学习之路(二)之gitlab部署 Gitlab+Jenkins学习之路(三)之g ...

  7. Java 中的 try catch 影响性能吗?

    前几天在 code review 时发现有一段代码中存在滥用try catch的现象.其实这种行为我们也许都经历过,刚参加工作想尽量避免出现崩溃问题,因此在很多地方都想着 try catch一下. 但 ...

  8. REST-framework快速构建API--四部曲

    代码目录结构: 一.使用原生APIView 使用rest-framework原生的APIView实现过程: 以url(r'^books/$', views.BookView.as_view(),nam ...

  9. docker之compose 编排项目

    一.docker-compose 的介绍 docker-compose是一种容器编排工具,可以将多个docker容器关联部署.通过yaml文件,可以描述应用的架构,如使用什么镜像.数据卷.网络.绑定服 ...

  10. numpy 初识(三)

    基本运算 exp: e sqrt:开放 floor:向下取整 ravel:矩阵拉成一个向 T:转置(行和列变换) 改变形状: resize: 更改其形状(返回值为None)a.resize(6,2) ...