如何解决,详见程序。
//#include "stdafx.h" //VS 默宴液认使用
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#pragma warning(disable: 4996)
//为了最大限晌友物度支持(目前)小数,整数定义成:__int64
//支持32位:以下分别为 int(__int32) 1E8
#define ZHENGSHU __int64
#define XIAOSHUDIANZUOYI 1E15 //double 最多只能确精确到15位有效数字
#define XIUZHENG 1E-16 //用于修正浮点数计算误差
int _tmain(int argc, _TCHAR* argv[])
{
//怎样用c语言写一个将小数转化成最简分数的算法
//开始
double xs = 0.000625//为了简便,直接指定。也可以让用户输入
//scanf("%lf", &xs)
ZHENGSHU fz, fm, zs, k, fh
//数学算法是:如 0.123456 = 123456/1000000 ,然后约分
//实现第一步转换
if (xs <0)
{
fh = -1
xs = -xs
}
else
{
fh = 1
}
zs = (ZHENGSHU)xs//支持假分数
//计算整数部分的有效数位数
k = 1
while (zs / k >0)
k *= 10
xs -= (double)zs
fm = XIAOSHUDIANZUOYI / k //分母
xs = (double)((ZHENGSHU)((xs + XIUZHENG * k) * fm)) / fm//修正浮点数计算误差
xs += XIUZHENG//修正浮点数计算误差
fz = (ZHENGSHU)(xs * fm)//分子
for () //这样循环,是因为 C/C++ 中,for 循环效率最高
{ //约去多乘的 10 的倍数
if ((fz % 10 == 0) &&(fm % 10 == 0))
{
fz /= 10
fm /= 10
}
else
{
break
}
}
//实现第二步:约分。以下看起来有点麻烦,目的是为了减少算法的时间复杂的
while ((fz % 2 == 0) &&(fm % 2 == 0))
{ //将公约数 2 约尽
fz /= 2
fm /= 2
}
k = 3
for ()
{
while ((fz % k == 0) &&(fm % k == 0))
{
fz /= k
fm /= k
}
k += 2
if (k >(ZHENGSHU)sqrt(fz))
break
}
//将整数部分加上,形成假分数。如告核果原数为不为0整数,则化成分母为 1 的假分数
//如果原数为 0,则化成分子为 0,分母为 1。
if (fz == 0)
{
fz = fh * zs
fm = 1
}
else
{
fz = fh * (fz + fm * zs)
}
//结束
//显示一下结果
printf("%lf = %lld/%lld\n", fh * (zs + xs), fz, fm)
//printf("%lf = %ld/%ld\n", fh * (zs + xs), fz, fm)//32位
system("pause")//防止窗口一闪而退
return 0
}
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)