高精度专辑（1码题库）
**10109 : 完美三角形**
近日，我们的陈老师在网上购买了一批木棍，这些木棍堪称完美，每一根长度大小都一样，且长度为1.陈老师买回了这些木棍之后就不管了，直接把任务交给了我们的恺成同学。任务是什么呢？需要我们的恺成同学编写一个程序，程序的功能就是每当用户查询一个数n,表示我们有有一个边长为n的等边三角形，如图所示.该程序需要计算出完美三角形的个数.完美三角形就是边长为1的三角形.
输入格式
输入一个整数n(0<n<=10^18),为等边三角形的边长;
输出格式
输出完美三角形的个数
样例输入
3
样例输出
9
这道题怎么做？？？很简单，就是高精度乘法>>>>>
代码如下：
#include<bits/stdc++.h>
using namespace std;
const int maxn=10005;
int na[maxn],nb[maxn],nc[maxn];
string mul(string a,string b)
{
int la=a.size(),lb=b.size();
for(int i=0;i<la;i++) na[la-1-i]=a[i]-'0';
for(int i=0;i<lb;i++) nb[lb-1-i]=b[i]-'0';
for(int i=0;i<la;i++)
for(int j=0;j<lb;j++)
nc[i+j]+=na[i]*nb[j];
for(int i=0;i<=la+lb;i++)
{
//cout<<nc[i];
nc[i+1]+=nc[i]/10;
nc[i]%=10;
}
// cout<<endl;
int len=la+lb;
string ans;
while(nc[len]==0&&len) len--;
for(int i=len;i>=0;i--) ans+=nc[i]+'0';
return ans;
}
int main()
{
string a,b;
cin>>a;
cout<<mul(a,a)<<endl;
return 0;
}
好了，解决了一道水提很难 很难 很难 很难 的题；
再来看下一题
**10176 : 大整数GCD**
问题描述
给出两个数a,b，求GCD(a,b)。
输入格式
输入两个数a,b(1<=a<=10100,1<=b<=10100)。
输出格式
输出GCD(a,b)。
样例输入
3 6
样例输出
3
是不是很水呢？直接上代码。
#pragma GCC optimize(2)
#include<iostream>
#include<cstdio>
#include<cstring>
#define inf 1000000000
using namespace std;
char ch1[10005],ch2[10005];
int la,lb,cnt;
struct data{int a[1205],l;}a,b;
bool com()
{
if(a.l<b.l)return 0;
if(a.l>b.l)return 1;
for(int i=a.l;i>0;--i)
if(a.a[i]>b.a[i])return 1;
else if(a.a[i]<b.a[i])return 0;
return 1;
}
void print(data a)
{
while(a.a[a.l]==0)a.l--;
for(int i=a.l;i>0;--i)
if(i==a.l)printf("%d",a.a[i]);
else printf("%09d",a.a[i]);
}
inline data sub(data a,data b)
{
int k;
data c;
for(int i=1;i<=1200;++i)
{
if(i<=b.l)c.a[i]=a.a[i]-b.a[i];
else if(i<=a.l)c.a[i]=a.a[i];
else c.a[i]=0;
if(c.a[i]<0)
{
c.a[i]+=inf;
a.a[i+1]--;
}
}
c.l=a.l;
while(c.a[c.l]==0&&c.l)c.l--;
return c;
}
void diva()
{
for(int i=1;i<=a.l;i++)
{
if(a.a[i]&1)a.a[i-1]+=inf/2;
a.a[i]>>=1;
}
if(!a.a[a.l])a.l--;
}
void divb()
{
for(int i=1;i<=b.l;i++)
{
if(b.a[i]&1)b.a[i-1]+=inf/2;
b.a[i]>>=1;
}
if(!b.a[b.l])b.l--;
}
void mul()
{
for(int i=a.l;i>0;i--)
{
a.a[i]<<=1;
a.a[i+1]+=a.a[i]/inf;
a.a[i]%=inf;
}
while(a.a[a.l]>0)a.l++;
for(int i=b.l;i>0;i--)
{
b.a[i]<<=1;
b.a[i+1]+=b.a[i]/inf;
b.a[i]%=inf;
}
while(b.a[b.l]>0)b.l++;
}
int main()
{
scanf("%s%s",ch1+1,ch2+1);
la=strlen(ch1+1);lb=strlen(ch2+1);
if(la%9)a.l=la/9+1;
else a.l=la/9;
if(lb%9)b.l=lb/9+1;
else b.l=lb/9;
for(int i=1;i<=a.l;++i)
{
int k1=max(1,la-i*9+1),k2=la-(i-1)*9;
for(int j=k1;j<=k2;++j)
a.a[i]=a.a[i]*10+ch1[j]-'0';
}
for(int i=1;i<=b.l;++i)
{
int k1=max(1,lb-i*9+1),k2=lb-(i-1)*9;
for(int j=k1;j<=k2;++j)
b.a[i]=b.a[i]*10+ch2[j]-'0';
}
while(1)
{
if((a.a[1]%2==0)&&(b.a[1]%2==0)){diva();divb();cnt++;}
else if((a.a[1]%2==0))diva();
else if((b.a[1]%2==0))divb();
if(com()){a=sub(a,b);if(!a.l){while(cnt--)mul();print(b);break;}}
else {b=sub(b,a);if(!b.l){while(cnt--)mul();print(a);break;}}
}
return 0;
}
可能我写的有些繁琐，
不过还是ac了，
还有本人码字速度慢，
请各位原谅。
10177 : 大整数取模1
又是一水提nananan题ing；
原题君~~~
问题描述
给出两个数a,b，输出a对b取模的结果。
输入格式
输入两个正整数a,b。a<=10100000，b<=109。
输出格式
输出a对b取模的结果。
样例输入
10 3
样例输出
1
因为是高精除低精，各大教科书上都有，就不多BB，上代码
#include<bits/stdc++.h>
using namespace std;
int main()
{
char n[100005];
int m;
scanf("%s%d",n,&m);
long long len = strlen(n);
long long ans = 0;
for (int i = 0; i < len; i++)
ans = (int)((((long long)ans * 10) % m + (n[i] - '0') % m) % m);
printf("%d\n",ans);
return 0;
}
小二一开始还TLE了呢！！！！
以下是错误代码
#include<iostream>
#include<cstring>
using namespace std;
int a[100001],b[1001],c[1001],d,i;
void init(int a[])
{ string s;
cin>>s;
a[0]=s.length();
for(i=1;i<=a[0];i++)
a[i]=s[a[0]-i]-'0';
}
void print(int a[])
{
int i;
if (a[0]==0){cout<<0<<endl;return;}
for(i=a[0];i>0;i--) cout<<a[i];
cout<<endl;
return ;
}
int compare (int a[],int b[])
{ int i;
if(a[0]>b[0]) return 1;
if(a[0]<b[0]) return -1;
for(i=a[0];i>0;i--)
{if (a[i]>b[i]) return 1;
if (a[i]<b[i]) return -1;
}
return 0;
}
void jian(int a[],int b[])
{ int flag,i;
flag=compare(a,b);
if (flag==0) {a[0]=0;return;}
if(flag==1)
{for(i=1;i<=a[0];i++)
{if(a[i]<b[i]){ a[i+1]--;a[i]+=10;}
a[i]-=b[i];
}
while(a[0]>0&&a[a[0]]==0) a[0]--;
return;
}
}
void numcpy(int p[],int q[],int det)
{
for (int i=1;i<=p[0];i++) q[i+det-1]=p[i];
q[0]=p[0]+det-1;
}
void chugao(int a[],int b[],int c[])
{
int i,tmp[101];
c[0]=a[0]-b[0]+1;
for (i=c[0];i>0;i--)
{
memset(tmp,0,sizeof(tmp));
numcpy(b,tmp,i);
while(compare(a,tmp)>=0){c[i]++;jian(a,tmp);}
}
while(c[0]>0&&c[c[0]]==0)c[0]--;
return ;
}
int main()
{
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
memset(c,0,sizeof(c));
init(a);init(b);
chugao(a,b,c);
print(a);
return 0;
}
这代码是我上课时老师讲的敲下来的，果然错了，还是得自己去思考呀！！！
