BZOJ5302: [Haoi2018]奇怪的背包

https://lydsy.com/JudgeOnline/problem.php?id=5302

分析:

  • 方程\(\sum\limits_{i=1}^nx_ia_i=y\)有整数解的条件是\(gcd|y\)。
  • 对于这道题,我们可以直接把\(P\)当成一个可以提供负值 or 可以抵消负值的存在。
  • 那么这道题的条件就是\(gcd(d,P)|q_i\),其中\(d\)是选出那些数的\(gcd\)。
  • 问题缩小到了\(P\)的约数这个范围,最多\(1400\)个。
  • 然后一通预处理就做完了,询问直接用\(gcd(P,q_i)\)对应答案来回答。

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
using namespace std;
#define N 1000050
#define mod 1000000007
typedef long long ll;
const int S=100000;
char buf[100000],*p1,*p2;
#define nc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++)
int rd() {
int x=0; char s=nc();
while(s<'0') s=nc();
while(s>='0') x=(((x<<2)+x)<<1)+s-'0',s=nc();
return x;
}
int n,m,a[N],P;
int id1[S+10],id2[S+10],w[1450],tot,b[N];
ll f[2][1450],ans[1450];
int gcd(int x,int y) {return y?gcd(y,x%y):x;}
inline void upd(ll &x,ll y) {x+=y; if(x>=mod)x-=mod;}
inline int pos(int x) {return x<=S?id1[x]:id2[P/x];}
inline ll qp(ll x,ll y) {
ll re=1;
for(;y;y>>=1,x=x*x%mod) if(y&1) re=re*x%mod; return re;
}
char pbuf[100000],*pp=pbuf;
void push(const char c) {
if(pp-pbuf==100000) fwrite(pbuf,1,100000,stdout),pp=pbuf;
*pp++=c;
}
int sta[20],tp;
void write(int x) {
do {
sta[++tp]=x%10;
x/=10;
}while(x);
while(tp) push(sta[tp--]+'0');
push('\n');
}
int main() {
n=rd(),m=rd(),P=rd();
int i,j;
for(i=1;i<=n;i++) a[i]=rd();
for(i=1;i*i<=P;i++) {
if(P%i==0) {
w[++tot]=i;
if(i*i!=P) w[++tot]=P/i;
}
}
sort(w+1,w+tot+1);
for(i=1;i<=tot;i++) {
if(w[i]<=S) id1[w[i]]=i;
else id2[P/w[i]]=i;
}
for(i=1;i<=n;i++) b[pos(gcd(a[i],P))]++;
for(i=1;i<=tot;i++) b[i]=qp(2,b[i])-1;
f[0][0]=1;
for(i=0;i<tot;i++) {
int i0=i&1,i1=(i+1)&1;
for(j=0;j<=i+1;j++) if(f[i0][j]) {
f[i1][j]=(f[i1][j]+f[i0][j])%mod;
int k=pos(gcd(w[j],w[i+1]));
f[i1][k]=(f[i1][k]+f[i0][j]*b[i+1])%mod;
}
memset(f[i0],0,sizeof(f[i0]));
}
int l=tot&1;
for(i=1;i<=tot;i++) {
for(j=1;j<=i;j++) if(w[i]%w[j]==0) {
ans[i]=(ans[i]+f[l][j]);
}
ans[i]%=mod;
}
while(m--) {
int x; x=rd();
write(ans[pos(gcd(x,P))]);
}
fwrite(pbuf,1,pp-pbuf,stdout);
}

BZOJ5302: [Haoi2018]奇怪的背包的更多相关文章

  1. BZOJ5302 [HAOI2018]奇怪的背包 【数论 + dp】

    题目 小 CC 非常擅长背包问题,他有一个奇怪的背包,这个背包有一个参数 PP ,当他 向这个背包内放入若干个物品后,背包的重量是物品总体积对 PP 取模后的结果. 现在小 CC 有 nn 种体积不同 ...

  2. BZOJ5302 HAOI2018奇怪的背包(动态规划)

    由裴蜀定理,子集S有解当且仅当gcd(S,P)|w. 一个显然的dp是设f[i][j]为前i个数gcd为j的选取方案.注意到这里的gcd一定是P的约数,所以状态数是n√P的.然后可以通过这个得到gcd ...

  3. [BZOJ5302][HAOI2018]奇怪的背包(DP)

    由裴蜀定理得,一个集合S能得到w当且仅当gcd(S+{P})|w. 于是f[i][j]表示前i个物品gcd为j的方案数,发现gcd一定是P的因数,故总复杂度$O(n\sqrt{P}\log P)$(需 ...

  4. 【BZOJ5302】[HAOI2018]奇怪的背包(动态规划,容斥原理)

    [BZOJ5302][HAOI2018]奇怪的背包(动态规划,容斥原理) 题面 BZOJ 洛谷 题解 为啥泥萌做法和我都不一样啊 一个重量为\(V_i\)的物品,可以放出所有\(gcd(V_i,P)\ ...

  5. [HAOI2018]奇怪的背包 (DP,数论)

    [HAOI2018]奇怪的背包 \(solution:\) 首先,这一道题目的描述很像完全背包,但它所说的背包总重量是在模P意义下的,所以肯定会用到数论.我们先分析一下,每一个物品可以放无数次,可以达 ...

  6. 洛谷 P4495 [HAOI2018]奇怪的背包 解题报告

    P4495 [HAOI2018]奇怪的背包 题目描述 小\(C\)非常擅长背包问题,他有一个奇怪的背包,这个背包有一个参数\(P\),当他 向这个背包内放入若干个物品后,背包的重量是物品总体积对\(P ...

  7. haoi2018奇怪的背包题解

    题目传送门:https://www.lydsy.com/JudgeOnline/problem.php?id=5302 对于一个物品,设它体积为v,那么,在背包参数为p的情况下,它能达到gcd(v,p ...

  8. bzoj 5302: [Haoi2018]奇怪的背包

    Description Solution 首先 \(v_1,v_2,v_3...v_n,P\) 能够构成的最小数是 \(gcd(P,v_1,v_2,v_3...v_n)\) 然后 \(gcd(P,v_ ...

  9. Luogu4495 [HAOI2018] 奇怪的背包 【扩展欧几里得算法】

    题目分析: 首先打个暴力求一下$10^9$以内因子最多的数的因子个数,发现只有$1344$个. 由于有$ax+by=k*(a,b)$和2017年noip的结论,所以我们可以发现对于任意多个数$a_1, ...

随机推荐

  1. spark学习(1)--ubuntu14.04集群搭建、配置(jdk)

    环境:ubuntu14.04 jdk-8u161-linux-x64.tar.gz 1.文本模式桌面模式切换 ctrl+alt+F6 切换到文本模式 ctrl + alt +F7 /输入命令start ...

  2. 31 整数中1出现的次数(从1到n整数中1出现的次数)

    题目描述 求出1~13的整数中1出现的次数,并算出100~1300的整数中1出现的次数?为此他特别数了一下1~13中包含1的数字有1.10.11.12.13因此共出现6次,但是对于后面问题他就没辙了. ...

  3. 026_默认的MapReduce Driver(最小驱动问题)

    1. 最小配置的MapReduce Driver 读取输入文件中的内容,输出到指定目录的输出文件中,此时文件中的内容为: Key---输入文件每行内容的起始位置. Value---输入文件每行的原始内 ...

  4. 基于SSM的单点登陆01

    使用SSM的Maven聚合项目 建立父项目market的pom文件 <?xml version="1.0" encoding="UTF-8"?> & ...

  5. jsonp跨域总结

    同源限制: 浏览器不允许一个域的脚本请求另一个域的文档,通俗来说,不允许访问不同协议.不同域名或不同端口的文档 跨域处理方法: 1.jsonp 前提: 浏览器虽然对ajax请求做出了限制,但scrip ...

  6. Go reflect反射

    Go语言中的反射非常强大,可以对string, int, struct, func...进行反射,使用起来也比较简单. 示例1:反射函数 package main import ( "fmt ...

  7. 【转载】poj 1276 Cash Machine 【凑钱数的问题】【枚举思路 或者 多重背包解决】

    转载地址:http://m.blog.csdn.net/blog/u010489766/9229011 题目链接:http://poj.org/problem?id=1276 题意:机器里面共有n种面 ...

  8. WEB开发中常见漏洞

    1.sql注入 SQL注入在黑客领域是一种非常常见的攻击手段,大家应该都听说过很多数据泄漏的案例,其中大部分都是采用SQL注入来获取数据的. SQL注入一般是前端向后台提交数据的时候,在数据中加入SQ ...

  9. java基础(8)-集合类

    增强for循环 /* *增强for循环 * for(元素类型 变量:数据或Collection集合){ * 使用变量即可,该变量就是元素 * } * 优点:简化了数组和集合的遍历 * 缺点:增强for ...

  10. MapReduce-边数据

    边数据 边数据(side data)是作业所需的额外的只读数据,以辅助处理主数据集.所面临的挑战在于如何使所有map或reduce任务(这些任务散布在集群内部)都能够方便而高效地使用边数据. 利用Jo ...