正の整数を読み込んで、素因数に分解するプログラムをお願いします。

[考え方]
　除数をｋとすると、まず、ｋに素因数となり得る最小の数2を設定する。
ｋがｎの約数かどうかを調べ、約数であれば、ｋの値を出力する。
このとき、n/ｋは、それまでに見つかっていない素因数の全ての積に
なるので、ｎの値をｎ/kの値で置き換える。
　同じ素因数がいくつもある場合があるので(上例の2および3）、ｋが
ｎの約数である限り、同じｋの値で、この処理を繰り返す。
　ｎがｋで割り切れなくなった時は、ｋの値を１増やし、以上に述べた処理を
繰り返す。このとき、3以上のｋの値については奇数だけを相手にすればよい。
また、一つのｎの値については、ｋが√ｎを越えない最大の整数になるまで
行えばよい（ｎの値は変わっていくことに注意しよう）。ただし、ｋが√ｎを越えた
ときには、ｎは１以外には√ｎ以下の約数をもたないことになるので、そのときの
ｎの値は素数であるか、１であるかのいずれかである。

sqrt((double)n)か(int)sqrt((double)n)を使うらしいのですが･･･。