跪求C++或Pascal高精度

#include<iostream>
using namespace std;
string a,b;
int c[100]={0};
int main()
{
int i,j,k,s=0;;
cin>>a>>b;
for(i=1;i<=a.size();i++)
for(j=1;j<=b.size();j++)
{
c[i+j-1]+=(a[i-1]-'0')*(b[j-1]-'0');
k=i+j-1;
s=max(s,k);
while(c[k]>9)
{
c[k+1]+=c[k]/10;
c[k]=c[k]%10;
k++;
s=max(s,k);
}
}
bool shit=false;
if(b[b.size()-1]=='0'||a[a.size()-1]=='0')
for(i=1;i<=s;i++)
{
shit=true;
cout<<c[i];
}
else
for(i=s;i>=1;i--)
{
shit=true;
cout<<c[i];
}
while(1);