これまでのiwaoさんの書き込みから判断すると、アセンブリ言語にも精通されている方のように思います。
アセンブリ言語を利用すれば、更に高速な処理を記述できます。
もしかしたら、既にご存知かもしれません。そのときはご容赦ください。
ＢＣＤをご存知でしょうか。Binary Coded Decimal（２進化１０進数）といいます。
これはメモリ効率が悪いデータ形式ですが、今回のような処理にはもってこいです。
というか、これ以外に使い道があるのかどうか、疑問です。だれかご存知の方いらっしゃいます？
初めてＢＣＤを知ったとき「○ンテルはどうしてこんな命令組み込んだんだろう」と首を傾げてしまいました（爆）
それは置いといて、、、
ＣＰＵの命令だから動作速度は最速（のハズ）です。また、除算を必要としません。
（ただし動作速度については検証しておりません。悪しからず。）
ＢＣＤは次のようなフォーマットです。
２進化１０進数（ＢＣＤ整数）は、範囲０～９の有効値を持つ符号なし４ビット整数である。
ＩＡ－３２アーキテクチャは、１つ以上の汎用レジスタ内または１つ以上のｘ８７ ＦＰＵレジスタ内にあるＢＣＤ整数の演算を定義している。
『IA-32 インテル アーキテクチャ・ソフトウェア・デベロッパーズ・マニュアル 上巻：基本アーキテクチャ（資料番号 253665-013J）第４章 データ型』より引用
簡単に言えば、１０進数の一桁を表現するために４ビット使う、ということです。
例えば、１０進数の５２をＢＣＤで表現すると、0x52となります。
ｘ８７ ＦＰＵレジスタは１つあたり１０バイトの容量を持ちます。
先頭のバイトは符号バイトとして扱われるので、実質９バイトが利用できます。
ビットに換算すると７２ビット、一桁４ビットだから、７２÷４＝１８桁の１０進数をＢＣＤで表現できます。
この命令を利用する意義は大きいはずです。
こんな感じでいいとおもいます。
コード: 全て選択
unsigned int i = 98765; // 元の整数
unsigned char bcd[10]; // ＢＣＤ整数を格納するバッファ
char d[18]; // 各桁の値を格納する配列 d[k]は 10 ^ k の位の値を保持する
int n;
__asm {
fild i
fbstp bcd
}
for(n = 0; n < 9; n++) {
d[n * 2] = bcd[n] & 0x0F;
d[n * 2 + 1] = (bcd[n] & 0xF0) >> 4;
}
FILDは整数をST(0)にロードする命令のニーモニック
FBSTPはST(0)をＢＣＤ形式に変換してストアしたあと、STをポップする命令のニーモニック
このアルゴリズムをＤＬＬ化するとよいとおもいます。
さすがに、コードの切れ端や実践コードモジュールに掲載できる程度を逸脱してしまいましたが・・・
［追記］
6月11日 12:18
上記コードは汎用形ですが、次のような無駄もあります。
- forループ内で乗算を２回行っている
- Long型の最大値は +2,147,483,647 であり、１８桁を下回る
このような観点から、次のように最適化できます。
コード: 全て選択
int n2;
for(n = 0; n < 5; n++) { // 下位１０桁を有効に扱えるよう、カウンタを調整
n2 = n << 1; // ２倍にはシフト演算を利用
d[n2] = bcd[n] & 0x0F;
n2++;
d[n2] = (bcd[n] & 0xF0) >> 4;
}