大数开根号

最新推荐文章于 2021-10-19 18:52:17 发布

Peson_Du

最新推荐文章于 2021-10-19 18:52:17 发布

阅读量1.7k

点赞数

分类专栏：数学+数论文章标签：精度

本文链接：https://blog.csdn.net/Peson_Du/article/details/98028728

版权

这篇博客探讨了如何在处理大数时实现高精度的开方运算，指出所采用的方法在精度方面表现出色。

摘要由CSDN通过智能技术生成

#include<bits/stdc++.h>
typedef long long ll;
const int K=3,N=40003/K,M=pow(10,K),p=998244353,g=3,gi=332748118;
int m,i,l,r,L,lim,re[N],f[11],a[N],b[N],inv;
char s[10003];
struct NUM{
    int t,a[N];
}x,A;
inline int pls(int x,int y){x+=y;return x>=p?x-p:x;}
inline int dec(int x,int y){x-=y;return x<0?x+p:x;}
inline int mul(int x,int y){return (ll)x*y%p;}
inline int pow(int x,int y){
    int ans=1;
    for (;y;y>>=1,x=mul(x,x))
        if (y&1) ans=mul(ans,x);
    return ans;
}
inline void ntt(int *A,int opt){
    for (int i=0;i<lim;i++)
        if (i<re[i]) A[i]^=A[re[i]],A[re[i]]^=A[i],A[i]^=A[re[i]];
    int pp=p-1>>1;
    for (int mid=1;mid<lim;mid<<=1,pp>>=1){
        int wn=pow((opt==1)?g:gi,pp);
        for (int R=mid<<1,j=0;j<lim;j+=R){
            int w=1;