当前位置: 首页 > news >正文

[杭电春季联赛5]1004 赛马

原帖地址:https://www.cnblogs.com/Reisentyan/p/19885859

[杭电春季联赛5]1004 赛马

我们将使用拉马努金瞪眼法解决这一题:
注意到,样例很有规律

考虑规律,当 \(r=x\) 时,对答案有多少贡献
手玩枚举发现:

r (十进制) 贡献
1 0
2 1
3 1
4 3
5 3
6 4
7 4
8 7

确实很有规律
注意到,题面中存在 \(log\),发现与二进制相关,于是将 \(r\) 转换成二进制:

r (二进制) 贡献
0001 0
0010 1
0011 1
0100 3
0101 3
0110 4
0111 4
1000 7

此时答案已经很清晰了,我们发现,当末尾突然多出几个后缀,做出的贡献就多出同样的数量
并且我们发现,做出的贡献等于:$$r - \text{popcount}(r)$$
也就是说\(1\)\(n\) 求和 减去 所有的二进制位 \(1\) 的数量 = \(ans\)

$ \sum_{i=1}^{n} i - \sum_{i=1}^{n} \text{popcount}(i) = \text{ans} $​

二进制位 \(1\) 的数量可以用数位DP算出

注意到,代码要写成这样:

void solve()
{ll n = q_;ll all = (n + 1) * n / 2;constexpr int N = 64;ll l, r, dp[N], mi[N];ll ans1[2] = { 0 }, ans2[N] = { 0 };int a[N];auto gett = [&](ll nn, ll* ans)->void{ll tmp = nn;int len = 0;while (nn) a[++len] = nn % 2, nn /= 2;for (int i = len; i >= 1; --i) {for (int j = 0; j < 2; j++) ans[j] += dp[i - 1] * a[i];for (int j = 0; j < a[i]; j++) ans[j] += mi[i - 1];tmp -= mi[i - 1] * a[i], ans[a[i]] += tmp + 1;ans[0] -= mi[i - 1];};};l = 1;r = n;mi[0] = 1ll;dp[0] = 0;for (int i = 1; i <= 62; ++i) {dp[i] = dp[i - 1] * 2 + mi[i - 1];mi[i] = 2ll * mi[i - 1];}gett(r, ans1), gett(l - 1, ans2);cout << all - ans1[1] + ans2[1] << endl;return;
}/*
⡀⠎⠀⠀⠀⠀⠀⠀⠀⣸⣿⣿⣿⣿⣄⠃⠈⣶⡛⠿⠭⣉⠛⠿⡿⠛⠉⣀⣠⣤⣭⡏⠴⢀⣴⣿⣿⣿⣿⣿⣿⠀⠀⠀⠀⠀⠀⠀⠙⣿⣿
⠀⠀⠀⠀⠀⠀⠀⠀⣿⣿⣿⣿⣿⣿⣷⣱⣬⠛⠉⠀⠀⢠⠀⠀⠀⢀⣀⠀⠉⠿⣿⣾⣿⣿⣿⣿⣿⣿⣿⣿⠀⠀⠀⠀⠀⠀⠀⠀⠈⡿
⠀⠀⠀⠀⠀⠀⠀⢀⢿⣿⣿⣿⣿⣿⣿⠋⠀⠀⠀⠀⠀⡏⠀⠀⠀⠀⠈⠳⠀⠀⠀⠻⣿⣿⣿⣿⣿⣿⠋⠀⣇⠀⠀⠀⠀⠀⠀⠀⠀⠈
⠀⠀⠀⠀⠀⠀⠀⣸⠀⣿⣿⣿⣿⠟⠀⠀⠀⠂⠀⠀⢠⠀⠀⠀⠀⠀⠀⠀⠈⡀⠀⠀⠀⠻⣿⣿⣿⣿⣷⡀⠘
⠀⠀⠀⠀⠀⠀⠀⣧⣿⣿⣿⣿⠋⠀⠀⠀⠀⠀⠀⠀⢸⠀⠀⠀⠀⠀⠀⠀⠀⠈⠀⠀⠀⠀⠙⣿⣿⣿⣿⣿⣄⣧
⠀⠀⠀⠀⠀⠀⣸⣿⣿⣿⣿⠁⠀⠀⠀⠀⠀⠀⠀⠀⣾⠀⠀⠀⠀⠀⠀⠀⠀⠀⢧⠀⠀⠀⠀⠈⢿⣿⣿⣿⣿⣿⣆
⠀⠀⠀⠀⠀⢀⣿⣿⣿⣿⠇⠀⠀⠀⠀⠀⠀⠀⠀⠀⢹⠀⠀⠀⠀⠀⠀⠀⠀⠀⠈⠀⠀⠀⠀⠀⢂⠻⣿⣿⣿⣿⣿⣄
⠀⠀⠀⠀⠀⣿⣿⣿⣿⣹⠀⠀⠀⠀⠀⢸⠀⠀⠀⠀⠸⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⣇⠀⠀⠀⠀⠀⡄⠈⢿⣿⣿⣿⣿⣆
⠀⠀⠀⠀⣿⣿⣿⣿⠁⡇⠀⠀⠀⠀⠀⢸⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⣿⠀⠀⠀⠀⠐⠸⠀⠀⠻⣿⣿⣿⣆⢦
⠀⠀⢠⣿⣿⣿⣿⠃⠀⠀⠀⠀⠀⠀⠀⣼⡇⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⡏⣧⠀⠀⠀⠀⠐⣇⠀⠀⠙⣿⣿⣿⡄⠙⣄
⠀⣴⣿⣿⣿⣿⠏⠀⢸⠀⠀⠀⠀⠀⠀⡿⢿⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⣃⣈⣦⠀⠀⠀⠀⢹⠀⠀⠀⠸⣿⣿⣿⠀⠀⠳⣀
⠋⣸⣿⣿⣿⡟⠀⠀⠀⡆⠀⠀⠀⠀⠀⡏⠙⡀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⢰⠀⢠⠀⠀⠀⢧⠀⠀⠀⠀⡇⠀⠀⠀⠘⣿⣿⣷⠀⠀⠘
⠀⣿⣿⣿⢩⠀⠀⠀⠀⣿⠀⠀⠀⠀⠀⣀⠀⢱⠀⡀⠀⠀⠀⠀⠀⠀⠀⠀⣿⠀⠂⢀⣴⣶⣿⣿⡀⠀⠀⢻⠀⠀⠀⠀⠹⣿⣿⡄
⢸⣿⣿⠃⠈⠀⠀⢸⠀⣿⣆⠀⠀⠀⠀⣿⣿⣿⠷⠘⡀⠀⠀⠀⠀⠀⠀⢠⢹⡀⠈⡿⠻⣿⣛⢿⣿⣷⡀⠈⠀⠀⠀⠀⠀⢻⣿⣿
⣿⣿⣿⠀⠀⠀⠀⢸⠀⡇⣼⣄⠀⠀⠀⢻⣿⡄⠑⠑⣿⡀⠀⠀⠀⢀⠀⠂⠇⠀⠀⠖⠛⢿⣿⣿⣌⢿⣿⣿⡆⠀⠀⠀⠀⠀⣿⣿⡀
⣿⣿⡇⠀⠀⠀⠀⢸⠀⣾⣿⣿⡷⠿⣷⣤⣿⣿⡄⠀⠀⠀⠑⠤⡀⠀⠃⠀⠀⠀⠀⣿⣶⣿⣿⣿⣿⣆⠙⣿⣧⠀⠀⠀⠀⠀⣿⣿⡇
⣿⣿⠁⠀⠀⠀⠀⠘⣾⣿⣿⠁⣴⣿⣿⣿⣿⣿⣇⠀⠀⠀⠀⠀⠀⠈⠀⠀⠀⠀⠀⠸⡏⠙⣿⠉⠻⣿⠀⠀⣿⠀⠀⠀⣄⠀⣿⢸⣷
⣿⣿⡇⠀⠀⠀⠀⠀⣿⣿⠁⠀⣿⣿⠋⣿⠏⠙⠇⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⢹⠀⢀⢻⠀⠀⢀⡟⢀⣿⣸⢃⠟
⣿⣿⣿⠀⡄⠀⠀⠀⠘⠻⡄⠀⢹⠁⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⡘⠀⢀⣿⠃⣿⣿⡗⠁
⣧⣿⣿⣧⢹⡀⠀⠀⠀⠱⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⢀⠀⣴⣿⣿⣾⣿⣿⣿
⢿⠘⣿⣿⣿⣿⣤⠀⠢⡀⠱⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀ ⠀⠀⡀⠀⠀⠀⠀⠀⠀⠀⠀⠀⣠⣵⣿⣿⣿⣿⣿⣿⣿⣿⣷
⠀⠉⣿⣿⣿⡿⣿⠻⣷⣬⣓⣬⣄⠀⠀⠀⠀⠀⠀⠀ ⠀⠀⠉⠈⠈⠈⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⢀⣾⠃⠼⢉⣿⣿⣿⣿⣿⣿⣿
⠀⠀⣿⣿⣿⣷⠀⠀⠀⠘⣿⣄⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⣠⣾⣿⡏⠀⠀⢸⠀⢻⢿⣿⣿⡏⣿
⠀⢸⣿⣿⣿⣿⠀⠀⠀⠀⢻⣿⣿⣤⣀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⠀⣀⣴⣾⣿⣿⣿⣿⠀⠀⠀⢸⠀⠀⢸⣿⣿⠘⡀
⢦⡿⣿⣿⣿⢿⠀⠀⠀⠀⢸⣿⣿⣿⣿⣿⣿⣿⣶⣶⣦⡄⠀⠀⠀⠀⠀⠀⠀⠀⣰⣿⣿⣿⣿⣿⣿⣿⣿⠀⠀⠀⠘⡄⠀⠈⣿⣿⡄⠱
⣴⠛⣾⣿⣿⢸⠀⠀⠀⠀⠀⣿⣿⣿⣿⣿⣿⣿⣿⣿⣿⠿⡄⠀⠀⠀⠀⠀⠀⠀⣯⠛⣿⣿⣿⣿⣿⣿⣿⠀⠀⠀⠀⣇⠀⠀⣿⣿⣿
⠿⠀⣿⣿⣿⠀⠀⠀⠀⠀⠀⣿⣿⣿⣿⣿⣿⣿⣿⠟⠰⡾⠃⠀⠀⠀⠀⠀⠀⠀⠙⡟⠀⢻⣿⣿⣿⣿⣿⡆⠀⠀⠀⠸⠀⠀⠸⣿⣿⣷
⠆⢳⣿⣿⡇⠀⠀⠀⠀⠀⠀⣿⣿⣿⠛⠿⠿⢿⡟⠀⠀⠉⠦⣀⡤⢶⠀⠖⠲⠶⠊⠀⠀⠀⢻⡛⠛⠛⣿⣿⠀⠀⠀⠀⠃⠀⠀⢿⣿⣿
*/
http://www.jsqmd.com/news/657961/

相关文章:

  • CMake实战指南:利用FetchContent优雅集成GitHub热门库
  • STM32LL库实战入门:从零搭建高效开发环境
  • gInk多显示器使用教程:如何在多个屏幕上完美标注
  • Hermes Agent横空出世!开源智能体新里程碑,轻松超越OpenClaw龙虾
  • 题解:AcWing 3646 分水果
  • 维普论文AI率60%怎么办?2026年这3款降AI工具帮你降到10%以下 - 我要发一区
  • Windows 10/11下FFmpeg调用NVIDIA显卡加速视频转码全攻略(含驱动版本检查)
  • Gumbo-Parser持续集成优化:测试时间缩短50%的终极指南
  • 别再用SonarQube跑规则了!2026奇点大会实测:LLM-native审查工具对逻辑漏洞识别率提升6.8倍(附12类业务逻辑缺陷特征库)
  • mysql如何通过Docker快速搭建_mysql容器化部署实践
  • puqk实名一个2025
  • 如何快速上手Kaf:从零开始的Kafka集群管理教程
  • Flutter ShadcnUI核心组件深度解析:30+精美UI元素一览
  • 2026长沙整装怎么选?权威选购指南与深度测评 - 品牌策略主理人
  • 别再让布线拖后腿!手把手教你用AXI Register Slice给Zynq设计提频(附Vivado配置避坑点)
  • 别再只用命令流了!用Workbench表格功能动态控制ANSYS流体渗透压力阈值
  • Redis 配置指南
  • RealWorld SvelteKit:终极全栈博客平台完整指南
  • NoSQL数据库Redis(二):Redis持久化详解
  • 01华夏之光永存:黄大年茶思屋榜文解法「第7期1题」OXC超快速切波技术·双路径解法
  • 互信息神经估计:从理论到实践的深度解析
  • 从PPT到产线:2026奇点大会AI重构建议的6步工业化落地路径,已验证缩短实施周期47%
  • 信号处理实战:用Python的SciPy库快速搞定傅里叶变换与拉普拉斯变换(附代码)
  • Linux 的 pwd 命令
  • 告别盲目调管子!用gm/ID方法在Cadence Virtuoso里搞定模拟IC设计(附SMIC 13nm工艺库仿真脚本)
  • 实测好用!Z-Image-Turbo-辉夜巫女快速体验,8步生成高质量辉夜巫女风格图
  • mcp-obsidian 最佳实践:7个实用技巧提升你的工作流效率
  • 终极指南:使用gumbo-parser轻松解析HTML5动态内容的10个技巧
  • 题解:洛谷 B2124 判断字符串是否为回文
  • TypeScriptToLua核心原理解析:深入理解AST转换与代码生成机制