POJ 1503 大整数
之前做的大整数,都是一位一位操作。
优化方案:压缩方案。
模板: + - * 操作符重载
#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<string>
#include<vector>
using namespace std; struct bigint{ // only positive number;
static const int BASE=;
static const int WIDTH=;
vector<int> s;
//value
bigint(long long num=){ *this = num;}
bigint operator = (long long num){
s.clear();
do{
s.push_back(num%BASE);
num/=BASE;
}while(num>);
return *this;
}
bigint operator = (const string& str){
s.clear() ;
int x,len = (str.length()-)/WIDTH + ;
for(int i=;i<len;i++){
int end = str.length() - i*WIDTH;
int start = max(,end - WIDTH);
sscanf(str.substr(start,end-start).c_str(),"%d",&x);
s.push_back(x);
}
return *this;
}
//input&output
friend ostream& operator << (ostream &out, const bigint& x){
out << x.s.back();
for(int i=x.s.size()-;i>=;i--){
char buf[];
sprintf(buf,"%08d",x.s[i]);
for(int j=;j<strlen(buf);j++) out << buf[j];
}
return out;
}
friend istream& operator >>(istream &in, bigint& x){
string s;
if(!(in>>s)) return in;
x=s;
return in;
}
//compare
bool operator < (const bigint& b) const {
if(s.size()!=b.s.size()) return s.size() < b.s.size();
for(int i=s.size()-;i>=;i++) if(s[i]!=b.s[i]) return s[i] < b.s[i];
return false;//equal
}
bool operator > (const bigint& b) const {return b < *this;}
bool operator <= (const bigint& b) const {return !(b < *this);}
bool operator >= (const bigint& b) const {return !(*this < b);}
bool operator != (const bigint& b) const {return b < *this || *this < b;}
bool operator == (const bigint& b) const {return !(b < *this) && !(*this < b);}
//calculate
bigint operator +(const bigint& b) const {
bigint c;
c.s.clear();
for(int i=,g=;;i++){
if(g== && i>=s.size() && i>=b.s.size()) break;
int x=g;
if(i<s.size()) x+=s[i];
if(i<b.s.size()) x+=b.s[i];
c.s.push_back(x%BASE);
g = x/BASE;
}
return c;
}
bigint operator +=(const bigint& b){
*this = *this + b;
return *this;
}
bigint operator -(const bigint& b) const {
bigint c;
c.s.clear();
for(int i=,g=;;i++){
if(g== && i>=s.size() && i>=b.s.size()) break;
int x=g;
if(i<s.size()) x+=s[i];
if(i<b.s.size()) x-=b.s[i];
x+=BASE;
c.s.push_back(x%BASE);
g = x/BASE - ;
}
return c;
}
bigint operator * (const bigint& b) const {
bigint c;
c.s.clear();
bigint g=;
for(int i=;;i++){
if(g.s.size()== && i>=s.size()+b.s.size()-) break;
bigint x;
x.s.clear() ;
for(int j=;j<g.s.size();j++) x.s.push_back(g.s[j]);
if(i<s.size()+b.s.size()-){
for(int j = max( , i-(int)s.size()+);j<=min(i,(int)b.s.size()-);j++){
bigint t = (long long)b.s[j]*s[i-j];
x += t;
}
}
c.s.push_back(x.s[]);
g.s.clear();
if(x.s.size()>) for(int j=;j<x.s.size();j++) g.s.push_back(x.s[j]);
}
return c;
}
}; int main()
{
//freopen("in.txt","r",stdin);
bigint a;
bigint b = ;
bigint sum = ; while(cin>>a) {
if(a==)
break;
sum +=a;
} cout<< sum <<endl; return ;
}
POJ 1503 大整数的更多相关文章
- POJ 1811 大整数素数判断 Miller_Rabin
#include <cstdio> #include <cstring> #include <cmath> #include <ctime> #incl ...
- POJ C++程序设计 编程题#1 大整数的加减乘除
编程题#4:大整数的加减乘除 来源: POJ (Coursera声明:在POJ上完成的习题将不会计入Coursera的最后成绩.) 注意: 总时间限制: 1000ms 内存限制: 65536kB 描述 ...
- POJ 1001 解题报告 高精度大整数乘法模版
题目是POJ1001 Exponentiation 虽然是小数的幂 最终还是转化为大整数的乘法 这道题要考虑的边界情况比较多 做这道题的时候,我分析了 网上的两个解题报告,发现都有错误,说明OJ对于 ...
- POJ 2506 Tiling(递推+大整数加法)
http://poj.org/problem?id=2506 题意: 思路:递推.a[i]=a[i-1]+2*a[i-2]. 计算的时候是大整数加法.错了好久,忘记考虑1了...晕倒. #includ ...
- POJ 1811 Prime Test (Pollard rho 大整数分解)
题意:给出一个N,若N为素数,输出Prime.若为合数,输出最小的素因子.思路:Pollard rho大整数分解,模板题 #include <iostream> #include < ...
- 2981:大整数加法-poj
2981:大整数加法 总时间限制: 1000ms 内存限制: 65536kB 描述 求两个不超过200位的非负整数的和. 输入 有两行,每行是一个不超过200位的非负整数,可能有多余的前导0. 输 ...
- POJ 1503
http://poj.org/problem?id=1503 对于这个题我也是醉了,因为最开始是有学长和我们说过这个题目的,我以为我记得题目是什么意思,也就没看题目,结果按案例去理解题意,结果WA了一 ...
- poj2389-Bull Math(大整数乘法)
一,题意: 大整数乘法模板题二,思路: 1,模拟乘法(注意"逢十进一") 2,倒序输出(注意首位0不输出) 三,步骤: 如:555 x 35 = 19425 5 5 5 5 5 ...
- AC日记——大整数的因子 openjudge 1.6 13
13:大整数的因子 总时间限制: 1000ms 内存限制: 65536kB 描述 已知正整数k满足2<=k<=9,现给出长度最大为30位的十进制非负整数c,求所有能整除c的k. 输入 ...
随机推荐
- Unable to connect to data source (DSN: shangjihuiclient; Network Address: ; Port Number: 53397). Cannot connect to TimesTen Server. Verify that the TimesTen Server is running or verify that your TCP PORT is set correctly.
安装完毕TimesTen后,在客户端连接时一直报错. 解决办法是:重启服务器,在连接没在报这个错误.
- Yii 自带的分页实例
yii自带的分页很好用,简单的几行代码就能把分页搞出来,唯一恼火的是只能写在controller中,所以有时候controller中的方法有点臃肿.废话少说,上代码上图. 一.代码实例: 1.控制器中 ...
- IDEA 导入cordova3.5工程目录注意事项
IDEA 导入cordova3.5工程目录注意事项 1 eclipse很不稳定,有很多小问题.平时我自己用idea,但是当用cordova3.5创建好工程目录是,用eclipse导入时没有问题的.但是 ...
- SpringBoot | 第三十章:Spring-data-jpa的集成和使用
前言 在前面的第九章:Mybatis-plus的集成和使用章节中,介绍了使用ORM框架mybatis-plus进行数据库的访问.今天,我们来简单学习下如何使用spring-data-jpa进行数据库的 ...
- 6、Angular Route 路由
1.没有嵌套路由 类似 ui-route 上述的html用红字标记的是必须导入的.这是因为:路由并不在ng2中,需要我们额外引入,另外我们需要设置base href,这是个什么东西呢?相当于我们后续所 ...
- Scrapy框架之基于RedisSpider实现的分布式爬虫
需求:爬取的是基于文字的网易新闻数据(国内.国际.军事.航空). 基于Scrapy框架代码实现数据爬取后,再将当前项目修改为基于RedisSpider的分布式爬虫形式. 一.基于Scrapy框架数据爬 ...
- css居中那些事
一.css垂直居中 1.line-height(适用于单行文本居中) eg: html:<p class="wordp">123</p>- css: .w ...
- Vue中使用eslint
.eslintrc.js module.exports = { root: true, parser: 'babel-eslint', "env": { "browser ...
- iOS中使用RNCryptor对资源文件加密
原文:http://blog.csdn.net/chenpolu/article/details/46277587 RNCryptor源码https://github.com/RNCryptor/RN ...
- 64位MATLAB和C混合编程以及联合调试
[环境说明] 电脑操作系统:win 7 旗舰版,64位 MATLAB版本:R2013b VS版本:Microsoft Visual Studio 2010 [操作说明] 1.在原来工程的基础上添加下列 ...