高精度阶乘的运算

阶乘：1x2x3x4.....N，仿照2的N次方的手法，只不过这次从前往后计算，得到的数从左往后，依次为个位十位百位...等等。
例如：021，表示120

#include <iostream>
using namespace std;
#define max 100
/*保存结果的数组长度*/
int f[max];
int main()
{
     int n=10; /*n表示n!，自由修改*/
     f[0] = 1;
     for(int i = 2;i <= n;i++)
     {
          int c = 0;/*进位值*/
          for(int j = 0;j < max;j++)/*整个数组乘一遍*/
              {
                   int s = f[j] * i + c;
                   f[j] = s % 10;
                   c = s / 10;
             }
     }
     for(int k = max - 1;k >= 0;k--) /*扫描清除右边的0*/
     if(f[k]) break;
     for(i = k;i >= 0;i--) /*从右往左打印结果*/
     cout<<f[i];
     return 0;
}

这个算法的一个缺点是，每次都要把数组乘一遍，有什么办法能优化它呢？