#前言
话说我真的只能颓一颓[普及-]难度的内容了…药丸。
而且以前还在被一本通毒害时期连高精度取模和高精除高精都不会…我是废柴啊…
因为考完月考后干劲为0,不想码程序,所以还是先从高精度下刀来恢复手感吧。(像abs和Angelina那样的大神都是从平衡树和数学开始的)
高精度加法
洛谷P1610
模拟竖式加法。
#include<iostream>
#include<cstdio>
#include<iomanip>
#include<algorithm>
#include<cstring>
using namespace std;
int a[501],b[501],c[505];
char sa[501],sb[501];
int main()
{
int i,j,na,nb,nc,x;
scanf("%s\n",sa);
scanf("%s",sb);
na=strlen(sa);nb=strlen(sb);
for(i=0;i<na;i++)a[na-i]=sa[i]-'0';
for(i=0;i<nb;i++)b[nb-i]=sb[i]-'0';
nc=1;
x=0;
while(nc<=na||nc<=nb)
{
c[nc]=a[nc]+b[nc]+x;
x=c[nc]/10;c[nc]=c[nc]%10;
nc++;
}
c[nc]=x;
while(c[nc]==0&&nc>1)nc--;
for(i=nc;i>=1;i--)printf("%d",c[i]);
return 0;
}
高精度减法
洛谷P2142
模拟竖式减法。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<iomanip>
#include<algorithm>
using namespace std;
char s1[1005],s2[1005],s3[1005];
int a[1005],b[1005],c[1005];
int main()
{
int i,j,la,lb,tmp;
scanf("%s",s1);scanf("%s",s2);
la=strlen(s1);lb=strlen(s2);
if(la<lb||(la==lb&&(strcmp(s1,s2)<0))){
printf("-");
strcpy(s3,s1);strcpy(s1,s2);strcpy(s2,s3);
tmp=la;la=lb;lb=tmp;
}
for(i=0;i<la;i++)a[la-i]=s1[i]-'0';
for(i=0;i<lb;i++)b[lb-i]=s2[i]-'0';
i=1;
while(i<=la&&i<=lb){
if(a[i]<b[i]){a[i]+=10;a[i+1]--;}
c[i]=a[i]-b[i];
i++;
}
while(i<=la){
if(a[i]<0){a[i]+=10;a[i+1]--;}
c[i]=a[i];i++;
}
while(!c[i]&&i>1)i--;
for(j=i;j>=1;j--)printf("%d",c[j]);
return 0;
}
高精乘高精
模拟竖式乘法。
洛谷P1303
#include<iostream>
#include<cstdio>
#include<cstring>
#include<iomanip>
#include<algorithm>
using namespace std;
char s1[2005],s2[2005];
int a[2005],b[2005],c[4005];
int main()
{
int i,j,la,lb,lc,x;
scanf("%s",s1);scanf("%s",s2);
la=strlen(s1);lb=strlen(s2);
for(i=0;i<la;i++)a[la-i]=s1[i]-'0';
for(i=0;i<lb;i++)b[lb-i]=s2[i]-'0';
for(i=1;i<=la;i++){
x=0;
for(j=1;j<=lb;j++){
c[i+j-1]+=a[i]*b[j]+x;
x=c[i+j-1]/10;c[i+j-1]%=10;
}
c[i+lb]+=x;//i+lb!
}
lc=la+lb;
while(!c[lc]&&lc>1)lc--;
for(i=lc;i>=1;i--)printf("%d",c[i]);
return 0;
}
高精除低精
这是很简单的模拟竖式除法
#include<iostream>
#include<cstdio>
#include<cstring>
#include<iomanip>
#include<algorithm>
using namespace std;
char s1[5005];
int a[5005],c[5005];
long long b;
int main()
{
int i,j,la,lc,x=0;
scanf("%s",s1);scanf("%lld",&b);
la=strlen(s1);
for(i=0;i<la;i++)a[i+1]=s1[i]-'0';
for(i=1;i<=la;i++){
a[i]+=x*10;c[i]=a[i]/b;x=a[i]%b;
}
lc=1;
while(!c[lc]&&lc<la)lc++;
for(i=lc;i<=la;i++)printf("%d",c[i]);
return 0;
}
高精除高精
本来以为很高级,最后发现是简单的模拟减法…
随便啦…还是写一写啦…
模拟减法还是有窍门的,一个一个的减会减到天荒地老,只要十倍百倍千倍…地减就会快很多,具体看代码吧。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<iomanip>
#include<algorithm>
using namespace std;
char s1[5005],s2[5005];
int a[10005],b[10005],c[10005],tmp[10005];
int com(int a[],int b[]){//比较大小,a<b -1,a>b 1,a=b 0
int i;
if(a[0]<b[0])return -1;
if(a[0]>b[0])return 1;
for(i=a[0];i>=1;i--)
if(a[i]<b[i])return -1;
else if(a[i]>b[i])return 1;
return 0;
}
void jian(int a[],int b[]){//高精度减法
int flag,i,x;
flag=com(a,b);
if(flag==0){a[0]=1;a[1]=0;return;}
for(i=1;i<=a[0];i++){
x=0;
if(a[i]<b[i]){a[i]+=10;a[i+1]--;}
a[i]-=b[i];
}
while(a[0]>1&&!a[a[0]])a[0]--;
}
void go(int a[],int b[],int c[]){
int i,j,lc,cnt=0;
lc=a[0]-b[0]+1;//因为从第一位开始多了很多多余的计算,所以我们可以豪迈点,直接十倍百倍千倍的减
for(i=lc;i>=1;i--){
memset(tmp,0,sizeof(tmp));
//复制b数组到tmp数组从i开始的地方
for(j=1;j<=b[0];j++)
tmp[j+i-1]=b[j];
tmp[0]=b[0]+i-1;
while(com(a,tmp)>=0){
c[i]++;
jian(a,tmp);cnt++;
}
}
while(lc>1&&!c[lc])lc--;
for(i=lc;i>=1;i--)printf("%d",c[i]);
}
int main()
{
int i,j;
scanf("%s",s1);scanf("%s",s2);
a[0]=strlen(s1);b[0]=strlen(s2);
for(i=0;i<a[0];i++)a[a[0]-i]=s1[i]-'0';
for(i=0;i<b[0];i++)b[b[0]-i]=s2[i]-'0';
go(a,b,c);
return 0;
}
高精度取模
如果是高精除高精就是最后减法过后A数组里剩下的东西。
如果是高精除低精…就…
for(int i=1;i<=len;++i) ans=ans*10%mod+a[i];