高精度（加 减 乘 除 快速幂）模板

高精度加法：

#include<bits/stdc++.h>

using namespace std ; 
string a , b ; 

vector<int> add(vector<int> &A , vector<int> &B )
{
	int t = 0 ; 
	vector<int> c ; 
	for(int i = 0 ; i < A.size() || i < B.size() ; i ++ )
	{
		if(i < A.size()) t += A[i] ; 
		if(i < B.size()) t += B[i] ; 
		c.push_back(t % 10 ) ; 
		t /= 10 ; 
	}
	if(t) c.push_back(t) ; 
	return c ; 
}
int main()
{
	vector<int> A , B ; 
	
	cin >> a >> b ; 
	for(int i = a.size() - 1 ; i >= 0 ; i -- ) A.push_back(a[i]-'0');
	for(int i = b.size() - 1 ; i >= 0 ; i -- ) B.push_back(b[i]-'0'); 
	vector<int> C = add(A,B) ; 
	for(int i = C.size() - 1 ; i >= 0; i -- ) cout << C[i] ;
	 
	return 0 ; 
}

高精度减法

#include<bits/stdc++.h>

using namespace std ; 
string a , b ; 

bool cmd(vector<int> &A , vector<int> &B)
{
	if(A.size() != B.size()) return A.size() > B.size() ; 
	else
	{
		for(int i = A.size() - 1 ; i >= 0 ;  i -- ) 
		{
			if(A[i] != B[i]) return A[i] > B[i] ; 
		}
	}
	return true ; 
}
vector<int> sub(vector<int> &A , vector<int> &B )
{
	int t = 0 ; 
	vector<int> C ; 
	for(int i = 0 ; i < A.size() ; i ++ )
	{
		t += A[i] ; 
		if(i < B.size()) t -= B[i] ; 
		C.push_back((t + 10)%10) ;
		if(t >= 0 ) t = 0 ;  
		else t = -1 ; 
	}
	while(C.size() > 1 && C.back() == 0 ) C.pop_back() ; 
	return C ; 
}
int main()
{
	vector<int> A , B ; 
	
	cin >> a >> b ; 
	for(int i = a.size() - 1 ; i >= 0 ; i -- ) A.push_back(a[i]-'0');
	for(int i = b.size() - 1 ; i >= 0 ; i -- ) B.push_back(b[i]-'0'); 
	if(cmd(A,B))
	{
		vector<int> C  = sub(A,B);
		for(int i = C.size() - 1 ; i >= 0 ; i -- ) cout << C[i] ; 
	}
	else 
	{
		vector<int> C = sub(B,A);
		cout <<"-";
		for(int i = C.size() - 1 ; i >= 0 ; i -- ) cout << C[i] ; 
	} 
	 
	return 0 ; 
}

高精度减法模板中 cmd函数是判断A是否大于B

高精度乘法

#include<bits/stdc++.h>

using namespace std ; 
const int N = 1e5 +5 ; 
string a , b ; 

vector<int> mul(const vector<int> &A, const vector<int> &B)
{
	vector<int> t(N); 
	for(int i = 0 ; i < A.size() ; i ++ )
	{
		for(int j = 0 ; j < B.size() ;  j ++ )
		{
			t[i + j ] += A[i] * B[j] ; 
			t[i + j + 1 ] += t[i + j ] /10 ; 
			t[i + j ] %= 10 ; 
		}
	}
	while(t.size() > 1 && t.back() == 0 ) t.pop_back() ; 
	return t ; 
} 
int main()
{
	
	cin >> a >> b ; 
	vector<int> A , B ; 
	
	for(int i = a.size() - 1 ; i >= 0 ; i -- ) A.push_back(a[i] - '0' ) ; 
	for(int i = b.size() - 1 ; i >= 0 ; i -- ) B.push_back(b[i] - '0' );
	
	vector<int> C = mul(A,B) ; 
	
	for(int i = C.size() - 1 ; i >= 0 ; i -- ) cout << C[i] ; 
	return 0 ; 
}

高精度除法

#include<bits/stdc++.h>

using namespace std ; 
string a ; 
int b ; 
int t = 0 ; 
vector<int> mul(vector<int> &A ,int b )
{
	
	vector<int> C ; 
	for(int i = A.size() - 1  ; i >= 0 ; i -- )
	{
		t = t *10 + A[i] ; 
		C.push_back(t / b) ; 
		t %= b ; ;
	}
	reverse(C.begin() , C.end()) ; 
	while(C.size() > 1 && C.back() == 0 ) C.pop_back() ; 
	return C ; 
}
int main()
{
	vector<int> A ; 
	
	cin>> a >> b  ; 
	for(int i = a.size() - 1 ; i >= 0 ; i -- ) A.push_back(a[i]-'0');
	
	vector<int> C = mul(A,b) ; 
	
	for(int i = C.size() - 1 ; i >= 0 ; i -- ) cout << C[i] ; 
	cout << endl <<t <<endl ; 
	return 0 ; 
}

高精度快速幂（洛谷p1045）

#include<bits/stdc++.h>

using namespace std ; 
const int N = 500 ; 
int p ; 
vector<int> ans(N) , a(N) ; 
vector<int> mul(vector<int> &A , vector<int> &B)
{
	vector<int> t(N*2) ; 
	
	for(int i = 0 ; i < N ; i ++ )
		for(int j = 0 ; j < N ; j ++ )
		{
			t[i + j ] += A[i] * B[j] ; 
			t[i + j + 1 ] += t[i + j] /10 ; 
			t[i + j ] %= 10 ; 
		}
		return t ; 
}
void  quickpow(int p )
{
	ans[0] = 1 , a[0] = 2 ; 
	
	while(p)
	{
		if(p & 1)
		ans = mul(ans,a) ; 
		a = mul(a,a) ; 
		p >>=1 ; 
	}
	ans[0] -- ; 
}
int main()
{
	cin >> p ; 
	cout<<ceil(p*log10(2))<<endl;
	quickpow(p); 
	int cnt = 0 ; 
	for(int i = N - 1  ; i >= 0  ; i -- ) 
	{
		cout << ans[i] ;
		cnt ++  ; 
		if(cnt == 50 ) 
		{
			cnt = 0 ; 
			cout << endl ; 
		}
	} 
	return 0 ; 
}

题目要求只需要输出后500位数，所以我们只需要初始ans和a的大小为500即可，但t的大小要弄成2*N，因为在乘法函数里面 0<=i+j <=2*N防止越界

这里还有个stl函数log10和ceil