从蛮力到策略—百鸡问题与完美立方

彼岸无生花

已于 2023-04-12 00:45:35 修改

阅读量578

点赞数 3

文章标签： c++ 算法

于 2023-04-12 00:20:11 首次发布

本文链接：https://blog.csdn.net/weixin_63530698/article/details/130096712

版权

（1）百鸡问题

“鸡翁一，值钱五；鸡母一，值钱三；鸡雏三，值钱一。百钱买百鸡，问鸡翁、母、雏各几何？”

：在不甚思考的情况下(每种鸡都可能有0-100只)凭直觉写出蛮力法求解百鸡问题的基本思路并编程实现，思考该算法复杂度的渐进表达式。

设鸡翁为x，鸡母为y，鸡雏为z，根据问题要求我们就能有以下两个公式：

5x+3y+1/3z=100

x+y+z=100

然后我们就能根据这个数学模型来编写代码

：缩小三种鸡数量可能的范围，优化算法，减少计算量，此时算法复杂度是否减少？

缩小鸡的范围，我们可以通过简单的数学计算得到：鸡翁百钱最多买20只，鸡母最多买33只，鸡雏最多100只，这样我们就能给代码加上条件，提升运算速度。

：进一步优化算法，能否降低算法的复杂度？请写出思路并编程实现。

进一步思考，我们发现z（鸡雏）可以用x和y表示（这类问题最好把非整数转换成其他数表示），这样我们就能得到：z=100-x-y

代码粘贴区：

1.#include<iostream>
using namespace std;


int main(){
	int x,y,z;
	for(x=0;x<=100;x++)
		for(y=0;y<=100;y++)
			for(z=0;z<=100;z++)
				if((x+y+z==100)&&(5*x+3*y+z/3)==100&&z%3==0)
					cout<<"x="<<x<<" "<<"y="<<y<<" "<<"z="<<z<<endl;							      					
	return 0;
}

2.#include<iostream>
using namespace std;

//鸡翁为x，鸡母为y，鸡雏为z。 
int main(){
	int x,y,z;
	for(x=0;x<=20;x++)
		for(y=0;y<=33;y++)
			for(z=0;z<=100;z++)
				if((x+y+z==100)&&(5*x+3*y+z/3)==100&&z%3==0)
					cout<<"x="<<x<<" "<<"y="<<y<<" "<<"z="<<z<<endl;							      					
	return 0;
}

3.#include<iostream>
using namespace std;

//鸡翁为x，鸡母为y，鸡雏为z。 
int main(){
	int x,y,z;
	for(x=0;x<=20;x++)
		for(y=0;y<=33;y++){
			z=100-x-y; 
			if((5*x+3*y+z/3)==100&&z%3==0)
			cout<<"x="<<x<<" "<<"y="<<y<<" "<<"z="<<z<<endl;
		} 
			
											      					
	return 0;
}

运行结果截图：

：调整输入规模(如1000鸡1000钱问题、2000鸡2000钱问题等)，观察并记录至少5组(每组多次运行取平均值)不同输入规模下任务1和任务3算法的运行时间，制作excel图表进行对比。

（2）掌握枚举算法的基本思想和实现方法。

任务描述：如a3= b3 + c3 + d3的等式被称为完美立方等式。例如 123= 63 + 83 + 103 。编写一个程序，对任给的正整数N (N≤100)，寻找所有的四元组(a, b, c, d)，使得等式成立，其中a,b,c,d 大于 1, 小于等于N，且 b<=c<=d。

输入：一个正整数N (N≤100)。 l

输出：每行输出一个完美立方。输出格式为：a^3=b^3+c^3+d^3。其中a,b,c,d所在位置分别用实际求出四元组值代入。

请按照a的值，从小到大依次输出。当两个完美立方等式中a值相同时，b值小的优先输出，若b仍相同，c值小的优先输出，c再相同时d值小的先输出。

样例输入

样例输出

6^3=3^3+4^3+5^3

12^3=6^3+8^3+10^3

18^3=2^3+12^3+16^3

18^3=9^3+12^3+15^3

19^3=3^3+10^3+18^3

20^3=7^3+14^3+17^3

24^3=12^3+16^3+20^3

代码粘贴区：

#include<iostream>
using namespace std;

int main(){
	int a,b,c,d;
	int N;
	cout<<"输入N";
	cin>>N;
	for(a=2;a<=N;a++){
		for(b=2;b<=a;b++)
			for(c=b;c<=a;c++)
				for(d=c;d<=a;d++){
					if(a*a*a==b*b*b+c*c*c+d*d*d)
					cout<<a<<"^3="<<b<<"^3+"<<c<<"^3+"<<d<<"^3"<<endl;
			}
			
		
	}
	return 0;
}

运行结果截图：