最新要闻
- 10万元到1万元 激光电视降价为了什么
- 全球新消息丨ASML:芯片供应链脱钩自主生产不可能 美国日本也不例外
- 《消失的她》票房破5亿:打破44项影史纪录 端午档最强国产悬疑电影诞生 全球最新
- 五年来最火端午假期正式结束!这些城市最热门 你去了吗?
- 鸿海加快iPhone 15在印度量产:与中国大陆时间差将缩减至一个月-天天速递
- 400万播放量 救不了东北偶像
- 标压i7+32GB仅2999元 这款迷你主机能入手吗?实测来了
- 党的最大优势是党执政后的最大危险是_党执政后的最大危险是
- 全球新消息丨韩国yg公司旗下艺人_韩国yg公司
- k360次列车硬卧分布图(k360次列车)|天天时讯
- 安徽唯一盲人考生 高考成绩出来了:311分
- 因脑溢血开颅3次的少年高考445分 母亲:努力有了结果
- garage两种读音_garage怎么读有几种读法是对的
- 观焦点:女子醉酒吐车上 20元都不肯赔 司机怒骂不要脸:网友力挺
- 环球快报:再获11亿美元投资:蔚来“长期主义”的底气
- 距离今年“召回王”又进一步 本田海外召回超百万辆汽车
手机
iphone11大小尺寸是多少?苹果iPhone11和iPhone13的区别是什么?
警方通报辅警执法直播中被撞飞:犯罪嫌疑人已投案
- iphone11大小尺寸是多少?苹果iPhone11和iPhone13的区别是什么?
- 警方通报辅警执法直播中被撞飞:犯罪嫌疑人已投案
- 男子被关545天申国赔:获赔18万多 驳回精神抚慰金
- 3天内26名本土感染者,辽宁确诊人数已超安徽
- 广西柳州一男子因纠纷杀害三人后自首
- 洱海坠机4名机组人员被批准为烈士 数千干部群众悼念
家电
『题解』BZOJ3462 DZY Loves Math II
(相关资料图)
前言
没啥前言,摆了摆了。
题面长这个样子
思路
没啥思路,摆了摆了。这题总的来说挺难想的,思考过程比较繁琐,我也就不辞辛劳列举一下。
显然,条件 \(2\) 和条件 \(3\) 很好说,放一边。
我们设一个质数数列 \(\{p_k\}\) ,假定这些质数是由 \(S\) 分解坤坤数
质因数得到的,若要满足条件 \(4\) (即 \(\mathbb{lcm}\{p_k\} = S\) ) ,则需要满足 \(\{p_k\}\) 中的元素两两互不相等,否则无解。\(\\\) 证明显然。 \(\\\) 设 \(sum = \sum_{i = 1}^{k}{p_i}\) ,我们讨论另一种无解的情况。当 \(n < sum\) 时无解。\(\\\) 证明如下:\(\\\) 因为 \(n\) 一定由 \(\{p_k\}\) 中的全部元素组成,当然,每个元素可以有多个。若 \(n < sum\) ,则 \(n\) 一定不能由 \(\{p_k\}\) 中的全部元素组成(因为在这种情况下 \(sum = \sum_{i = 1}^{k}{p_i} > n\) ,就算每个元素只选一个,这些选的元素相加也会比 \(n\) 大),所以便不能满足条件 \(4\) 。下面我们只考虑有解的情况。乍一看,这个题很像多重背包。\(\\\) 你说得对,但是这题数据范围很大,只用多重背包会炸掉。\(\\\) 所以,我们便要采取一些优化手段。\(\\\) 设 \(S = p_i \times k_i\) , 则 \(k_i\) 的含义便是除了坤坤子 \(p_i\) 外所有坤坤子的乘积。\(\\\) 设 \(n = \sum_{i = 1}^{k}{p_i \times c_i}\) ,\(c_i\) 是坤坤子 \(p_i\) 出现的个数,对于每一项 \(p_i \times c_i\) ,我们可以写成 \(p_i \times (x_i \times k_i + y_i),(y_i < k_i)\) ,拆开括号得 \(p_i \times k_i \times x_i + p_i \times y_i\) ,即 \(S \times x_i + p_i \times y_i , (p_i \times y_i < S)\) 。把每一项加在一起便可得 \(n = S \times \sum_{i = 1}^{k}{x_i} + \sum_{i = 1}^{k}{p_i \times y_i} , (\sum_{i = 1}^{k}{p_i \times y_i} < k \times S)\) 。简化一下(就不设了,知道什么跟什么对应就行)可得 \(n = a \times S + b,(b < k \times S)\) 。\(\\\) 有了这个式子(别管为什么能这么得到这个式子),我们便用加号前面的跑组合数,加号后面的跑
多重背包就行。组合数怎么跑?对于 \(a \times S\) ,因为 \(\{p_k\}\) 中每一个元素都是 \(S\) 的约数,所以每个 \(S\) 都可以用任何一个 \(p_i\) 表示。因为有 \(k\) 个坤坤数,所以 \(S\) 最多能表示成 \(k\) 种形式(有的可以不用表示),一共有 \(a\) 个这样的数,所以转化一下就是插板法求组合数,(根据条件 \(3\) 可知是有序的)把 \(a\) 个数分成 \(k - 1\) 个可空的部分,答案显然是 \(\mathrm{C}_{a + k - 1}^{k - 1}\)。
背包怎么跑?我们不知道 \(a \times S\) 和 \(b\) 中是否都存在 \(p_i\) ,所以我们用 \(2\) 的 \(sum\) ,让 \(b = b - sum\) ,保证每一个 \(p_i\) 都存在,然后我们就可以跑完全背包了。这题让求方案数,
稍微改一下就行了。我们先枚举每一个坤坤子 \(p_i\) ,先都加上,然后再把多的减去即可(这里我犯懒了Orz)。
完结撒代码
#include #define LL long longusing namespace std;const int MOD(1e9 + 7), maxn(2000005);inline LL read() { LL f(1), x(0); char c = getchar(); for (; !isdigit(c); c = getchar()) if (c == "-") f = -1; for (; isdigit(c); c = getchar()) x = (x << 1) + (x << 3) + (c & 15); return f * x;}LL s, q, n, cnt, sum, M, ans;LL inv[maxn], prime[maxn], v[10], dp[7 * maxn];void DP() { dp[0] = 1; for (int i = 1; i <= cnt; ++i) { for (int j = 0; j + v[i] <= M; ++j) { if (dp[j]) dp[j + v[i]] = (dp[j] + dp[j + v[i]] + MOD) % MOD; } for (int j = M - s; j >= 0; --j) { dp[j + s] = (dp[j + s] - dp[j] + MOD) % MOD; } }}bool divide(LL n) { LL x = n; for (int i = 2; i <= sqrt(x); ++i) { if (!(x % i)) { v[++cnt] = i; sum += i; } while (!(x % i)) { ++prime[i]; if (prime[i] == 2) return false; x /= i; } } if (x > 1) { ++prime[x]; v[++cnt] = x; sum += x; } return true;}LL QuickPow(LL a, LL b, LL p) { LL res(1); for (; b; b >>= 1) { if (b & 1) res = res * a % p; a = a * a % p; } return res % p;}int main() { s = read(), q = read(); if (!divide(s)) { while (q--) { n = read(); printf("0\n"); } return 0; } M = cnt * s; DP(); inv[1] = 1; for (int i = 2; i <= 8; ++i) { inv[i] = inv[i - 1] % MOD * QuickPow(i, MOD - 2, MOD) % MOD; } while (q--) { n = read(); if (n < sum) { printf("0\n"); continue; } n -= sum; for (int i = 0; i < cnt && i <= n / s; ++i) { LL a = n / s - i; LL b = n % s + i * s; int res1(1), res2(0); res1 = inv[cnt - 1] % MOD; for (int j = 1; j < cnt; ++j) { res1 = res1 * ((a + j) % MOD) % MOD; } res1 %= MOD; res2 = res1 % MOD * dp[b] % MOD; ans = (ans + res2) % MOD; } ans = (ans % MOD + MOD) % MOD; printf("%lld\n", ans); ans = 0; }}
关键词:
-
读发布!设计与部署稳定的分布式系统(第2版)笔记11_无限长的结果集
![](https: img2023 cnblogs com blog 3076680 202306 3076680-202306
来源: 『题解』BZOJ3462 DZY Loves Math II
读发布!设计与部署稳定的分布式系统(第2版)笔记11_无限长的结果集
10万元到1万元 激光电视降价为了什么
全球新消息丨ASML:芯片供应链脱钩自主生产不可能 美国日本也不例外
全球速讯:k8s 深入篇———— pod 深入实战[七]
CAP定理
【环球速看料】强化学习从基础到进阶-常见问题和面试必知必答[4]::深度Q网络-DQN、double DQN、经验回放、rainbow、分布式DQN
《消失的她》票房破5亿:打破44项影史纪录 端午档最强国产悬疑电影诞生 全球最新
五年来最火端午假期正式结束!这些城市最热门 你去了吗?
鸿海加快iPhone 15在印度量产:与中国大陆时间差将缩减至一个月-天天速递
400万播放量 救不了东北偶像
标压i7+32GB仅2999元 这款迷你主机能入手吗?实测来了
党的最大优势是党执政后的最大危险是_党执政后的最大危险是
全球新消息丨韩国yg公司旗下艺人_韩国yg公司
每日动态!使用libavcodec将mp3音频文件解码为pcm音频采样数据【[mp3float @ 0x561c1ec49940] Header missing
小试Blazor——实现Ant Design Blazor动态表单
焦点快看:【解决方法】按键精灵 实现 狂野飙车9 自动领取每日广告卡牌包 (1)
k360次列车硬卧分布图(k360次列车)|天天时讯
安徽唯一盲人考生 高考成绩出来了:311分
因脑溢血开颅3次的少年高考445分 母亲:努力有了结果
k8s 深入篇———— pod 实战[六]
garage两种读音_garage怎么读有几种读法是对的
观焦点:女子醉酒吐车上 20元都不肯赔 司机怒骂不要脸:网友力挺
环球快报:再获11亿美元投资:蔚来“长期主义”的底气
距离今年“召回王”又进一步 本田海外召回超百万辆汽车
法媒:罗马准备听取对贝洛蒂的报价,意甲英超西甲多队有意
腾讯两大国民APP账号又打通了!QQ悄然支持微信登陆 有了微信为啥仍离不开QQ?
2023年端午小长假收官:1.06亿人次出游 满血超越2019-当前热讯
【世界独家】word安装(ppt、excel),以及word的常用设置
【当前独家】Blazor开发小游戏?趁热打铁上!!!
当前热门:火狐主页diy设置以及常用的扩展
DZY Loves Math|全球即时看
女生发烧考出618分:一蹦三尺高 挨个房间报喜
每日消息!ChatGPT漏洞 讲故事送Window11激活Key!
世界快播:刘慈欣谈ChatGPT:人类的无能反而是人类最后的屏障
今日热议:关于高考志愿填报,这些热点问题需要关注
缓存一致性如何保障
推出长达7年:任天堂股东质疑Switch已逼近极限
实时焦点:苹果前总监炮轰App Store存在灰色地带 标准随心所欲
三乙醇胺油酸皂商品报价动态(2023-06-24)
救命一声吼!山洪暴发女子大喊提醒救下多名游客 世界讯息
索尼PS5串流掌机价格曝光:最高2100元能接受么?
全球视点!泰坦号观光潜艇“打破常规留名后世”,老板一语成谶片受热议
DLang 与 C 语言交互
Apollo2.1.0+Springboot使用OpenApI
邮箱:微信企业域名邮箱给gmail或hotmail等域外邮箱发邮件被退回问题如何解决? 环球观焦点
今日热门!Stable Diffusion模型发布新版本:生成图像以假乱真
钻石价格,突发“跳水”!未来还会更便宜?
使用python发送sip协议的OPTIONS 热门
k8s 深入篇———— k8s 的pod[五]-全球播资讯
8. Java-AOP 面向切面编程
文心一言 VS 讯飞星火 VS chatgpt (46)-- 算法导论6.1 4题|全球热点评
泰坦号事故后:加拿大将展开事故调查
车主自曝差点被闷死在特斯拉Model X里 车门锁死 原因揭晓
蔡徐坤巡演新加坡站开票 《Hug me(remix版)》同日上线
来一打自建IP Proxy玩玩之Majora
kafka学习之五_多个磁盘的性能验证 世界快看点
Go——常用函数
每日速递:卷福的十年同学会
每日简讯:4亿票房端午黑马:《消失的她》官宣海外定档
冲入球场拥抱梅西小伙获释后道歉:我真不是没素质的人
世界观天下!新会绿美生态园票价(新会绿美生态园票价多少钱)
腾讯两大国民APP账号又打通了?QQ悄然支持微信登陆 环球精选
环球快报:调查称安卓更易上手:iPhone用户遇到问题概率高出58%
java 异常处理,事务管理,事务共用,事务传递 天天微头条
Go-闭包和defer|最新资讯
环球热资讯!Zen3清库存?突然冒出个很特别的锐龙5 5700
虚幻5打造!腾讯动漫《斗罗大陆2》今日两集首播 霍雨浩初入星斗大森林|每日观察
“空中出租车”亮相巴黎航展:可降落空间直径仅需15米-环球速递
古力娜扎曾遭换脸视频威胁勒索:不给钱就毁了你!
男生单曲循环《好运来》查出593分大哭:比平时多出50分 超常发挥
全球钻石价格较峰值暴跌18%:人造钻石市场规模不断扩大_全球视点
环球最新:甲亢遇到异食癖:法国男子一顿吃15人份 急了还吃石头木塞
全球热点评!公鸡突然从背后“偷袭”萌娃 飞起两脚踹倒在地 第二天端午节就被炖了
泰国和美国两地大量鱼类死亡 或与海洋升温有关 快看点
【技术积累】C语言中的指针【一】_世界百事通
Go-自定义数据类型(函数类型)详解
行业风险管理需求强烈
焦点信息:AMD RX 7800被逼急了!硬塞进去个“大胖子”
世界视点!一考生查分 全家一起喊出“666”:打算冲击复旦、交大
Kafka学习之四_Grafana监控相关的学习
一天吃透MySQL面试八股文 环球微速讯
什么是大模型? 每日热讯
内马尔在足球界的地位_内马尔的盘带水平在足球史上处于什么地位 全球要闻
【环球财经】伦敦金属交易所基本金属23日多数下跌_全球信息
白玉兰奖完整名单出炉 年初大热电视剧《狂飙》挂零陪跑-天天实时
中国高空开伞试验运载器发射连续成功:木星、天王星我们来了! 今日热讯
男子微信回了个“OK”表情 结果竟成被告!一点都不冤 速读
环球报道:尼康Z8新故障导致无法锁定镜头:官方承诺免费维修
【天天新要闻】读发布!设计与部署稳定的分布式系统(第2版)笔记10_自动化和缓慢的响应
无线路由器怎么连接电视(无线路由器怎么连接)
【环球新要闻】要考北大!汶川“敬礼娃娃”郎铮高考637分:15年前被埋20小时
关注:OPPO突然放弃自研芯片 真是因为没钱了?3000哲库人不信
看热讯:四川学霸女生高考712分查完分就睡觉、汶川“敬礼娃娃”郎铮637分
微软承认输掉“主机战争”:Xbox难以与竞争对手抗衡 每日头条
环球微动态丨特斯拉AI账号悄然上线:Dojo超级计算机下月开始生产
pro e
无牌产品硬刚国际大牌 就因为带货主播们买地建厂?
复兴号开进青藏铁路 提速至160公里/时 全程不到6小时
微软终于认怂!重新恢复Win11文件管理器经典功能