Description

Input

第一行,两个正整数 S 和 q,q 表示询问数量。
接下来 q 行,每行一个正整数 n。

Output

输出共 q 行,分别为每个询问的答案。

Sample Input

30 3
9
29
1000000000000000000

Sample Output

0
9
450000036

HINT

感谢the Loser协助更正数据

对于100%的数据,2<=S<=2*10^6,1<=n<=10^18,1<=q<=10^5

这个题面让人很费解啊

其实题意还是挺简洁的

首先对S分解质因数

如果有相同的因数出现了多次,那么$lcm(p_1,p_2,p_3...p_k)=S$一定不能被满足

此时直接特判输出Q个0即可

(这10分这么好拿还想什么正解hhh)

接下来要做的,是取一定数量的每个$p_i$,使他们的和为$n$

假设$n=\sum{p_i*c_i}$,即$p_i$对应取$c_i$个

由$p_i$是S的因数,可得$p_i*c_i$一定可以表示成$X*S+Y*p_i$

将$p_i$除过去,得到$c_i=X*\frac{S}{p_i}+Y$

那么就可以把$c_i$分成+号前后两部分,不同的$X$和不同的$Y$都会导致方案不同

设前半部分为$a_i$,后半部分为$b_i$

将$c_i\ mod\ \frac{S}{p_i}$得到后半部分,这部分的方案数用多重背包解决

则背包大小不超过$k*S$(每个物品不超过$\frac{S}{p_i}$),

之后我们考虑$a_i$部分,分配每一个$X$等价于将$\frac{S}{p_i}$个物品放到k个盒子内(可空),

利用隔板法得到方案数为$C_{\frac{S}{p_i}+k-1}^{k-1}$

直接打表求k范围内的阶乘逆元,暴力求排列再$*(m-1)$的阶乘逆元即可

$a_i\ b_i$的贡献是相对独立的,所以可以直接相乘

但这样可能包含了几倍的S,

而通过枚举一个倍数m( $p_1*b_1+p_2*b_2+...+p_k*b_k=n-m*S$ )即可得到结果

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<vector>
using namespace std;
typedef long long ll;
const int N=;
const ll mod=1e9+;
int S,Q,tot=,sumf;
int fact[];
int sz,f[],g[],mu[],pr[];
bool v[];
ll n,side,inv[],fac[];
void dp()
{
f[]=;
for(int i=;i<=sz;i++)
{
for(int j=;j<fact[i];j++)
{
ll tmp=;
for(int k=j;k<=side;k+=fact[i])
{
tmp=(f[k]+tmp)%mod;
if(k>=S)tmp=(tmp-f[k-S])%mod;
g[k]=tmp;
}
}
for(ll j=;j<=side;j++)f[j]=g[j];
}
}
ll qpow(ll a,ll b)
{
ll res=;
for( ;b;b>>=,a=a*a%mod)
if(b&)res=res*a%mod;
return res;
}
ll C(ll x,ll y)
{
ll res=;
for(ll i=x;i>=x-y+;i--)res=res%mod*(i%mod)%mod;
res=res%mod*(inv[y]%mod)%mod;
return res;
}
void prime()
{
mu[]=;int num=;
for(int i=;i<=S;i++)
{
if(!v[i])pr[++num]=i,mu[i]=-;
for(int j=;j<=num;j++)
{
if(i*pr[j]>S)break;
v[i*pr[j]]=;
if(i%pr[j])mu[i*pr[j]]=-mu[i];
else
{
mu[i*pr[j]]=;
break;
}
}
}
for(int i=;i<=num;i++)
if(S%pr[i]==)fact[++sz]=pr[i],sumf+=pr[i];
}
int main()
{
scanf("%d%d",&S,&Q);
prime();
if(!mu[S])
{
while(Q--)puts("");
return ;
}
side=S*sz;
inv[]=fac[]=;
for(ll i=;i<=sz;i++)
fac[i]=fac[i-]*i%mod;
inv[sz]=qpow(fac[sz],mod-);
for(int i=sz-;i;i--)
inv[i]=inv[i+]*(i+)%mod;
dp();
while(Q--)
{
scanf("%lld",&n);
if(n<sumf)
{
puts("");
continue ;
}
n-=sumf;
ll tm=n/S,ans=;
for(ll i=;i<=sz;i++)
{
if(n-i*S<)break;
ll now=i*S+n%S;
if(now>side)continue;
ans=(ans+(ll)f[now]*C(tm-i+sz-,sz-))%mod;
}
printf("%lld\n",(ans+mod)%mod);
}
return ;
}

[bzoj3462]DZY Loves Math II (美妙数学+背包dp)的更多相关文章

  1. BZOJ3462 DZY Loves Math II 【多重背包 + 组合数】

    题目 输入格式 第一行,两个正整数 S 和 q,q 表示询问数量. 接下来 q 行,每行一个正整数 n. 输出格式 输出共 q 行,分别为每个询问的答案. 输入样例 30 3 9 29 1000000 ...

  2. BZOJ3462 DZY Loves Math II(动态规划+组合数学)

    容易发现这是一个有各种玄妙性质的完全背包计数. 对于每个质数,将其选取个数写成ax+b的形式,其中x=S/pi,0<b<x.那么可以枚举b的部分提供了多少贡献,多重背包计算,a的部分直接组 ...

  3. bzoj3462: DZY Loves Math II

    状态很差脑子不清醒了,柿子一直在推错.... ... 不难发现这个题实际上是一个完全背包 问题在于n太大了,相应的有质数的数量不会超过7个 假设要求sigema(1~plen)i pi*ci=n 的方 ...

  4. bzoj 3462: DZY Loves Math II

    3462: DZY Loves Math II Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 211  Solved: 103[Submit][Sta ...

  5. DZY Loves Math II:多重背包dp+组合数

    Description Input 第一行,两个正整数 S 和 q,q 表示询问数量.接下来 q 行,每行一个正整数 n. Output 输出共 q 行,分别为每个询问的答案. Sample Inpu ...

  6. BZOJ 3462 DZY Loves Math II ——动态规划 组合数

    好题. 首先发现$p$是互质的数. 然后我们要求$\sum_{i=1}^{k} pi*xi=n$的方案数. 然后由于$p$不相同,可以而$S$比较小,都是$S$的质因数 可以考虑围绕$S$进行动态规划 ...

  7. DZY Loves Math II

    简要题面 对于正整数 \(S, n\),求满足如下条件的素数数列 \((p_1,p_2,\cdots,p_k)\)(\(k\) 为任意正整数) 的个数: \(p_1\le p_2\le\cdots\l ...

  8. DZY Loves Math系列

    link 好久没写数学题了,再这样下去吃枣药丸啊. 找一套应该还比较有意思的数学题来做. [bzoj3309]DZY Loves Math 简单推一下. \[\sum_{i=1}^n\sum_{j=1 ...

  9. [BZOJ] DZY Loves Math 系列 I && II

    为了让自己看起来有点事干 ,做个套题吧..不然老是东翻翻西翻翻也不知道在干嘛... \(\bf 3309: DZY \ Loves \ Math\) 令 \(h=f*\mu\) 很明显题目要求的就是\ ...

随机推荐

  1. Vagrant 入门 - 项目设置

    原文地址 配置 Vagrant 项目的第一步是创建 Vagrantfile 文件.Vagrantfile 文件的目的有两个: 设置项目的根目录.Vagrant 中的许多配置选项是相对于这个根目录的. ...

  2. python 按二维数组的某行或列排序 (numpy lexsort)

    lexsort支持对数组按指定行或列的顺序排序:是间接排序,lexsort不修改原数组,返回索引. (对应lexsort 一维数组的是argsort a.argsort()这么使用就可以:argsor ...

  3. [LeetCode] 342. Power of Four(位操作)

    传送门 Description Given an integer (signed 32 bits), write a function to check whether it is a power o ...

  4. 如何通过xmake进行交叉编译

    xmake 提供了方便灵活的交叉编译支持,大部分情况下,都不需要配置很复杂的toolchains前缀,例如:arm-linux- 什么的 只要这个toolchains目录满足如下结构(大部分的交叉工具 ...

  5. JNDI配置笔记

    先在tomcat Context.xml配置文件中配置 <Resource name="jdbc/elifecrm" type="javax.sql.DataSou ...

  6. 读取交货单拣配数量PIKMG(转)

    原文链接:https://www.591sap.com/thread-953-1-1.html SAP交货单交货数量在lips中直接读取,但是拣配数量lfimg,只存在vbfa中,且如果基本计量单位和 ...

  7. 关于Object.create()与原型链的面试题?

    原文地址 https://segmentfault.com/q/1010000004670616?utm_source=weekly&utm_medium=email&utm_camp ...

  8. Codeforces 1156D 0-1-Tree ( 并查集 || 树形DP )

    <题目链接> 题目大意: 给定一颗无向树,树的边权只要0/1两种情况,现在问你这棵树上存在多少对有序对<u,v>,满足u-->v的路径上,如果出现边权为1的边之后,就不能 ...

  9. jQuery easyUI 使用 datagrid 表格

    获取后台数据依旧是使用一般处理程序(ashx) ,分页上添加一个函数(pagerFilter(data)) 前端代码: <%@ Page Language="C#" Auto ...

  10. 任务21 :了解ASP.NET Core 依赖注入,看这篇就够了

    DI在.NET Core里面被提到了一个非常重要的位置, 这篇文章主要再给大家普及一下关于依赖注入的概念,身边有工作六七年的同事还个东西搞不清楚.另外再介绍一下.NET  Core的DI实现以及对实例 ...