18444 分数拆分
时间限制:2500MS 代码长度限制:10KB
题型: 编程题 语言: G++;GCC;VC
Description
输入正整数k(k<=1000),将1/k变为不少于2项,但不多于3项的1/(xi)之和,xi为正整数,且i表示序号。
注:请使用long long。
输入格式
多case,一行一个整数k。
最后一行是0。
输出格式
对每一个case,按等式最右边一项分母,由小到大排序输出满足条件的等式,最右边一项分母相同,则按最右边第二项,依次类推。
每一个case完成后,输出一个空行(没有满足的等式时,也要输出该空行)。
输入样例
2
3
4
0
输出样例
1/2=1/6+1/3
1/2=1/42+1/7+1/3
1/2=1/24+1/8+1/3
1/2=1/18+1/9+1/3
1/2=1/15+1/10+1/3
1/2=1/12+1/12+1/3
1/2=1/4+1/4
1/2=1/20+1/5+1/4
1/2=1/12+1/6+1/4
1/2=1/8+1/8+1/4
1/2=1/10+1/5+1/5
1/2=1/6+1/6+1/6
1/3=1/12+1/4
1/3=1/156+1/13+1/4
1/3=1/84+1/14+1/4
1/3=1/60+1/15+1/4
1/3=1/48+1/16+1/4
1/3=1/36+1/18+1/4
1/3=1/30+1/20+1/4
1/3=1/28+1/21+1/4
1/3=1/24+1/24+1/4
1/3=1/120+1/8+1/5
1/3=1/45+1/9+1/5
1/3=1/30+1/10+1/5
1/3=1/20+1/12+1/5
1/3=1/15+1/15+1/5
1/3=1/6+1/6
1/3=1/42+1/7+1/6
1/3=1/24+1/8+1/6
1/3=1/18+1/9+1/6
1/3=1/15+1/10+1/6
1/3=1/12+1/12+1/6
1/3=1/21+1/7+1/7
1/3=1/12+1/8+1/8
1/3=1/9+1/9+1/9
1/4=1/20+1/5
1/4=1/420+1/21+1/5
1/4=1/220+1/22+1/5
1/4=1/120+1/24+1/5
1/4=1/100+1/25+1/5
1/4=1/70+1/28+1/5
1/4=1/60+1/30+1/5
1/4=1/45+1/36+1/5
1/4=1/40+1/40+1/5
1/4=1/12+1/6
1/4=1/156+1/13+1/6
1/4=1/84+1/14+1/6
1/4=1/60+1/15+1/6
1/4=1/48+1/16+1/6
1/4=1/36+1/18+1/6
1/4=1/30+1/20+1/6
1/4=1/28+1/21+1/6
1/4=1/24+1/24+1/6
1/4=1/140+1/10+1/7
1/4=1/42+1/12+1/7
1/4=1/28+1/14+1/7
1/4=1/8+1/8
1/4=1/72+1/9+1/8
1/4=1/40+1/10+1/8
1/4=1/24+1/12+1/8
1/4=1/16+1/16+1/8
1/4=1/36+1/9+1/9
1/4=1/18+1/12+1/9
1/4=1/20+1/10+1/10
1/4=1/15+1/12+1/10
1/4=1/12+1/12+1/12
问题分析:
对于拆分后的形式:1/A=1/C+1/D、1/A=1/B+1/C+1/D进行变形,以得到其中C或B的表达式(前式:C=(A×D)÷(D-A);后式:B=(A×C×D)÷(C×D-A×C-A×D))。若该表达式值为整数,即说明该数符合要求。其中先对前者(由两个分数组成的情况)进行检验,再判断三个分数的情况,并注意分母的大小次序问题。
代码实现:
#include<cstdio>
#include<iostream>
using namespace std;
int main()
{
long long a;
while(cin>>a&&a)
{
long long b,c,d;
for(d=a+1;d<=3*a;d++)
{
if((a*d)%(d-a)==0)//判断1/a=1/c+1/d中的c是否为整数
{
c=(a*d)/(d-a);
if(c>=d)
{
printf("1/%lld=1/%lld+1/%lld\n",a,c,d);
}
}
long long up=(a*d*2)/(d-a)+1;//这个上限可能是经验值
for(c=(a*d)/(d-a)+1;c<=up;c++)
{
if((a*c*d)%(c*d-a*c-a*d)==0)//判断1/a=1/b+1/c+1/d中的b是否为整数
{
b=(a*c*d)/(c*d-a*c-a*d);
if(b>=c&&c>=d)
{
printf("1/%lld=1/%lld+1/%lld+1/%lld\n",a,b,c,d);
}
}
}
}
cout<<endl;
}
return 0;
}