多倍長演算
数値にはunsigned short a[N+1];のような配列を用いる.ここでa[0]が整数部分、
a[1]からa[N]までが小数部分で、32768(=32768)進法で表す.たとえば３という値は
a[0]=3,a[1]=...=a[N]=0で表す.log[10]32768≒4.5154であるからｍ配列要素1個分が
10進4.5桁に相当する.したがって、10進で小数第1000位まで求めるには、N≒1000/4.5≒
222とする.（〜〜うんたらかんたら）
とか書いてあります。
「Ｃ言語による最新アルゴリズム事典」(P.179)より抜粋
本読めよ。