第十二届蓝桥杯C/C++ B组试题 D: 货物摆放

试题 D: 货物摆放
本题总分：10 分

【问题描述】
小蓝有一个超大的仓库，可以摆放很多货物。
现在，小蓝有 n 箱货物要摆放在仓库，每箱货物都是规则的正方体。小蓝
规定了长、宽、高三个互相垂直的方向，每箱货物的边都必须严格平行于长、
宽、高。
小蓝希望所有的货物最终摆成一个大的立方体。即在长、宽、高的方向上
分别堆 L、W、H 的货物，满足 n = L × W × H。
给定 n，请问有多少种堆放货物的方案满足要求。
例如，当 n = 4 时，有以下 6 种方案：1×1×4、1×2×2、1×4×1、2×1×2、
2 × 2 × 1、4 × 1 × 1。
请问，当 n = 2021041820210418 （注意有 16 位数字）时，总共有多少种
方案？
提示：建议使用计算机编程解决问题。
【答案提交】
这是一道结果填空的题，你只需要算出结果后提交即可。本题的结果为一
个整数，在提交答案时只填写这个整数，填写多余的内容将无法得分。

答案：2430
思路：先获得2021041820210418所有质因数，再通过质因数去组合从而获得所有的正约数，最后只需在所有的正约数找3个乘积为2021041820210418就行。
code

#include<iostream>
#include<set>
using namespace std;
const long long N = 2021041820210418; 

long long a[100],len = 0; //a数组存储质因数 len存储个数 
set<long long> s;//存储正约数 
long long ans = 0;//答案 

void get_primes(long long n){//获得N所有的质因数 
	for(int i = 2; i <= n/i; i++){
		while(n % i == 0) n/=i,a[len++] = i;
	}
	if(n > 1) a[len++] = n;
}

//需要选n个不同位置的数 已选了m个数 
//上一层循环在k位置 j为 m 个数的乘积 
void fun(int n,int m,int k,long long j){ 
	if(m == n + 1){ //已选n个数 
		s.insert(j);
		return;
	}  
	for(int i = k+1; i < len; i++)
		fun(n, m+1, i, j*a[i]);
}

int main(){
	//2021041820210418
	get_primes(N);
	cout<<"所有质因数："; 
	for(auto x : a) if(x != 0) cout<<x<<" "; 
	
	//选i个质因数组合正约数 获得所有的正约数使用set去重 
	for(int i = 1; i <= len; i++) fun(i,1,-1,1); 
	s.insert(1);
	
	cout<<"所有正约数："<<s.size()<<endl;
	for(auto x: s)	cout<<x<<endl;
	
	cout<<"正约数个数："<<s.size()<<endl;
	
	//开始组合正约数 
	for(auto a: s)
		for(auto b: s)
			for(auto c: s)
				if(a*b*c == N) ans++;		
	cout<<"答案："<<ans;
	return 0;
} 

更简单的做法，我们可以发现直接求正因数就行了，不用上面那样拆成质因数再组成正因数。

#include<iostream>
#include<algorithm>
#include<set>
#include<unordered_set> 
#include<vector>
#include<cmath>
using namespace std;
const int INF = 0x3f3f3f3f;
int main(){
	long long N = 2021041820210418;
	int res = 0;
	set<long long> ans;
	for(long long  i = 1; i*i <= N; i++){
		if(N % i == 0) {
			ans.insert(i);
			ans.insert(N/i);
		}
	} 
	for(auto i:ans)
		for(auto j:ans)
			for(auto k:ans)
				if(i*j*k == N) res++;
	cout<<res;			
	return 0;
}