二分+高精度——BZOJ1213/Luogu2293 [HNOI2004]高精度开根

http://www.lydsy.com/JudgeOnline/problem.php?id=1213
https://www.luogu.org/problem/show?pid=2293
首先说一下，这又是我们某次膜你赛的题目。。。

---

先上题解：这题只要求一个数非负根取整后的结果，没有小数部分，所以思路很简单
我们可以直接二分答案，二分出某个值然后判断这个数的m次幂与原数判断：

• mid^m==n，答案直接输出
• mid^m < n，l=mid+1
• mid^m > n，r=mid-1

答案是r

    l.l=1;l.c[l.l]=1;r.l=k;
    r.c[k]=10000000;
    kkk.l=1;kkk.c[1]=5;kk.l=1;kk.c[1]=1;
    while(bj(l,r)==1){
        m=gjj(l,r);
        m=gjc(m,kkk);
        m=kd(m);
        rp=check(m,n);
        if(rp==0){r=m;break;}
        if(rp==1)r=gj(m);
        else l=gjj(m,kk);
    }    

其中gjj是加法，gjc是乘法，/2这个操作我有点暴力：先乘5，然后暴力向后退一位（代码中的kd函数）
使用Python的同学下面应该也不用看了。。。（BZOJ肯定过luogu没有py233）

---

没压位自测时发现还是TLE了（评测每点6s），无奈中压了9位
番外：我居然发现我的程序不用快速幂居然比用了快速幂还快。。。
压9位以后本机上是跑过了，交到BZOJ也是轻松无压力AC（总时限50s）
然而Luogu时限没开大T了n发（不能再压了），和管理员大大交涉了没多久开大了时限总算过了。。。
能压位尽量呀啊。。。

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<string>
#define ll long long
using namespace std;
struct bigint{ll c[1301],l;}x,y,p,tmp,jzq,l,r,m,kkk,kk,abc;
char c[10001];
int n,rp;
inline int abcdefg(int x){
    if(x>=10000000)return 0;
    if(x>=1000000)return 1;
    if(x>=100000)return 2;
    if(x>=10000)return 3;
    if(x>=1000)return 4;
    if(x>=100)return 5;
    if(x>=10)return 6;
    return 7;
}
inline void shuchu(bigint a){
    int rp;
    for(int i=a.l;i>0;i--){
        if(i<a.l){rp=abcdefg(a.c[i]);for(int j=1;j<=rp;j++)printf("0");}
        printf("%d",a.c[i]);
    }
}
inline int bj(bigint a,bigint b){
    if(a.l<b.l)return 1;
    if(a.l>b.l)return 0;
    for(int i=a.l;i>0;i--)if(a.c[i]<b.c[i])return 1;
    else if(a.c[i]>b.c[i])return 0;
    return 2;
}
inline bigint gjj(bigint a,bigint b){
    memset(p.c,0,sizeof p.c);
    p.l=max(a.l,b.l);
    for(int i=1;i<=p.l;i++){
        if(i<=a.l)p.c[i]+=a.c[i];
        if(i<=b.l)p.c[i]+=b.c[i];
        if(p.c[i]>99999999){
            if(i==p.l)p.l++;
            p.c[i]-=100000000;
            p.c[i+1]++;
        }
    }
    return p;
}
inline bigint gjc(bigint a,bigint b){
    int i,j;
    memset(tmp.c,0,sizeof tmp.c);
    for(i=1;i<=a.l;i++)
        for(j=1,tmp.l=i-1;j<=b.l;j++)tmp.c[++tmp.l]+=b.c[j]*a.c[i];
    for(i=1;i<=tmp.l;i++)if(tmp.c[i]>99999999){
        if(tmp.c[tmp.l]>99999999)tmp.l++;
        tmp.c[i+1]+=tmp.c[i]/100000000;
        tmp.c[i]%=100000000;
    }
    return tmp;
}
inline bigint mi(bigint a,int b){
    x=a;
    int rprp;
    for(rprp=1;rprp<b;rprp++)x=gjc(x,a);
    return x;
}
inline int check(bigint a,int b){
    abc=mi(a,b);
    int rp=bj(abc,jzq);
    if(rp==2)return 0;
    if(rp==0)return 1;
    return 2;
}
inline bigint gj(bigint p){
    int jsg=1;
    for(int i=1;i<=p.l;i++,jsg++)if(p.c[i]){
        p.c[i]--;
        if(i==p.l&&p.c[i]==0)p.l--;
        break;
    }
    for(int i=1;i<jsg;i++)p.c[i]=99999999;
    return p;
}
inline bigint kd(bigint p){
    for(int i=1;i<=p.l;i++){
        p.c[i]/=10;
        if(i==p.l&&p.c[i]==0)p.l--;
        p.c[i]+=(p.c[i+1]%10)*10000000;
    }
    return p;
}
inline ll e(int rp){
    if(rp==0)return 1;
    if(rp==1)return 10;
    if(rp==2)return 100;
    if(rp==3)return 1000;
    if(rp==4)return 10000;
    if(rp==5)return 100000;
    if(rp==6)return 1000000;
    if(rp==7)return 10000000;
}
int main()
{
    scanf("%d",&n);scanf("%s",&c);
    if(n==1){
        printf("%s",c);
        return 0;
    }
    if(c[0]=='0'){
        printf("0");
        return 0;
    }
    int le=strlen(c);
    jzq.l=1;rp=0;
    for(int i=le-1;i>=0;i--)if(c[i]>='0'&&c[i]<='9'){
        jzq.c[jzq.l]+=(int)(c[i]-'0')*e(rp);
        rp=(rp+1)%8;
        if(rp==0)jzq.l++;
    }
    if(jzq.c[jzq.l]==0)jzq.l--;
    int k=jzq.l/n+2;
    memset(l.c,0,sizeof l.c);
    memset(r.c,0,sizeof r.c);
    memset(m.c,0,sizeof m.c);
    l.l=1;l.c[l.l]=1;r.l=k;
    r.c[k]=10000000;
    kkk.l=1;kkk.c[1]=5;kk.l=1;kk.c[1]=1;
    while(bj(l,r)==1){
        m=gjj(l,r);
        m=gjc(m,kkk);
        m=kd(m);
        rp=check(m,n);
        if(rp==0){r=m;break;}
        if(rp==1)r=gj(m);
        else l=gjj(m,kk);
    }    
    if(bj(mi(r,n),jzq)==0)r=gj(r);
    shuchu(r);
    return 0;
}