【模板】vector + __int128压位高精度大全

/*

包含操作:
快读            rd
快输            print
高精度读入      iird
高精度输出      iiprint
高精度比较      iicmp
高精度加法      pls
高精度减法      sub
高精乘低精      mul
高精乘高精      iimul
高精整除以低精  div
对低精取模      mod
高精度gcd       iigcd
高精度快速幂    iipow
高精度开n次方根 root
// 注:也可以专门写 ++ 和 -- 的函数
// 文末附有高精度整除以高精度 和 负数高精度的修改方法

*/

// 笔记:
// 补齐前导零
// 是否需要x存储进位?
// 如果有x,注意处理最高位剩余的x
// 去除前导零

#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
// I love OI!

typedef __int128 ii;
typedef long long ll;
#define V vector<ii>
#define il inline 
#define pu push_back
#define sz size() // 改代码使用时强转为int
const ii base = 1e18;
const int baselen = 18;
// 下标从0开始
// 若只需加减运算,压位可以适当压多一点,保证两数相加不溢出即可
// 若需乘法运算,压位一般只压一半大小,保证两数相乘不溢出
// 该代码默认没有负数

// 快读
template <typename TT>
void inline rd(TT &x){
    x = 0;
    char c = getchar();
    for ( ; c < '0' || c > '9' ; c = getchar());
    for ( ; c >= '0' && c <= '9' ; c = getchar()) x = (x << 3) + (x << 1) + (c ^ 48);
}

// 快输
template <typename tt>
void il print(tt x){
	if (x > 9) print(x / 10);
	putchar(x % 10 + '0');
}

// 高精度读入
void il iird(V &A){
	A.clear();
	string S;
	cin >> S;
	reverse(S.begin(), S.end());
	for (int i = 0 ; i < (int)S.sz ; i += baselen){
		ii x = 0;
		for (int j = min((int)S.sz, i + baselen) - 1 ; j >= i ; j --)
			x = (x << 3) + (x << 1) + (S[j] ^ 48);
		A.pu(x);
	}
}

// 高精度输出
void il iiprint(V A){
    print(A[A.sz - 1]);
	for (int i = A.sz - 2 ; i >= 0 ; i --) 
		printf("%018lld", (ll)A[i]); // 改变baselen时需要同步修改
}

// 笔记:
// 补齐前导零
// 是否需要x存储进位?
// 如果有x,注意处理最高位剩余的x
// 去除前导零

// 高精度比较:若A < B返回1, A > B返回-1, A = B返回0
int il iicmp(V A, V B){
	if (A.sz != B.sz) return A.sz < B.sz ? 1 : -1;
	for (int i = A.sz - 1 ; i >= 0 ; i --)
		if (A[i] != B[i]) return A[i] < B[i] ? 1 : -1;
	return 0;
}

// 高精度加法
V il pls(V A, V B){
	while(A.sz < B.sz) A.pu(0); // 补齐前导零
	while(A.sz > B.sz) B.pu(0);
	ii x = 0; 
	V C;
	for (int i = 0, len = A.sz ; i < len ; i ++){
		x += A[i] + B[i];
		C.pu(x % base);
		x /= base;
	}
	if (x) C.pu(x); // 处理最高位剩余的x
	return C;
}

// 高精度减法
// 请提前用iicmp判断正负,这里默认A > B
V il sub(V A, V B){
	V C;
	while(B.sz < A.sz) B.pu(0); // 补齐前导零
	for (int i = 0, len = A.sz ; i < len ; i ++){
		if (A[i] < B[i]) {
			A[i + 1] --;
			A[i] += base;
		}
		C.pu(A[i] - B[i]);
	}
	while(C.sz > 1 && !C.back()) C.pop_back(); // 去除前导零
	return C;
}

// 高精乘低精
il V mul(V A, ii B){
	ii x = 0;
	V C;
	for (int i = 0, len = A.sz ; i < len ; i ++){
		x += A[i] * B;
		C.pu(x % base);
		x /= base;
	}
	while(x) C.pu(x % base), x /= base; // 处理最高位剩余的x
	while(C.sz > 1 && !C.back()) C.pop_back(); // 去除前导零
	return C;
}

// 高精乘高精
// 优化:FFT
il V iimul(V A, V B){
	V C;
	while(C.sz < A.sz + B.sz) C.pu(0);
	for (int i = 0, lenA = A.sz ; i < lenA ; i ++){
		ii x = 0;
		for (int j = 0, lenB = B.sz ; j < lenB ; j ++){
			x += C[i + j] + A[i] * B[j]; // 由于下标从0开始,注意这里是i + j
										 // 如果下标从1开始,则为 i + j - 1 
			C[i + j] = x % base;
			x /= base;
		}
		C[i + B.sz] = x; // 处理最高位剩余的x
	}
	while(C.sz > 1 && !C.back()) C.pop_back(); // 去除前导零
	return C;
}

// 高精整除以低精
V il div(V A, ii B){
	V C;
	ii x = 0;
	// 注意除法从最高位开始
	for (int i = A.sz - 1 ; i >= 0 ; i --){
		x = x * base + A[i];
		C.pu(x / B);
		x %= B;
	}
	reverse(C.begin(), C.end()); // 反转
	while(C.sz > 1 && !C.back()) C.pop_back(); // 去除前导零
	return C;
}

// 对低精取模
il ii mod(V A, ii B){
	ii x = 0;
	for (int i = A.sz - 1 ; i >= 0 ; i --)
		x = (x * base + A[i]) % B;
	return x;
}

// 高精度gcd
V il iigcd(V A, V B){
	ii x = 0;
	V C;
	C.pu(1);
	while(!iicmp(B, C)){
		if (iicmp(A, B)) swap(A, B);
		if (iicmp(B, C)) break;
		if (A[0] % 2 == 0 && B[0] % 2 == 0)
			A = div(A, 2), B = div(B, 2), x ++;
		else if (A[0] % 2 == 0) A = div(A, 2);
		else if (B[0] % 2 == 0) B = div(B, 2);
		else A = sub(A, B);
	}
	while(x --) A = mul(A, 2);
	return A;
}

// 高精度快速幂
V il iipow(V A, ii k){
	if (k == 1) return A;
	V C;
	C.pu(1);
	for ( ; k ; k >>= 1, A = iimul(A, A))
		if (k & 1) C = iimul(C, A);
	return C;
}

// 高精度开n次方根(二分法+快速幂)
// 其他做法:基于倍增的牛顿迭代法
// https://www.luogu.com.cn/problem/solution/P2293
V il root(V A, ii k){
	if (k == 1) return A;
	V L, R = A, mid, _1, ans;
	mid.pu(0);
	if (A.empty() || (A.sz == 1 && A[0] == 0)) return mid;
	L.pu(1);
	_1.pu(1);
	while(iicmp(L, R) != -1){
		mid = div(pls(L, R), 2);
		if (mid.sz * k - k + 1 > A.sz) R = sub(mid, _1);
		else 
		if (iicmp(A, iipow(mid, k)) != 1) {
		    if (iicmp(mid, ans) != 1) ans = mid;
		    L = pls(mid, _1);
		} else R = sub(mid, _1);
	}
	return ans;
}


int main(){
	// 注意初始化,如A.pu(0),乘法运算前应先A.pu(1)
    V x, y;
	iird(x); 
	iird(y);
    x = pls(x, y);
	iiprint(x);
}

高精整除以高精,模拟竖式除法并转化为减法

void numcpy(int p[],int q[],int det){//将p数组拷贝到q数组中,从q的det位置开始 
	for(int i=1;i<=p[0];i++)q[i+det-1]=p[i];
	q[0]=p[0]+det-1;//更新q的位数 
}
void chu(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 ;
} 

其它做法:
如每次不是减1次,而是减10的幂次
或转化用二分枚举商,再用乘法判断
https://blog.csdn.net/lybc2019/article/details/103712905

高精度有负数的修改方式

iird中:略
iiprint中:if(f) putchar('-'),f=0;
pls中:

if(f1^f2){
	if(f1) f1^=1, mnu(b,a), f1^=1;
	if(f2) f2^=1, mnu(a,b), f2^=1; //加负数等效于减正数
	return;
}
if(f1&f2){
	f1=f2=0, f^=1, pls(a,b);
	f1=f2=1;
	return;
}

sub中:

if(f1^f2){
	if(f1) f1^=1,f^=1,pls(a,b),f1^=1;
	if(f2) f2^=1,pls(a,b),f2^=1;//减负数等效于加正数
	return;
}
if(f1&f2){
	f1=f2=0,mnu(b,a);
	f1=f2=1;
	return;
}
if(ABScmp(a,b) == 1){ // 注意这里是比较绝对值,|a| < |b|
	f^=1;mnu(b,a);return; // (a - b) = - (b - a)
}

mul中:略
iimul中:if(f1^f2) f^=1;
iidiv中:

if(f1^f2){
	if(f1) f1^=1,f^=1,div(a,b),f1^=1;
	if(f2) f2^=1,f^=1,div(a,b),f2^=1;
	return;
}

iicmp中:略
其它函数的修改:略
参考博客:https://www.luogu.com.cn/blog/AH2002/solution-p2005

利用reserve和resize提高vector的效率

由于需要不断进行push_back,故不用resize
https://blog.csdn.net/qq_37037492/article/details/86568290

memcpy的用法:https://blog.csdn.net/qq_35040828/article/details/71123521

posted @ 2021-09-23 13:34  Randolph、  阅读(161)  评论(0编辑  收藏  举报