两个长整数相乘

#include <iostream>

#include <iomanip>

using namespace std;

 

 

#define MAX_USHORT 65536

#define MAX_UCHAR  256

 

void Foo(unsigned long fry,unsigned char *szVer)

{

   //f＝2814749×fry＋4×fry/5 fry为unsigned long型，

   // fry范围为100至30000000；

   //2^16=65536

   unsigned long val[4]={0};

   unsigned short tFry[2],tCONST[2];

   unsigned long temp=4*fry/5;

   int i,j;

   tCONST[1]=2814749/MAX_USHORT;

   tCONST[0]=2814749%MAX_USHORT;

   tFry[1]=fry/MAX_USHORT;

   tFry[0]=fry%MAX_USHORT; 

  

   val[1]=temp/MAX_USHORT;

   val[0]=temp%MAX_USHORT;

   for (i=0;i<2;++i)

   {

      for (j=0;j<2;++j)

      {

         temp=tFry[i]*tCONST[j];

            val[i+j]+=(temp%MAX_USHORT);

            val[i+j+1]+=(temp/MAX_USHORT);

 

         val[i+j+1]+=val[i+j]/MAX_USHORT;

         val[i+j]%=MAX_USHORT;

      }

   }

   j=5;

   for (i=0;i<3;++i)

   {

      szVer[j]=val[i]%MAX_UCHAR;

      --j;

      szVer[j]=val[i]/MAX_UCHAR;

      --j;

   }

 

}

 

int main()

{ 

   unsigned  char szVer[6]={0};

   unsigned long fry;

   cin>>fry;

   unsigned long f=2814749*fry+4*fry/5;

   unsigned char*p=(unsigned char *)&f;

   for(int i=0;i<4;++i)

   {

      cout<<setbase(16)<<setw(5)<<(unsigned long)p[i];

   }

   cout<<endl;

   Foo(fry,szVer);

   for (int i=0;i<6;i++)

   {

      cout<<setbase(16)<<setw(5)<<(unsigned long)szVer[i];

   }

   cout<<endl;

   return 0;

}