C++实现大正整数及其相关运算(长期更新)
/**
只考虑正数[1, +∞);
“-”运算只允许大数减小数;
小端存储;
*/
typedef struct BigInteger0 {
vector<int> v;
BigInteger0(int len) {
v.resize(len);
}
BigInteger0(const vector<int>& nv) {
v.assign(nv.begin(), nv.end());
}
BigInteger0(const char str[]) {
int len = strlen(str);
v.resize(len);
; i < len; ++i)
v[len - - i] = str[i] - ';
}
BigInteger0 operator + (const BigInteger0& obj) const {
BigInteger0 res(v);
int len = max(v.size(), obj.v.size());
res.v.resize(len + );
; i < len; ++i) {
if(i < obj.v.size())res.v[i] += obj.v[i];
) {
res.v[i + ] += res.v[i] / ;
res.v[i] %= ;
}
}
] == )res.v.resize(res.v.size() - );
return res;
}
/*默认v > obj,即v.size() > obj.v.size()*/
BigInteger0 operator - (const BigInteger0& obj) const {
BigInteger0 res(v);
; i < v.size(); ++i) {
if(i < obj.v.size())res.v[i] -= obj.v[i];
) {
res.v[i + ]--;
res.v[i] += ;
}
}
int len = res.v.size();
] == ; --len);
res.v.resize(len);
return res;
}
BigInteger0 operator * (const BigInteger0& obj) const {
BigInteger0 res(v.size() + obj.v.size());
, len1 = v.size(); i < len1; ++i) {
, len2 = obj.v.size(); j < len2; ++j) {
res.v[i + j] += v[i] * obj.v[j];
) {
res.v[i + j + ] += res.v[i + j] / ;
res.v[i + j] %= ;
}
}
}
int len = res.v.size();
] == ; --len);
res.v.resize(len);
return res;
}
BigInteger0 operator / (BigInteger0& obj) const {
");
");
else {
BigInteger0 res("), tmp(v);
while(tmp > obj) {
int lendif = tmp.v.size() - obj.v.size();
BigInteger0 b1 = ten ^ lendif;
BigInteger0 b2 = b1 * obj;
if(tmp < b2) {
b1 = ten ^ (lendif - );
b2 = b1 * obj;
}
while(tmp >= b2) {
tmp = tmp - b2;
res = res + b1;
}
}
return res;
}
}
/**power*/
BigInteger0 operator ^ (int n) const {
BigInteger0 res("), tmp(v);
; n >>= ) {
)res = res * tmp;
tmp = tmp * tmp;
}
return res;
}
BigInteger0 sqrt() {
) / ;
BigInteger0 res(len), ten(");
; i >= ; --i) {
, high = , mid, bit;
) {
mid = (low + high) >> ;
BigInteger0 b = res + (BigInteger0(vector<, mid)) * (ten ^ i));
if(b * b <= *this) {
low = mid;
bit = mid;
}
else high = mid;
}
res.v[i] = bit;
}
] == ; --len);
res.v.resize(len);
return res;
}
bool operator < (const BigInteger0& obj) const {
if(v.size() < obj.v.size())return true;
else if(v.size() > obj.v.size())return false;
else {
;
&& v[i] == obj.v[i]; --i);
)return false;
else return v[i] < obj.v[i];
}
}
bool operator > (const BigInteger0& obj) const {
return obj < *this;
}
bool operator == (const BigInteger0& obj) const {
if(v.size() != obj.v.size())return false;
else {
; i < v.size(); ++i) {
if(v[i] != obj.v[i])return false;
}
return true;
}
}
bool operator >= (const BigInteger0& obj) const {
return *this > obj || *this == obj;
}
bool operator <= (const BigInteger0& obj) const {
return *this < obj || *this == obj;
}
string value() {
string res;
; i >= ; --i)res.push_back(v[i] + ');
)res.push_back(');
return res;
}
} BigInteger;
V2.0
/**
只考虑正数[1, +∞);
“-”运算只允许大数减小数;
小端存储;
*/
typedef struct BigInteger0 {
typedef long long LL;
, BOUND = int(1e9);
vector<LL> v;
void trim() {
;
&& v[i] == ; --i);
)v.resize(i + );
}
string value() {
];
string res;
trim();
)res.push_back(');
else {
; i >= ; --i) {
)sprintf(tmp, "%I64d", v[i]);
else sprintf(tmp, "%09I64d", v[i]);
res += string(tmp);
}
}
return res;
}
bool isEven(){
] % == ;
}
BigInteger0(const vector<LL>& nv): v(nv) {}
BigInteger0(int len) {
v.resize(len);
}
BigInteger0(const char str[]) {
int len = strlen(str);
, cnt = , last = len;
v.resize(cap);
, len - BIT); i >= ; i = max(, i - BIT)) {
LL tmp = ;
for(int j = i; j < last; ++j) {
tmp = tmp * + str[j] - ';
}
v[cnt++] = tmp;
last = i;
)break;
}
trim();
}
BigInteger0 operator + (const BigInteger0& obj) const {
BigInteger0 res(v);
int len = max(v.size(), obj.v.size());
res.v.resize(len + );
; i < len; ++i) {
if(i < obj.v.size())res.v[i] += obj.v[i];
if(res.v[i] >= BOUND) {
res.v[i + ] += res.v[i] / BOUND;
res.v[i] %= BOUND;
}
}
res.trim();
return res;
}
/*默认v >= obj,即v.size() >= obj.v.size()*/
BigInteger0 operator - (const BigInteger0& obj) const {
BigInteger0 res(v);
; i < v.size(); ++i) {
if(i < obj.v.size())res.v[i] -= obj.v[i];
) {
res.v[i + ]--;
res.v[i] += BOUND;
}
}
res.trim();
return res;
}
BigInteger0 operator * (const BigInteger0& obj) const {
BigInteger0 res(v.size() + obj.v.size());
, len1 = v.size(); i < len1; ++i) {
, len2 = obj.v.size(); j < len2; ++j) {
res.v[i + j] += v[i] * obj.v[j];
if(res.v[i + j] >= BOUND) {
res.v[i + j + ] += res.v[i + j] / BOUND;
res.v[i + j] %= BOUND;
}
}
}
res.trim();
return res;
}
BigInteger0 divide(const BigInteger0& obj, BigInteger0& reminder){
BigInteger0 zero(");
if(*this < obj){reminder = *this; return zero;}
else if(*this == obj){reminder = zero; return one;}
else{
BigInteger0 tmp(v), p1 = one, p2 = one, res = zero;
while(tmp >= obj){
p1 = one;
while(p1 * obj <= tmp){
p2 = p1;
p1 = p1 * two;
}
res = res + p2;
tmp = tmp - (obj * p2);
}
reminder = tmp;
return res;
}
}
BigInteger0 operator / (const BigInteger0& obj){
BigInteger0 rem(");
return this->divide(obj, rem);
}
BigInteger0 operator % (const BigInteger0& obj){
BigInteger0 rem(");
this->divide(obj, rem);
return rem;
}
/**power*/
BigInteger0 operator ^ (BigInteger0 n) const {
BigInteger0 zero(");
BigInteger0 res("), tmp(v);
for(; n > zero; n = n / two) {
if(!n.isEven())res = res * tmp;
tmp = tmp * tmp;
}
return res;
}
bool operator < (const BigInteger0& obj) const {
if(v.size() < obj.v.size())return true;
else if(v.size() > obj.v.size())return false;
else {
;
&& v[i] == obj.v[i]; --i);
)return false;
else return v[i] < obj.v[i];
}
}
bool operator > (const BigInteger0& obj) const {
return obj < *this;
}
bool operator == (const BigInteger0& obj) const {
if(v.size() != obj.v.size())return false;
else {
; i < v.size(); ++i) {
if(v[i] != obj.v[i])return false;
}
return true;
}
}
bool operator >= (const BigInteger0& obj) const {
return *this > obj || *this == obj;
}
bool operator <= (const BigInteger0& obj) const {
return *this < obj || *this == obj;
}
BigInteger0& operator = (const BigInteger0& obj) {
v.assign(obj.v.begin(), obj.v.end());
return *this;
}
} BigInteger;
C++实现大正整数及其相关运算(长期更新)的更多相关文章
- js 时间函数 及相关运算大全
js 时间函数 及相关运算大全 var myDate = new Date(); myDate.getYear(); //获取当前年份(2位) myDate.getFullYear(); ...
- Problem B: 大整数的加法运算
Problem B: 大整数的加法运算 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 112 Solved: 57[Submit][Status][W ...
- 大整数加减运算的C语言实现
目录 大整数加减运算的C语言实现 一. 问题提出 二. 代码实现 三. 效果验证 大整数加减运算的C语言实现 标签: 大整数加减 C 一. 问题提出 培训老师给出一个题目:用C语言实现一个大整数计算器 ...
- Python 整数和浮点数运算
和数学运算不同的地方是,Python的整数运算结果仍然是整数,浮点数运算结果仍然是浮点数:1 + 2 # ==> 整数 31.0 + 2.0 # ==> 浮点数 3.0 整数和浮点数混合运 ...
- Python:Base1(数据类型,print语句,变量,定义字符串,raw字符串与多行字符串,Unicode字符串,整数和浮点数运算,布尔类型运算)
1.Python中数据类型: 计算机顾名思义就是可以做数学计算的机器,因此,计算机程序理所当然地可以处理各种数值.但是,计算机能处理的远不止数值,还可以处理文本.图形.音频.视频.网页等各种各样的数据 ...
- Java代码优化(长期更新)
前言 2016年3月修改,结合自己的工作和平时学习的体验重新谈一下为什么要进行代码优化.在修改之前,我的说法是这样的: 就像鲸鱼吃虾米一样,也许吃一个两个虾米对于鲸鱼来说作用不大,但是吃的虾米多了,鲸 ...
- iOS-工作经验+资料分享(长期更新)
在此记录工作中的一些经验和技术资料 长期更新 欢迎各位业内朋友指正.交流技术上的问题 0.苹果开发联盟电话 4006 701855 1.轻易不用使用tableViewController,因为改变他自 ...
- Ubuntu16.04系统美化、常用软件安装等,长期更新
Ubuntu16.04系统美化.常用软件安装等,长期更新 IT之家啊 18-09-0915:00 因为我个人偏向于玩VPS.服务器之类的东西,所以一般我都是用CentOS.不过对于桌面版的Linux, ...
- Spring MVC & Boot & Cloud 技术教程汇总(长期更新)
昨天我们发布了Java成神之路上的知识汇总,今天继续. Java成神之路技术整理(长期更新) 以下是Java技术栈微信公众号发布的关于 Spring/ Spring MVC/ Spring Boot/ ...
随机推荐
- Win10累积更新KB3081424安装失败的解决方法
我在个人的PC上多次更新KB3081424都失败了,后来按照网上给出的删除注册表残留项的方法后更新成功了(删除注册表残留项后我还执行了从微软帐户切换到本地帐户的操作,有兴趣的朋友可以试试不切换能否更新 ...
- Codeforces Round #401 (Div. 2) A,B,C,D,E
A. Shell Game time limit per test 0.5 seconds memory limit per test 256 megabytes input standard inp ...
- [翻译]将智能指针用于C++的类成员
http://stackoverflow.com/questions/15648844/using-smart-pointers-for-class-members Question: I'm hav ...
- 线程池ThreadPoolExecutor里面4种拒绝策略
ThreadPoolExecutor类实现了ExecutorService接口和Executor接口,可以设置线程池corePoolSize,最大线程池大小,AliveTime,拒绝策略等.常用构造方 ...
- redis的使用及方法
一.redis (1).redis是一个key-value存储系统.和Memcached类似,它支持存储的value类型相对更多,包括string(字符串).list(链表).set(集合).zset ...
- Java回顾之JDBC
这篇文章里,我们来讨论一些和JDBC相关的话题. 概述 尽管在实际开发过程中,我们一般使用ORM框架来代替传统的JDBC,例如Hibernate或者iBatis,但JDBC是Java用来实现数据访问的 ...
- Java网络编程学习A轮_03_抓包分析TCP四次挥手
参考资料: http://www.jellythink.com/archives/705 示例代码: https://github.com/gordonklg/study,socket module ...
- Mongo Plugin插件(编辑器PyCharm的Mongo插件安装与使用)
博主接触到MongoDB数据库.用普通的Navicat工具 是不支持的 正准备重新安装一款对应的可视化工具.刚好发现在PyCharm编辑中有连接mongoDB数据的插件 Mongo Plugin 这里 ...
- UVA-10600 ACM Contest and Blackout (次小生成树)
题目大意:给一张无向图,找出最小生成树和次小生成树. 题目分析:模板题...方法就是枚举所有的比最小生成树中两端点之间的最长边还要长的边,用它替换,再取一个最小的值便是次小生成树了. 代码如下: # ...
- Oracle linux安装Oracle 11G
Oracle linux安装Oracle 11G 系统环境 Oracle linux 5.8.19.6(64位) Oracle 11.2.0.1(64位) 本文档详细介绍Oracle ...