POJ 1906 Three powers G++ 找规律 背

#include <iostream>
#include <cstring>
#include <string>
#include <algorithm>
using namespace std;
//英语 看博友分析  抄博友程序  大数乘法  找规律 背 
string mul(string a, string b)
{
	reverse(a.begin(),a.end());
	reverse(b.begin(),b.end());
	int ta[200];
	int tb[200];
	int tc[200];
	memset(ta,0,sizeof(ta));
	memset(tb,0,sizeof(tb));
	memset(tc,0,sizeof(tc));
	for(int i=0;i<a.size();i++)
	{
		ta[i]=a[i]-'0';
	}
	for(int i=0;i<b.size();i++)
	{
		tb[i]=b[i]-'0';
	}
	for(int i=0;i<a.size();i++)
	{
		for(int j=0;j<b.size();j++)
		{
			tc[i+j]=tc[i+j]+ta[i]*tb[j];
		}
	}
	int jw=0;
	for(int i=0;i<199;i++)
	{
		int t=tc[i];
		tc[i]=(t+jw)%10;
		jw=(t+jw)/10;
	}
	int i;
	for(i=198;i>0;i--)
	{
		if(tc[i]!=0)
		{
			break;
		}
	}
	string jg;
	for(;i>=0;i--)
	{
		jg=jg+(char)(tc[i]+'0');
	}
	return jg;
}
int main()
{	
	string da[80];//抄博友程序 
	da[0]="1";
	for(int i=1;i<=70;i++)
	{
		da[i]=mul(da[i-1],"3");
		//cout<<da[i]<<endl;
	}
	unsigned long long a;//  unsigned __int64

	while(cin>>a)
	{
		if(a==0)
		{
			break;
		}
		a=a-1;		
		cout<<"{ ";
		int ws=0;
		while(a!=0)
		{
			if(a&1)
			{
				cout<<da[ws];
				if(a>>1)
				{
					cout<<", ";
				}
			}
			ws++;
			a=a>>1;
		}
		cout<<" }"<<endl;
	}	
	return 0;
}