题目

首先\(G,L\)肯定会满足\(G|L\),否则直接全部输出\(0\)

之后我们考虑一下能用到的质因数最多只有\(8\)个

同时我们能选择的数\(x\)肯定是\(L\)的约数,还得是\(G\)的倍数,还不能超过\(N\),感性理解一下这样的\(x\)显然不多,我们直接\(dfs\)出来

对于每一个数我们把它压成一个\(2\times 8\)的二进制数,第\(2i-1,2i\)位分别表示第\(i\)个质因子是否到上界/下界,我们直接一个\(dp[i][j]\)表示前\(i\)个数选择状态为\(j\)的方案数,同时我们还能处理出一个后缀的\(dp\)值来

对于一个必须选择的\(x\),我们用前缀和后缀做一下或\(FWT\),求所有和这个数的状态或起来为全集的位置的和就好了

复杂度看起来不是很科学,但是常数比较小,就跑过去了

代码

#include <bits/stdc++.h>
#define re register
const int mod = 1e9 + 7;
inline int read() {
char c = getchar();
int x = 0;
while (c < '0' || c > '9') c = getchar();
while (c >= '0' && c <= '9') x = (x << 3) + (x << 1) + c - 48, c = getchar();
return x;
}
std::map<int, int> ma;
int N, G, L, Q, M, n, len;
int f[50005], p[20005];
int cnt[2], b[10], a[10][2], l[10], r[10], beg[10];
int c[850], st[850], ans[850], vis[850];
int dp[850][65537], tp[850][65537], d[65537];
void dfs(int x, int now, int S) {
if (x == cnt[0] + 1) {
c[++n] = now, st[n] = S;
return;
}
int t = beg[x];
for (re int i = l[x]; i <= r[x]; ++i)
if (now <= N / t) {
int s = S;
if (i == l[x])
s |= (1 << (2 * x - 2));
if (i == r[x])
s |= (1 << (2 * x - 1));
dfs(x + 1, now * t, s);
t = t * a[x][0];
} else
break;
}
inline int ksm(int a, int b) {
int S = 1;
for (; b; b >>= 1, a = a * a)
if (b & 1)
S = S * a;
return S;
}
inline int qm(int x) { return x < 0 ? x + mod : x % mod; }
inline void Fwt(int *f, int o) {
for (re int i = 2; i <= len; i <<= 1)
for (re int ln = i >> 1, l = 0; l < len; l += i)
for (re int x = l; x < l + ln; ++x) f[x + ln] = qm(f[x + ln] + o * f[x]);
}
int main() {
N = read(), G = read(), L = read(), Q = read();
M = std::ceil(std::sqrt(L));
if (L % G) {
for (re int i = 1; i <= Q; i++) puts("0");
return 0;
}
for (re int i = 2; i <= M; i++) {
if (!f[i])
p[++p[0]] = i;
for (re int j = 1; j <= p[0] && p[j] * i <= M; j++) {
f[p[j] * i] = 1;
if (i % p[j] == 0)
break;
}
}
std::swap(G, L);
for (re int i = 1; i <= p[0]; i++) {
if (G % p[i])
continue;
a[++cnt[0]][0] = p[i];
while (G % p[i] == 0) G /= p[i], ++r[cnt[0]];
}
if (G > 1)
a[++cnt[0]][0] = G, r[cnt[0]] = 1;
for (re int i = 1; i <= p[0]; i++) {
if (L % p[i])
continue;
a[++cnt[1]][1] = p[i];
while (L % p[i] == 0) L /= p[i], ++b[cnt[1]];
}
if (L > 1)
a[++cnt[1]][1] = L, b[cnt[1]] = 1;
for (re int i = 1; i <= cnt[0]; i++)
for (re int j = 1; j <= cnt[1]; j++)
if (a[j][1] == a[i][0]) {
l[i] = b[j];
break;
}
for (re int i = 1; i <= cnt[0]; ++i) beg[i] = ksm(a[i][0], l[i]);
dfs(1, 1, 0);
for (re int i = 1; i <= n; i++) ma[c[i]] = i;
dp[0][0] = 1;
len = 1 << (2 * cnt[0]);
for (re int i = 1; i <= n; i++)
for (re int j = 0; j < len; j++) {
if (!dp[i - 1][j])
continue;
dp[i][j] = (dp[i][j] + dp[i - 1][j]) % mod;
dp[i][j | st[i]] = (dp[i][j | st[i]] + dp[i - 1][j]) % mod;
}
tp[n + 1][0] = 1;
for (re int i = n; i; --i)
for (re int j = 0; j < len; j++) {
if (!tp[i + 1][j])
continue;
tp[i][j] = (tp[i][j] + tp[i + 1][j]) % mod;
tp[i][j | st[i]] = (tp[i][j | st[i]] + tp[i + 1][j]) % mod;
}
while (Q--) {
int x = read(), pos = ma[x];
if (!pos) {
puts("0");
continue;
}
if (vis[pos]) {
printf("%d\n", ans[pos]);
continue;
}
Fwt(dp[pos - 1], 1), Fwt(tp[pos + 1], 1);
for (re int i = 0; i < len; i++) d[i] = 1ll * dp[pos - 1][i] * tp[pos + 1][i] % mod;
Fwt(d, -1);
for (re int i = 0; i < len; i++)
if ((i | st[pos]) + 1 == len)
ans[pos] = (ans[pos] + d[i]) % mod;
vis[pos] = 1;
printf("%d\n", ans[pos]);
}
return 0;
}

[SNOI2017]遗失的答案的更多相关文章

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

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

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

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

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

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

  4. BZOJ5019 SNOI2017遗失的答案(容斥原理)

    显然存在方案的数一定是L的因数,考虑对其因子预处理答案,O(1)回答. 考虑每个质因子,设其在g中有x个,l中有y个,则要求所有选中的数该质因子个数都在[x,y]中,且存在数的质因子个数为x.y.对于 ...

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

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

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

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

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

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

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

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

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

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

随机推荐

  1. springboot连接mysql数据库,JdbcTemplate和spring JPA方式

    SQL部分 CREATE TABLE test( id ) primary key, name ) not null, age ), address ) ); ,,'bj'); ,,'sh'); ,, ...

  2. bzoj1096题解

    [解题思路] 预处理spi=∑pj(j∈[1,i]),si=si-1+(xi-xi-1)*spi-1表示把工厂1~i-1的产品都运到工厂i的花费.于是把工厂j+1~i的产品都运到工厂i的花费为si-s ...

  3. CF601C Kleofáš and the n-thlon(期望+前缀和优化dp)

    传送门 解题思路 要求这个人的排名,我们可以先求出某个人比他排名靠前的概率,然后再乘上\(m-1\)即为答案.求某个人比他排名靠前可以用\(dp\),设\(f[i][j]\)表示前\(i\)场比赛某人 ...

  4. NOIp2018集训test-9-18

    T1.Conjugate 只能选没选过的点,看成如果选了选过的堆的点就不管它继续选.如果第一次选到某一堆的点在第一次选到第一堆的点之前,这一堆对答案就会有1的贡献.那么a[i]有贡献的概率是a[i]和 ...

  5. 机械手相机9点坐标标定-基于C#+EmguCV

    很多初学者,都对标定概念模糊不清,分不清坐标系之间的关系,搞不清相机标定和机械手相机标定有什么关系,想当初自己也是一个人摸索了很久,本文将尽量给大家解释. 我们通常所说的相机标定分为两种,一种是相机参 ...

  6. opencv-角点检测之Harris角点检测

    转自:https://blog.csdn.net/poem_qianmo/article/details/29356187 先看看程序运行截图:   一.引言:关于兴趣点(interest point ...

  7. spark-sql性能优化之——动态实现多个列应用同一个函数

    在对一个dataframe的多个列实现应用同一个函数时,是否能动态的指定? 例如: 对A,B,C三列实现分组统计 1.初始化spark,构建DF val spark = SparkSession.bu ...

  8. oracle 管理表空间

    表空间:是oracle数据库中最大的逻辑存储结构,与操作系统中的数据文件相对应,用于存储数据库中用户创建的所有内容 表空间>数据文件 4.1创建一个永久性表空间myspace create ta ...

  9. this.$router.go()和this.$router.push()的差别

    1.this.$router.go(val) => 在history记录中前进或者后退val步,当val为0时刷新当前页面. 2.this.$router.push(path) => 在h ...

  10. 金三银四铜五铁六,Offer收到手软!

    作者:鲁班大师 来源:cnblogs.com/zhuoqingsen/p/interview.html 文中的鲁班简称LB 据说,金三银四,截止今天为止面试黄金时间已经过去十之八九,而LB恰逢是这批面 ...