bzoj 1213 //1213: [HNOI2004]高精度开根二分+高精度+结构体+运算符重载

原创于 2019-11-05 17:23:53 发布 · 375 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

跟着大佬学算法

本文深入探讨了高精度开根算法的实现细节，包括二分法、高精度加减乘除运算、结构体和运算符重载等关键概念。通过具体题目实例，展示了如何在代码中有效应用这些技术，解决高精度数学问题。

bzoj 1213 //1213: [HNOI2004]高精度开根 二分+高精度+结构体+运算符重载

bzoj 1213 //1213: [HNOI2004]高精度开根 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1213

//在线测评地址https://www.luogu.org/problem/P2293

更多题解，详见https://blog.csdn.net/mrcrack/article/details/90228694BZOJ刷题记录

高精度运算中，多用加减，少用乘除，若一定要用乘除，那就尽可能多用乘，少用除，否则出处处超时。2019-11-6 21:36

没有重载运算符的代码，编得凌乱，且得分不易。

//1213: [HNOI2004]高精度开根
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1213
//以高精度算法经验，预感此题会超越100行代码
//估计8位压位是少不了了。考虑相乘，数据可能会溢出，采用4位压位。
//估计需编写，比较大小的函数，高精度乘高精度，高精度除低精度，高精度加高精度。
//编写高精度算法，很重要的一条，每编写一个函数都要测试一遍。
//样例通过，bzoj里提交，
//Runtime_Error   1016 kb   16 ms   C++/Edit   2384 B   2019-11-5
//在https://www.luogu.org/problem/P2293提交20分。20分代码如下，尽力了。
#include <stdio.h>
#include <string.h>
#define BASE 10000
#define BIT 4
#define maxn 10020
int a[maxn/BIT],b[maxn/BIT],c[maxn/BIT],m,left[maxn],right[maxn],mid[maxn],left1[maxn],one[10];
char s[maxn];
int cmp(int *x,int *y){//x>y 1;x<y -1;x==y 0
   int i;
   if(x[0]>y[0])return 1;
   if(x[0]<y[0])return -1;
   for(i=x[0];i>=1;i--)
       if(x[i]>y[i])return 1;
       else if(x[i]<y[i])return -1;
   return 0;
}
void read(int *x){//读取，需要多编多练
   int len,i,tmp;
   char c;
   scanf("%s",s+1);
   len=strlen(s+1);
   for(i=1;i<=len/2;i++)c=s[i],s[i]=s[len-i+1],s[len-i+1]=c;//逆序
   x[0]=(len-1)/BIT+1;
   for(i=1;i<=len;i++){//123456789
       if((i-1)%BIT==0)x[(i-1)/BIT+1]=0,tmp=1;
       x[(i-1)/BIT+1]+=(s[i]-'0')*tmp;
       tmp*=10;
   }
}
void print(int *x){
   int i;
   printf("%d",x[x[0]]);
   for(i=x[0]-1;i>=1;i--)printf("%04d",x[i]);
}
void div(int *x,int *y){//y=x/2
   int i;
   x[x[0]+1]=0;
   for(i=x[0];i>=1;i--)y[i]=(x[i+1]%2*BASE+x[i])/2;
   i=x[0];
   while(i&&!y[i])i--;//保证最后结果i>=1
   y[0]=i;
}
void add(int *x,int *y,int *z){//z=x+y
   int i;
   z[0]=x[0]>y[0]?x[0]:y[0];
   for(i=1;i<=z[0];i++){
       z[i]+=x[i]+y[i];//此处错写成z[i]=x[i]+y[i];
       z[i+1]=z[i]/BASE;//此处错写成z[i+1]+=z[i]/BASE;
       z[i]%=BASE;//此处错写成z[i]/=BASE;
   }
   if(z[i])z[0]=i;
}
void mul(int *x,int *y,int *z){//z=x*y
   int i,j;
   z[0]=x[0]+y[0];
   for(i=1;i<=x[0];i++)
       for(j=1;j<=y[0];j++){
           z[i+j-1]+=x[i]*y[j];
           z[i+j-1+1]+=z[i+j-1]/BASE;//此处错写成z[i+j-1+1]=z[i+j-1]/BASE;
           z[i+j-1]%=BASE;
       }
   i=z[0];
   while(i&&!z[i])i--;
   z[0]=i;
}
int judge(int *x){
   int f[maxn/BIT],g[maxn/BIT],i;
   memset(f,0,sizeof(f)),f[0]=1,f[1]=1;
   for(i=1;i<=m;i++){
       memset(g,0,sizeof(g));
       mul(x,f,g);
       memcpy(f,g,sizeof(g));
   }
   if(cmp(g,a)>=0)return 1;
   else return 0;
}
int main(){
   scanf("%d",&m);
   read(a),left[0]=1,left[1]=0,memcpy(right,a,sizeof(a)),one[0]=1,one[1]=1;
   while(cmp(left1,right)==-1){//left+1<right
       memset(b,0,sizeof(b));
       add(left,right,b);
       div(b,mid);
       if(judge(mid))memcpy(right,mid,sizeof(mid));
       else memcpy(left,mid,sizeof(mid));
       memset(left1,0,sizeof(left1));
       add(left,one,left1);
   }
   print(right);
}

AC代码如下

Accepted

1072 kb

21752 ms

C++/Edit

3010 B

//此文https://www.cnblogs.com/lxyyyy/p/10889765.html代码，与本人比较接近，进行对照研究，
//高精度算法，编写过程中种种的不便，决定还是要重载一些运算符
//样例通过，提交AC.2019-11-6 21:25
//看懂是一回事，编写又是另一回事，struct省去了许多初始化的麻烦
#include <cstdio>
#include <cstring>
#define N 10100
#define BIT 4
#define BASE 10000
using namespace std;
int m;
struct num{
    int a[N*2/BIT];//因涉及平方，故*2少不了
    num(){
        memset(a,0,sizeof(a));
    }
    void read(){
        char s[N],c;
        int len,i,tmp;
        scanf("%s",s+1);
        len=strlen(s+1);
        for(i=1;i<=len/2;i++)c=s[i],s[i]=s[len-i+1],s[len-i+1]=c;//逆序
        a[0]=(len-1)/BIT+1;
        for(i=1;i<=len;i++){
            if((i-1)%BIT==0)a[(i-1)/BIT+1]=0,tmp=1;
            a[(i-1)/BIT+1]+=(s[i]-'0')*tmp,tmp*=10;
        }
    }
    void print(){
        int i;
        printf("%d",a[a[0]]);
        for(i=a[0]-1;i>=1;i--)printf("%0*d",BIT,a[i]);
    }
}one,two,ten,l,r,mid,b;
num operator +(const num &p,const num &q){
    num c;
    int i;
    c.a[0]=p.a[0]>q.a[0]?p.a[0]:q.a[0];
    for(i=1;i<=c.a[0];i++)c.a[i]=p.a[i]+q.a[i];
    for(i=1;i<=c.a[0];i++)c.a[i+1]+=c.a[i]/BASE,c.a[i]%=BASE;
    if(c.a[c.a[0]+1])c.a[0]++;
    return c;
}
num operator *(const num &p,const num &q){
    num c;
    int i,j;
    c.a[0]=p.a[0]+q.a[0]-1;
    for(i=1;i<=p.a[0];i++)
        for(j=1;j<=q.a[0];j++)
            c.a[i+j-1]+=p.a[i]*q.a[j];//此处错写成c.a[i+j-1]=p.a[i]*q.a[i];
    for(i=1;i<=c.a[0];i++)c.a[i+1]+=c.a[i]/BASE,c.a[i]%=BASE;
    if(c.a[c.a[0]+1])c.a[0]++;//此处错写成if(c.a[r.a[0]+1])c.a[0]++;
    return c;//漏了此句
}
num operator /(const num &p,int q){
   num c;
   int y=0,i;
   c.a[0]=p.a[0];
   for(i=p.a[0];i>=1;i--){
       y*=BASE,y+=p.a[i];
       c.a[i]=y/q,y%=q;
   }
   while(c.a[0]&&!c.a[c.a[0]])c.a[0]--;//去除前导０
   return c;
}
bool operator <(const num &p,const num &q){
   int i;
   if(p.a[0]!=q.a[0])return p.a[0]<q.a[0];
   for(i=p.a[0];i>=1;i--)
       if(p.a[i]!=q.a[i])return p.a[i]<q.a[i];
   return 0;
}
bool operator ==(const num &p,const num &q){
   int i;
   if(p.a[0]!=q.a[0])return 0;
   for(i=p.a[0];i>=1;i--)
       if(p.a[i]!=q.a[i])return 0;
   return 1;
}
num qpow(num p,int n){//快速幂
   num c=one;
   while(n){
       if(n&1)c=c*p;
       p=p*p,n>>=1;
   }
   return c;
}
int main(){
    one.a[0]=1,one.a[1]=1,two.a[0]=1,two.a[1]=2,ten.a[0]=1,ten.a[1]=10;
    scanf("%d",&m);
    b.read();
    if(m==1){
       b.print();
       return 0;
    }
    r=one;//漏了此句
    while(qpow(r,m)<b)l=r,r=r*two;//此处错写成while(qpow(r,2)<b)l=r,r=r*two;
    b=b*qpow(ten,m),l=l*ten,r=r*ten;//因精度问题，将整体扩大10倍,之后，再将整体缩小为1/10，变回原来。//此处错写成b=b*ten,l=l*ten,r=r*ten;
    while(l+one<r){
       mid=(l+r)/2;
       if(qpow(mid,m)<b)l=mid;
       else r=mid;
    }
    if(qpow(r,m)==b)r=r/10,r.print();//此处错写成b=b/10,l=l/10,r=r/10;
    else l=l/10,l.print();
    return 0;
}