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

洛谷题单指南-进阶数论-P3861 拆分

原题链接:https://www.luogu.com.cn/problem/P3861

题意解读:将整数n拆分成不同因数之积的方案数,不含1*n的情况。

解题思路:

1、背景知识-超级合数

n的数据范围最大是10^12,尽管n很大,但是n以内的整数的约数个数最多是多少呢?

在数论中通常可以查询超级合数表来计算:超级合数介绍

查表可以,不超过10^12的超级合数是:

963761198400 2*2*2*2*2*2*3*3*3*3*5*5*7*11*13*17*19*23 6720

6720就是963761198400的约数个数,这个就是最多的约数个数。

有此基础,问题就好办了!

2、问题分析

对于一个整数n,先分解出所有的约数(包括1和n),存入数组a[],一共是cnt个,并对约数从小到大排序;

再定义一个数组pos[],用来存每个约数对应a中的下标,需要注意约数中在1~sqrt(n)内的不超空间,sqrt(t)~n范围内可能超空间,可以分开两个数组,对sqrt(t)~n范围的数进行hash处理,这里直接统一用map处理。

接下来用线性动态规划来解决,

状态定义:设f[i][j]表示将a[i]拆解成a[1]~a[j]之积的方案数,根据约数个数最大6720,数组开f[10000][10000]即可

状态转移:

对于f[i][j],主要看最后一个a[j]是否包含在a[i]拆解的约数之中

如果不包含a[j],必然不选a[j],则有f[i][j] = f[i][j - 1]

如果包含a[j],可以选也可以不选a[j],则有f[i][j] = f[i][j - 1] + f[pos[a[i] / a[j]]][j - 1]

初始化:f[1][1] = 1

结果:f[cnt][cnt] - 1,减一是因为排除1 * n的情况

100分代码:

#include <bits/stdc++.h>
using namespace std;typedef long long LL;
const int N = 10000, M = 1000000;
LL a[M], cnt;
unordered_map<LL, int> pos;
int f[N][N];
int t;int main()
{cin >> t;while(t--){cnt = 0;memset(f, 0, sizeof(f));LL n;cin >> n;for(int i = 1; i <= n / i; i++){if(n % i == 0){a[++cnt] = i;if(i != n / i) a[++cnt] = n / i;}}sort(a, a + cnt + 1);for(int i = 1; i <= cnt; i++) pos[a[i]] = i;f[1][1] = 1;for(int i = 1; i <= cnt; i++){for(int j = 2; j <= cnt; j++){f[i][j] = f[i][j - 1];if(a[i] % a[j] == 0) f[i][j] += f[pos[a[i] / a[j]]][j - 1];}}cout << f[cnt][cnt] - 1 << endl;}return 0;
}

 

http://www.gsyq.cn/news/12369.html

相关文章:

  • 线程的状态流转
  • AI工作流详解以及应用场景(AI)
  • 非结构网格中计算场梯度的手段比较
  • 前端模块化——彻底搞懂AMD、CMD、ESM和CommonJS
  • 实用指南:Java基础(十四):枚举类详解
  • VGGT: Visual Geometry Grounded Transformer
  • 微信小程序使用地图map 实现定位和实时绘画轨迹
  • 嵌入式入门,基于keil5用stm32寄存器和标准库实现LED流水灯
  • 小人鱼的数学题 - Li
  • CentOS将软件源修改为国内源
  • 【C++进阶】C++11 的新特性 | lambda | 包装器 - 实践
  • orcad放置器件时器件不在预览框中心
  • 从零开始:VirtualBox 虚拟机安装与 CentOS 7 部署 + 双网卡网络配置指南
  • 【光照】[物理模型]中的[BRDF]是什么?
  • 《Linux Robust锁》
  • Manim实现气泡特效
  • C# Inno Setup
  • CF2139虚拟游记
  • 融合多元定位技术,帮助应用破解精准定位难题
  • hutool主要内容list
  • Kurt-Blender零基础教程:第2章:建模篇——第3节:陈列/父子级/蒙皮/置换修改器与小狐狸角色建模 - 教程
  • 学习:uniapp全栈微信小程序vue3后台(26) - 指南
  • HTML5介绍(HTML5特性、HTML5功能) - 指南
  • 读书笔记:Oracle 自动索引:让数据库自己管索引?
  • 故障处理:Oracle RAC集群CTSS时钟同步故障案例分析与解决
  • PostgreSQL技术大讲堂 - 第106讲:分区表索引优化
  • AI智能体:从认知到实践
  • vue3小坑之-为什么把ref定义的数组赋值给数组对象后取值为空数组?
  • 【C++STL详解】带头双向循环结构 + 双向迭代器,核心接口 + 排序效率 + 避坑指南 - 教程
  • VBA ETH功能应用 | “0”代码构建SOME/IP节点