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. 值得关注的10个Python语言学习博客

    大家好,还记得我当时学习python的时候,我一直努力地寻找关于python的博客,但我发现它们的数量很少.这也是我建立这个博客的原因,向大家分享我自己学到的新知识.今天我向大家推荐10个值得我们关注 ...

  2. Kattis - sortofsorting 【排序】

    题意 给出一系列字符串,然后要排序 排序规则 只按前两位按字典序来排序,如果前两位完全一样,则按输入的顺序来排 思路 要用 冒泡排序 不能用STL里面的 SORT 因为它不稳定 AC代码 #inclu ...

  3. PHPExcel读写封装

    <? require_once ('inc/PHPExcel/PHPExcel/IOFactory.php'); /** * @author lgl * 使用实例 * $fieldMap=['昵 ...

  4. windows7下手工搭建Apache2.2 php5.3 Mysql5.5开发环境

    Apache2.2(apache_2.2.2-win32-x86-no_ssl)php5.3.5(php-5.3.5-Win32-VC6-x86,请注意选择VC6版本,否则无法加载php5apache ...

  5. iOS_触摸事件与手势识别

    目  录: 一.触摸事件 1.1iOS的输入事件 1.2 触摸事件的处理 1.3 UITouch类中包含五个属性 1.4 UITouch类中包含两个成员函数 1.5响应者链 二.手势识别 2.1使用手 ...

  6. 同类型元素,只有一个被选中js

    <div class="wrap-box flex_row"> <div class="wrap-block"> <div cla ...

  7. Linux基本常用命令

    说到Linux,它就是基于POSIX和UNIX的多用户,多任务,支持多线程和多CPU的操作系统.它能运行主要的UNIX的工具软件,应用程序和网络协议.它支持32位和64位硬件.linux继承Unix以 ...

  8. HTML5 画布canvas

    SVG的<defs> <symbols> 元素用于预定义一个元素使其能够在SVG图像中重复使用 <svg xmlns="http://www.w3.org/20 ...

  9. Windows 7安装PHP运行环境和开发环境

    1. 安装Apache 下载地址:http://www.apache.org/dyn/closer.cgi/httpd/binaries/win32 如需更改端口:打开Apache安装目录下conf目 ...

  10. vs2010下创建webservice ----第一天(建立项目,以及不连数据库进行加减乘除)

    Visual Studio 2010默认采用的框架为.NET Framework4,在这个框架中已找不到直接创建WebService的模板方式了.但VS2010可以创建WebService是毋庸置疑的 ...