USACO算法系列二——fracdec

    题目在：http://www.nocow.cn/index.php/Translate:USACO/fracdec

    题目大意是，给出两个数，求出它们的商，小数中循环部分用括号括起来。

    刚拿到这道题的时候，没有很好的想法，觉得如果直接double表示商，然后提取每位小数位判断是否循环。但是这样会出现以下问题：由于在C语言中,双精度(double)型占8 个字节（64位）内存空间，其数值范围为1.7E-308～1.7E+308，可提供16位有效数字。 也就是他的小数位最多只有16位。

    因此，只能用长除法来解决。所谓长除法其实跟我们小学学的除法相似。

    我们知道只有当出现了曾经出现过的余数时，小数部分才会出现重复。重复的部分就是自从我们上次见到同样的余数之后计算出的部分。 
    因此，我们先读入并打印整数部分。接下来，我们对剩下的真分数部分进行长除直到我们发现了重复的余数或余数变为0。如果我们发现了重复的余数，即出现了循环节，就分别恰当地打印重复的部分和不重复的部分。如果余数变为0，即已经除尽，就打印整个小数部分。如果小数位根本没有被生成，那么打印一个0就是正确答案了。代码如下：

#include <iostream> #include <fstream> #include <math.h> #define MAX 100001 using namespace std; int n,d; int mod[MAX]={0}; int length = 0; int divide[MAX] = {0}; int idiv; int index = -1; char result[MAX]; int rl; int isMod(int m) { for (int i=0; i < length; i ++) { if (m == mod[i]) { return i; } } return -1; } int main() { ifstream fin("fracdec.in"); ofstream fout("fracdec.out"); fin >> n >> d; idiv = n / d; n = n % d; mod[0] = n; length ++; while (true) { int temp = mod[length-1] * 10; divide[length-1] = temp / d; int m = temp % d; if (m == 0) { index = length; break; } index= isMod(m); if ( index == -1) { mod[length] = m; length ++; } else { break; } } //打印结果 int t = 0; int temp = idiv; while (temp >= 10) { t ++; temp = temp / 10; } for (int i=t; i >= 0; i -- ) { int p = pow((double)10,i); result[rl] = idiv / p + '0'; idiv = idiv % p; rl ++; } result[rl] = '.'; rl ++; if (length==0) { result[rl] ='0'; rl ++; } else { int i=0; for (; i <index; i ++) { result[rl] = divide[i] + '0'; rl ++; } if (i < length) { result[rl] = '('; rl ++; for (; i < length; i++) { result[rl] = divide[i] + '0'; rl ++; } result[rl] = ')'; rl ++; } } for (int i=0; i <rl; i ++) { fout << result[i]; if (((i+1) % 76 )== 0 ) { fout << endl; } } fout <<endl; return 0; }    

    它的想法是计算循环开始前的小数位数，这样你甚至无需保存各个小数位和余数，程序的空间花费将大幅减小，而运行速度也能有所提高。我们知道2和5的幂是仅有的两种不导致循环的数，因此要找到循环前的各数位，我们只需分别找到分子分母中包含的因子2和5个数的差，再取两者的最大值（详见代码片段）。然后我们仅使用第一个余数，在计算时输出各个数位即可。

     刚开始觉得没有道理，但是通过自己亲自动手验证一下，发现确实是可以通过数学证明的。（为什么是2和5呢？因为小数可以当做一个数除以10得到的）。代码如下：

#include <iostream> #include <fstream> using namespace std; ofstream out("fracdec.out"); int colcount=0; int numBeforeRepeat(int n, int d) { int c2=0, c5=0; if (n == 0) return 1; while (d%2==0) { d/=2; c2++; } while (d%5==0) { d/=5; c5++; } while (n%2==0) { n/=2; c2--; } /* 可以抵消 */ while (n%5==0) { n/=5; c5--; } /* 可以抵消 */ if (c2>c5) if (c2>0) return c2; else return 0; else if (c5>0) return c5; else return 0; } void print (char c) { if (colcount==76) { out<<endl; colcount=0; } out<<c; colcount++; } void print (int n) { if (n>=10) print (n/10); print ((char)('0'+(n%10))); } int main() { int n, d; ifstream in("fracdec.in"); in>>n>>d; in.close(); print (n/d); print ('.'); n=n%d; int m=numBeforeRepeat(n,d); for(int i=0;i<m;i++) { n*=10; print (n/d); n%=d; } int r=n; if(r!=0) { print ('('); do { n*=10; print (n/d); n%=d; } while (n!=r); print (')'); } out<<endl; out.close(); exit (0); }