题目

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

i

​ ,每种物品都有无限个. 他会进行 qq 次询问,每次询问给出重量 w_iw

i

​ ,你需要回答有多少种放入物品的方 案,能将一个初始为空的背包的重量变为 w_iw

i

​ .注意,两种方案被认为是不同的, 当且仅当放入物品的种类不同,而与每种物品放入的个数无关.不难发现总的方 案数为 2^n2

n

. 由于答案可能很大,你只需要输出答案对1e9+7取模的结果.

输入格式

从文件 knapsack.inknapsack.in 中读入数据. 第一行三个整数 nn , qq , PP ,含义见问题描述. 接下来一行 nn 个整数表示 V_iV

i

​ . 接下来一行 qq 个整数表示 w_iw

i

​ .

输出格式

输出到文件 knapsack.outknapsack.out 中. 输出 qq 行,每行一个整数表示答案.

输入样例

3 3 6

1 3 4

5 2 3

输出样例

5

6

6

提示

题解

考虑数论中的\(ax + by\)这样线性组合的形式,其只能表示\(gcd(a,b)\)的倍数

\(ax \mod P\)可以写成\(ax - Py\),也是\(a\)与\(P\)的线性组合

所以一个物品只能表示\(gcd(V[i],P)\)的倍数

我们只需要求出所有\(P\)的约数的方案数,对于询问\(W[i]\),我们只需要查询\(gcd(W[i],P)\)的倍数的方案数

设\(f[i][j]\)表示前\(i\)个数凑出以约数\(j\)为线性组合中的最小值的方案数【即\(j\)对其倍数的贡献】,转移即可

\(10^9\)以内约数最大不超过\(10^3\)个

建议别用STL

#include<iostream>
#include<cstdio>
#include<cmath>
#include<map>
#include<cstring>
#include<algorithm>
#define LL long long int
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define BUG(s,n) for (int i = 1; i <= (n); i++) cout<<s[i]<<' '; puts("");
#define MAP map<int,int>
#define res register
using namespace std;
const int maxn = 1000005,maxm = 100005,INF = 1000000000,md = 1e9 + 7;
inline int read(){
res int out = 0,flag = 1; res char c = getchar();
while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();}
while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}
return out * flag;
}
inline void write(int x){if (x >= 10) write(x / 10); putchar(x % 10 + '0');}
MAP f[2],Ans,A;
MAP::iterator it,IT;
int n,m,P;
inline int add2(int x,int y){x += y; if (x >= md) x -= md; return x;}
inline void add(int& x,int y){x += y; if (x >= md) x -= md;}
inline int gcd(int a,int b){return b ? gcd(b,a % b) : a;}
void init(){
int tmp,p = 1; f[0][P] = 1;
for (IT = A.begin(); IT != A.end(); IT++,p ^= 1){
f[p].clear();
for (it = f[p ^ 1].begin(); it != f[p ^ 1].end(); it++){
add(f[p][it->first],it->second);
tmp = gcd(IT->first,it->first);
add(f[p][tmp],1ll * add2(IT->second,md - 1) * it->second % md);
}
}
p ^= 1;
for (res int i = 1; 1ll * i * i <= P; i++) Ans[i] = Ans[P / i] = 0;
for (it = Ans.begin(); it != Ans.end(); it++)
for (IT = f[p].begin(); IT != f[p].end(); IT++)
if (it->first % IT->first == 0)
add(it->second,IT->second);
}
void solve(){
while (m--) write(Ans[gcd(read(),P)]),puts("");
}
int main(){
n = read(); m = read(); P = read();
int x;
for (res int i = 1; i <= n; i++){
x = gcd(read(),P);
if (!A.count(x)) A[x] = 2;
else A[x] = A[x] * 2 % md;
}
init();
solve();
return 0;
}

BZOJ5302 [HAOI2018]奇怪的背包 【数论 + dp】的更多相关文章

  1. BZOJ5302: [Haoi2018]奇怪的背包

    BZOJ5302: [Haoi2018]奇怪的背包 https://lydsy.com/JudgeOnline/problem.php?id=5302 分析: 方程\(\sum\limits_{i=1 ...

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

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

  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奇怪的背包(动态规划)

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

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

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

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

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

  7. 洛谷P4495 [HAOI2018]奇怪的背包(数论)

    题面 传送门 题解 好神仙的思路啊--orzyyb 因为不限次数,所以一个体积为\(V_i\)的物品可以表示出所有重量为\(\gcd(V_i,P)\)的倍数的物品,而所有物品的总和就是这些所有的\(\ ...

  8. haoi2018奇怪的背包题解

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

  9. [HAOI2018]奇怪的背包

    题目 暴力\(dp\)好有道理啊 于是我们来个反演吧 考虑一个体积序列\(\{v_1,v_2,...v_n\}\)能凑成\(w\)的条件 显然是 \[v_1x_1+v_2x_2+...+v_nx_n\ ...

随机推荐

  1. python_71_json序列化1

    #序列化:序列化 (Serialization)将对象的状态信息转换为可以存储或传输的形式的过程. #本例把字典数据类型存成字符串存在硬盘 #文件只能存字符串和二进制码,字典之类的不可以 info={ ...

  2. C#继承机制 多级继承

    一些面向对象语言允许一个类从多个基类中继承,而另一些面向对象语言只允许从一个类继承,但可以随意从几个接口或纯抽象类中继承. 只有C++支持多级继承,许多程序员对此褒贬不一.多级继承常会引起继承来的类之 ...

  3. 图片url转base64

    var xhr = new XMLHttpRequest() // 配置的代理,解决跨域问题 xhr.open('GET', url.replace('http://xxx.com', '/img') ...

  4. vue学习之路 - 2.基本操作(上)

    基本操作(上) 本章节简介: vue的安装 vue实例创建 数据绑定渲染 表单数据双向绑定 事件处理 安装 安装方式有三种: 一.vue官网直接下载 http://vuejs.org/js/vue.m ...

  5. vue.js 图表chart.js使用

    在使用这个chart.js之前,自己写过一个饼图,总之碰到的问题不少,所以能用现成的插件就用,能节省不少时间 这里不打算介绍chart.js里面详细的参数意义和各个参数的用法,只作为首次使用chart ...

  6. MongDB之各种修改操作

    接口IMongDaoUpdate: package com.net.test.mongdb.dao; import com.net.test.mongdb.entity.User; public in ...

  7. 如何使用pyinstaller打包32位的exe

    说明:原来安装的python为64位,故安装的pyinstaller和打包后的exe都为64位.而64位的exe文件在32位的win7操作系统下是无法执行的,显示不兼容.网上查询发现,简单(可能不方便 ...

  8. 模块numpy的用法

    https://blog.csdn.net/qq351469076/article/details/78817378 机器学习三剑客之Numpy: https://www.jianshu.com/p/ ...

  9. 安装 ubuntu 后,使用 sed 更换国内源

    cd /etc/aptsed -i "s/archive.ubuntu.com/mirrors.aliyun.com/g" /etc/apt/sources.list也可以使用 1 ...

  10. Disharmony Trees HDU - 3015

    Disharmony Trees HDU - 3015 One day Sophia finds a very big square. There are n trees in the square. ...