题目

首先\(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. jQuery的两把利器

    1 jQuery核心函数 * 简称: jQuery函数($/jQuery) * jQuery库向外直接暴露的就是$/jQuery * 引入jQuery库后, 直接使用$即可 * 当函数用: $(xxx ...

  2. paper 16 : 计算机视觉领域博客资源

    这是收录的图像视觉领域的博客资源的第一部分,包含:中国内地.香港.台湾 这些名人大家一般都熟悉,本文仅收录了包含较多资料的个人博客,并且有不少更新,还有些名人由于分享的paper.code或者数据集不 ...

  3. lazyload懒加载和swiper轮播懒加载的用法

    对于有较多的图片的网页,使用图片延迟加载,能有效的提高页面加载速度,比如商城网页. lazyload使用方法: 载入 JavaScript 文件: <script src="jquer ...

  4. 去除springboot内置tomcat

    /** * @author zx * @title: ServletInitializer * @projectName activiti * @description: 解决内置tomcat * @ ...

  5. 知识整理:字符串hash

    字符串hash唯一用途是快速判断两字符串是否相等,但存在极小概率假阳性(本来不相等,但算法返回相等). 根本思想是把一个字符串转换为一个整数,要求相同的字符串,对应的这个整数相同,不同的字符串,对应的 ...

  6. 从yjz那里偷来的fread读入挂

    struct fastio{ char s[100005]; int it,len; fastio(){it=len=0;} inline char get(){ if(it<len)retur ...

  7. SPSS如何调用已建立的数据文件

    SPSS如何调用已建立的数据文件 调用已建立的数据文件 SPSS可以调用SPSS(*.sav),Excel(*.xls),dBASE(*.dbf),ASCII(*.dat,*.txt)等数据文件. 2 ...

  8. 新手git遇到的问题

    1. 如何撤销git add,不小心执行了git add . 操作,但是又不能提交所有的文件,因为对应不同的分支,现在怎么样可以将git add 撤销回来? 解决:git reset HEAD(暂时尝 ...

  9. codeforces round#524 D - Olya and magical square /// 大概算是数学规律题?

    题目大意: t 个测试用例  (1≤t≤103) 给定n k  (1≤n≤10^9,1≤k≤10^18) 表示有一个边长为2^n的正方形格子 每次操作只能将一个格子切割为左上左下右上右下的四等分格子 ...

  10. oracle union 和 union all

    java.sql.SQLSyntaxErrorException: ORA-01789: 查询块具有不正确的结果列数 原因: 发现是sql语句用union时的 两个语句查询的字段不一致 解决:将 2个 ...