我写的高精度乘法有bug,请大家看看

οしαУ′颴

#include<iostream>
#include<cstring>
using namespace std;
char a[100]={0};
char b[100]={0};
int e;
int main()
{
while(cin>>a)
{
int c[100]={0};
int alen,blen,clen;
int temp,j,k;

alen=strlen(a)-1;
for(int i=0;i<alen/2;i++)//倒置
{
temp=a[i];
a[i]=a[alen-i];
a[alen-i]=temp;
}
strcpy(b,a);
blen=alen;
for( j=0;j<=alen;j++)//核心部分
for( k=0;k<=blen;k++)
{
c[j+k]+=(a[k]-48)*(b[j]-48);
c[k+j+1]+=c[j+k]/10;
c[j+k]%=10; //结果保存在c数组中
}
clen=j+k+1;
while(c[clen]==0)
clen--;
if(clen<0)
{
c[0]=0;
clen=0;
}
for(int t=0;t<=clen/2;t++)//c倒置
{
temp=c[t];
c[t]=c[clen-t];
c[clen-t]=temp;
}
for(int t1=0;t1<=clen;t1++)
cout<<c[t1];
cout<<endl;

memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
memset(c,0,sizeof(c));
}
system("pause");
return 0;
}

我想实现的是高精度平方运算..
bug:
不管是用VC还是用DEV C++编译运行,在计算12*12,15*15...13*13...70*70之类的都不能得到正确的结果..
但是有些数..例如三位数,或者是123456789之类都可以
很头痛...我已经把代码注释都写了..请高手帮我看看...告诉我是哪里出错了.
请帮我看看,错误在哪...我头都快爆了..
另,如果你们有比较清晰明了的高精度的算法教程(自然语言描述的,不要给我代码)请传给我一份
谢谢大家