MyException - 我的异常网
当前位置:我的异常网» C++ » (维护)高精模板(余除法)

(维护)高精模板(余除法)

www.MyException.Cn  网友分享于:2013-08-22  浏览:0次
(维护)高精模板(缺除法)
#include<cstdio>
#include<cstring>
#include<algorithm>
typedef long long B_INT;
const char p_c[]="%08lld";
const char i_c[]="%lld";
struct Bigint
{
/*
基本类型(char,int,float,double等)的静态成员可以在类里面声明并初始化,
非基本类(char[],string ,自定义等)必须在类里面声明,类外初始化。
*/
	static const B_INT p=8;//压p位,int压4位加乘法可能导致溢出
	static const B_INT base=100000000;//压p位,最大数为10^p-1
	static const int maxl=3000;
	B_INT a[maxl];//a[0]->len,a[i]->从后往前数第i个p位
	Bigint()
	{
		memset(a,0,sizeof(a));
	}
//	Bigint(char *str)
//	{
//		memset(a,0,sizeof(a));
//		B_INT k=0,p=1;
//		char *str1=str+strlen(str)-1;
//		while(str1>=str)
//		{
//			k=k+p*(*str1-48);
//			if(p==base)
//			{
//				a[++a[0]]=k%base;
//				k/=base;
//				p=1;
//			}
//			str1--;
//			p*=10;
//		}
//		a[++a[0]]=k;
//	}
	Bigint(const Bigint& b)
	{
		memcpy(a,b.a,sizeof(a));
	}
	Bigint& operator=(const Bigint& b)
	{
		memcpy(a,b.a,sizeof(a));
	}
	Bigint& operator=(char *str)
	{
		memset(a,0,sizeof(a));
		B_INT k=0,p=1;
		char *str1=str+strlen(str)-1;
		while(str1>=str)
		{
			k=k+p*(*str1-48);
			if(p==base)
			{
				a[++a[0]]=k%base;
				k/=base;
				p=1;
			}
			str1--;
			p*=10;
		}
		a[++a[0]]=k;
		return *this;
	}
	Bigint operator+(const Bigint &b) const
	{
		Bigint c;
		B_INT i;
		c.a[0]=std::max(a[0],b.a[0]);
		for(i=1;i<=c.a[0];i++)
		{
			c.a[i]+=a[i]+b.a[i];
			c.a[i+1]=c.a[i]/base;
			c.a[i]%=base;
		}
		if(c.a[c.a[0]+1]>0)
			c.a[0]++;
		return c;
	}
	Bigint operator*(const Bigint &b) const
	{
		Bigint c;
		B_INT i,j;
		for(i=1;i<=a[0];i++)
			for(j=1;j<=b.a[0];j++)
				c.a[i+j-1]+=a[i]*b.a[j];
		c.a[0]=a[0]+b.a[0]-1;
		for(i=1;i<=c.a[0];i++)
		{
			c.a[i+1]+=c.a[i]/base;
			c.a[i]%=base;
		}
		if(c.a[c.a[0]+1]>0)
			c.a[0]++;
		return c;
	}
	Bigint operator-(const Bigint &b) const//要求保证减数小于被减数 
	{
		Bigint c;
		B_INT i;
		c.a[0]=std::max(a[0],b.a[0]);
		for(i=1;i<=c.a[0];i++)
			c.a[i]=a[i]-b.a[i];
		for(i=1;i<=c.a[0];i++)
			if(c.a[i]<0)
			{
				c.a[i]+=base;
				c.a[i+1]--;
			}
		while(c.a[c.a[0]]==0&&c.a[0]>0)
			c.a[0]--;
		return c;
	}
	Bigint& operator+=(const Bigint &b)
	{
		*this=*this+b;
		return *this;
	}
	Bigint& operator-=(const Bigint &b)
	{
		*this=*this-b;
		return *this;
	}
	Bigint& operator*=(const Bigint &b)
	{
		*this=(*this)*b;
		return *this;
	}
	bool operator<(const Bigint &b) const
	{
		if(a[0]!=b.a[0])
			return a[0]<b.a[0];
		for(B_INT i=a[0];i>0;i--)
			if(a[i]!=b.a[i])
				return a[i]<b.a[i];
		return false;//相等
	}
	/*
	非静态成员函数后面加const(加到非成员函数或静态成员后面会产生编译错误),
	表示成员函数隐含传入的this指针为 const指针,
	决定了在该成员函数中,
	任意修改它所在的类的成员的操作都是不允许的
	(因为隐含了对this指针的const引用);
	唯一的例外是对于 mutable修饰的成员。
	加了const的成员函数可以被非const对象和const对象调用,
	但不加const的成员函数只能被非const对象调用。
	下方b是const,const函数不能修改其数据成员
	*/
	bool operator>(const Bigint &b) const
	{
		return b<*this;
		/*
		if(a[0]!=b.a[0])
			return a[0]>b.a[0];
		for(int i=a[0];i>0;i--)
			if(a[i]!=b.a[i])
				return a[i]>b.a[i];
		return false;//相等
		*/
	}
	bool operator<=(const Bigint &b) const
	{
		return !(b<*this);
		/*
		if(a[0]!=b.a[0])
			return a[0]>b.a[0];
		for(int i=a[0];i>0;i--)
			if(a[i]!=b.a[i])
				return a[i]>b.a[i];
		return true;//相等
		*/
	}
	bool operator>=(const Bigint &b) const
	{
		return !(*this<b);
		/*
		if(a[0]!=b.a[0])
			return a[0]>b.a[0];
		for(int i=a[0];i>0;i--)
			if(a[i]!=b.a[i])
				return a[i]>b.a[i];
		return true;//相等
		*/
	}
	bool operator==(const Bigint &b) const
	{
		if(a[0]!=b.a[0])
			return false;
		for(B_INT i=a[0];i>0;i--)
			if(a[i]!=b.a[i])
				return false;
		return true;
	}
	bool operator!=(const Bigint &b) const
	{
		return !(*this==b);
	}
	void print()
	{
		printf(i_c,a[a[0]]);
		for(B_INT i=a[0]-1;i>0;i--)
			printf(p_c,a[i]);
	}
}x[3];
char str[10010];
int main()
{
	scanf("%s",str);
	x[0]=str;
	scanf("%s",str);
	x[1]=str;
//	if(x[0]<x[1])
//	{
//		printf("-");
//		x[2]=x[0];
//		x[0]=x[1];
//		x[1]=x[2];
//	}
//	(x[0]-x[1]).print();
//	(x[0]*x[1]).print();
	//printf("%d",x[0]>x[1]);
	(x[1]=x[1]).print();
//	puts("");
//	x[0].print();
//	puts("");
//	x[1].print();
	return 0;
}

  

文章评论

编程语言是女人
编程语言是女人
聊聊HTTPS和SSL/TLS协议
聊聊HTTPS和SSL/TLS协议
十大编程算法助程序员走上高手之路
十大编程算法助程序员走上高手之路
当下全球最炙手可热的八位少年创业者
当下全球最炙手可热的八位少年创业者
10个帮程序员减压放松的网站
10个帮程序员减压放松的网站
程序员最害怕的5件事 你中招了吗?
程序员最害怕的5件事 你中招了吗?
鲜为人知的编程真相
鲜为人知的编程真相
如何区分一个程序员是“老手“还是“新手“?
如何区分一个程序员是“老手“还是“新手“?
一个程序员的时间管理
一个程序员的时间管理
程序员应该关注的一些事儿
程序员应该关注的一些事儿
看13位CEO、创始人和高管如何提高工作效率
看13位CEO、创始人和高管如何提高工作效率
5款最佳正则表达式编辑调试器
5款最佳正则表达式编辑调试器
初级 vs 高级开发者 哪个性价比更高?
初级 vs 高级开发者 哪个性价比更高?
旅行,写作,编程
旅行,写作,编程
每天工作4小时的程序员
每天工作4小时的程序员
亲爱的项目经理,我恨你
亲爱的项目经理,我恨你
Java 与 .NET 的平台发展之争
Java 与 .NET 的平台发展之争
总结2014中国互联网十大段子
总结2014中国互联网十大段子
要嫁就嫁程序猿—钱多话少死的早
要嫁就嫁程序猿—钱多话少死的早
Web开发者需具备的8个好习惯
Web开发者需具备的8个好习惯
程序员必看的十大电影
程序员必看的十大电影
漫画:程序员的工作
漫画:程序员的工作
那些争议最大的编程观点
那些争议最大的编程观点
做程序猿的老婆应该注意的一些事情
做程序猿的老婆应该注意的一些事情
为啥Android手机总会越用越慢?
为啥Android手机总会越用越慢?
2013年美国开发者薪资调查报告
2013年美国开发者薪资调查报告
我跳槽是因为他们的显示器更大
我跳槽是因为他们的显示器更大
60个开发者不容错过的免费资源库
60个开发者不容错过的免费资源库
程序员和编码员之间的区别
程序员和编码员之间的区别
老程序员的下场
老程序员的下场
10个调试和排错的小建议
10个调试和排错的小建议
如何成为一名黑客
如何成为一名黑客
我的丈夫是个程序员
我的丈夫是个程序员
程序员的鄙视链
程序员的鄙视链
程序员都该阅读的书
程序员都该阅读的书
为什么程序员都是夜猫子
为什么程序员都是夜猫子
程序员的一天:一寸光阴一寸金
程序员的一天:一寸光阴一寸金
我是如何打败拖延症的
我是如何打败拖延症的
那些性感的让人尖叫的程序员
那些性感的让人尖叫的程序员
写给自己也写给你 自己到底该何去何从
写给自己也写给你 自己到底该何去何从
什么才是优秀的用户界面设计
什么才是优秀的用户界面设计
Web开发人员为什么越来越懒了?
Web开发人员为什么越来越懒了?
Google伦敦新总部 犹如星级庄园
Google伦敦新总部 犹如星级庄园
程序员眼里IE浏览器是什么样的
程序员眼里IE浏览器是什么样的
 程序员的样子
程序员的样子
程序员周末都喜欢做什么?
程序员周末都喜欢做什么?
中美印日四国程序员比较
中美印日四国程序员比较
“懒”出效率是程序员的美德
“懒”出效率是程序员的美德
不懂技术不要对懂技术的人说这很容易实现
不懂技术不要对懂技术的人说这很容易实现
软件开发程序错误异常ExceptionCopyright © 2009-2015 MyException 版权所有