热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

[大整数乘法]java代码实现

本文介绍了使用java代码实现大整数乘法的过程,同时也涉及到大整数加法和大整数减法的计算方法。通过分治算法来提高计算效率,并对算法的时间复杂度进行了研究。详细代码实现请参考文章链接。

上一篇写的“[大整数乘法]分治算法的时间复杂度研究”,这一篇是基于上一篇思想的代码实现,以下是该文章的连接:

http://www.cnblogs.com/McQueen1987/p/3348426.html

代码主要实现大整数乘法,过程中也涉及到[大整数加法] 和 [大整数减法] 的计算,代码如下:

 

类1

————————————————————————————————————————————————————————————

package bigIntNum;

public class NumDividEqual { public char[] A;
public char[] B;
int n;
/**
* 将数组均分为两份,分别存入数组A和数组B中;
* @param input
*/
public NumDividEqual(char[] input){
n = input.length/2;
A = new char[n];
B = new char[n];
for(int i = 0; i A[i] = input[i];
}
for(int i = 0; i B[i] = input[i + n];
}
}

public static void main(String[] args) {
// TODO Auto-generated method stub
}
}

  

类2

————————————————————————————————————————————————————————————

package bigIntNum;

import java.util.Arrays;

public class bigIntMult {
/**
* 将字符数组倒序排列
*
@param input
*
@return
*/
public char[] reverse(char[] input) {
char[] output = new char[input.length];
for (int i = 0; i ) {
output[i] = input[input.length - 1 - i];
}
return output;
}

/**
* 将大整数平均分成两部分
*
@param input
*
@return
*/
public NumDividEqual partition(char[] input) {
return new NumDividEqual(input);
}

/**
* 求两数组中较大数组的长度,如果其长度为奇数则+1变偶
*
@param num1
*
@param num2
*
@return
*/
public int calLength(char[] num1, char[] num2) {
int len = num1.length > num2.length ? num1.length : num2.length;
if (len == 1)
return 1;
len
+= len & 1;
return len;
}

/**
* 除去数字前面多余的0
*
@param input
*
@return
*/
public static char[] trimPrefix(char[] input) {
char[] ret = null;
for (int i = 0; i ) {
if (ret == null && input[i] == '0')
continue;
else {
if (ret == null) {
ret
= new char[input.length - i];//出去数字前面多余的0
}
ret[i
- (input.length - ret.length)] = input[i];
}
}
if (ret == null)
return new char[] { '0' };
return ret;
}

/**
* 数组如果长度不足n,则在数组前面补0,使长度为n。
*
@param input 输入数组要求数字的最高位存放在数组下标最小位置
*
@param n
*
@return
*/
public static char[] format(char[] input, int n) {//
if (input.length >= n) {
return input;
}
char[] ret = new char[n];
for (int i = 0; i ) {
ret[i] = '0';
}
for (int i = 0; i ) {
ret[n - input.length + i] = input[i];
}
return ret;
}

/**
* 大整数尾部补0。相当于移位,扩大倍数
*
@param input
*
@param n
*
@return
*/
public char[] addTail(char[] input, int n) {//
char[] ret = new char[input.length + n];
for (int i = 0; i ) {
ret[i] = input[i];
}
for (int i = input.length; i ) {
ret[i] = '0';
}
return ret;
}

/**
* 大整数加法
*
@param num1
*
@param num2
*
@return
*/
public char[] add(char[] num1, char[] num2) {
int len = num2.length > num1.length ? num2.length : num1.length;
int carry = 0;//进位标识
num1 = format(num1, len);
num2
= format(num2, len);
char[] ret = new char[len + 1];

for (int i = len - 1; i >= 0; i--) {
int tmp = num1[i] + num2[i] - 96;
tmp
+= carry;
if (tmp >= 10) {
carry
= 1;
tmp
= tmp - 10;
}
else {
carry
= 0;
}
ret[len
- i - 1] = (char) (tmp + 48);
}
ret[len]
= (char) (carry + 48);//最后一次,最高位的进位
return trimPrefix(reverse(ret));
}
/**
* 大整数减法:
*
@param num1 被减数,大整数乘法中只有一个减法(A+B)(C+D)-(AC+BD)=AC+BC>0,因此參數num1>num2且都为正
*
@param num2 减数
*
@return
*/
public static char[] sub(char[] num1, char[] num2) {
int lenMax = num1.length > num2.length ? num1.length : num2.length;
char[] newNum1 = Arrays.copyOf(format(num1, lenMax), lenMax);//字符串前面补0,使两串长度相同
char[] newNum2 = Arrays.copyOf(format(num2, lenMax), lenMax);

for(int i=0;i//when num1-num2<0 return
if((newNum1[i]=='0' && newNum1[i]=='0') || newNum1[i] == newNum2[i]){//newNum1 is bigger;
continue;
}
else if(newNum1[i] //不滿足參數num1>num2;
System.out.println("The Parameter in sub(A,B).A MUST Bigger Than B!");
System.exit(
0);
}
else break;
}

for(int i=lenMax-1;i>=0;i--){
if(newNum1[i] //result <0
newNum1[i] = (char) (newNum1[i] + '0' + 10 - newNum2[i]);
newNum1[i
-1] = (char) (newNum1[i-1] - 1);
}
else{
newNum1[i]
= (char) (newNum1[i] + '0' - newNum2[i]);
}
}
return trimPrefix(newNum1);
}
/**
* 大整数乘法
*
@param num1
*
@param num2
*
@return
*/
public char[] mult(char[] num1, char[] num2) {
char[] A, B, C, D, AC, BD, AjB, CjD, ACjBD, AjBcCjD, SUM;
int N = calLength(num1, num2);//求两数组中较大数组的长度,如果长度为奇数则+1变偶,方便二分成两部分
num1 = format(num1, N);//数组高位存整数的高位数;数字前面补0,使长度为n;
num2 = format(num2, N);
if (num1.length > 1) {
NumDividEqual nu1
= partition(num1);//将大整数平均分成两部分
NumDividEqual nu2 = partition(num2);
A
= nu1.A;
B
= nu1.B;
C
= nu2.A;
D
= nu2.B;
AC
= mult(A, C);//分治求大整数乘法
BD = mult(B, D);
AjB
= add(A,B);
CjD
= add(C,D);
ACjBD
= add(AC,BD);
AjBcCjD
= mult(AjB, CjD);

char[] tmp1 = addTail(sub(AjBcCjD, ACjBD), N / 2);//尾部补0,相当于移位
char[] tmp2 = add(addTail(AC, N), BD);
SUM
= add(tmp1, tmp2);
char[] test = trimPrefix(SUM);//除去结果前面多余的0
return test;
}
else {
Integer ret
= (num1[0] - 48) * (num2[0] - 48);
return ret.toString().toCharArray();
}
}

public static void main(String[] args) {
String st1
= "168746315641347979798";
String st2
= "164681654767446887797451316158";
char[] a = st1.toCharArray();
char[] b = st2.toCharArray();
bigIntMult bg
= new bigIntMult();
char[] ret = bg.mult(a, b);
System.out.println(ret);
}
}

 

 

代码优化:

1.可以写个hash表,存储一些常见的乘法,从而避免每次都重复计算。比如9999x9999999,里面有重复的9x9计算,可以考虑hash表存储这些计算的结果,用到时直接查询结果,从而避免重复计算。

 

声明:代码是本人脑力劳动结果,请尊重他人劳动。转载注明出处!

 

 


推荐阅读
  • 本文介绍了南邮ctf-web的writeup,包括签到题和md5 collision。在CTF比赛和渗透测试中,可以通过查看源代码、代码注释、页面隐藏元素、超链接和HTTP响应头部来寻找flag或提示信息。利用PHP弱类型,可以发现md5('QNKCDZO')='0e830400451993494058024219903391'和md5('240610708')='0e462097431906509019562988736854'。 ... [详细]
  • 本文介绍了如何在给定的有序字符序列中插入新字符,并保持序列的有序性。通过示例代码演示了插入过程,以及插入后的字符序列。 ... [详细]
  • 《数据结构》学习笔记3——串匹配算法性能评估
    本文主要讨论串匹配算法的性能评估,包括模式匹配、字符种类数量、算法复杂度等内容。通过借助C++中的头文件和库,可以实现对串的匹配操作。其中蛮力算法的复杂度为O(m*n),通过随机取出长度为m的子串作为模式P,在文本T中进行匹配,统计平均复杂度。对于成功和失败的匹配分别进行测试,分析其平均复杂度。详情请参考相关学习资源。 ... [详细]
  • 关键词:Golang, Cookie, 跟踪位置, net/http/cookiejar, package main, golang.org/x/net/publicsuffix, io/ioutil, log, net/http, net/http/cookiejar ... [详细]
  • 动态规划算法的基本步骤及最长递增子序列问题详解
    本文详细介绍了动态规划算法的基本步骤,包括划分阶段、选择状态、决策和状态转移方程,并以最长递增子序列问题为例进行了详细解析。动态规划算法的有效性依赖于问题本身所具有的最优子结构性质和子问题重叠性质。通过将子问题的解保存在一个表中,在以后尽可能多地利用这些子问题的解,从而提高算法的效率。 ... [详细]
  • 个人学习使用:谨慎参考1Client类importcom.thoughtworks.gauge.Step;importcom.thoughtworks.gauge.T ... [详细]
  • CF:3D City Model(小思维)问题解析和代码实现
    本文通过解析CF:3D City Model问题,介绍了问题的背景和要求,并给出了相应的代码实现。该问题涉及到在一个矩形的网格上建造城市的情景,每个网格单元可以作为建筑的基础,建筑由多个立方体叠加而成。文章详细讲解了问题的解决思路,并给出了相应的代码实现供读者参考。 ... [详细]
  • 猜字母游戏
    猜字母游戏猜字母游戏——设计数据结构猜字母游戏——设计程序结构猜字母游戏——实现字母生成方法猜字母游戏——实现字母检测方法猜字母游戏——实现主方法1猜字母游戏——设计数据结构1.1 ... [详细]
  • Java太阳系小游戏分析和源码详解
    本文介绍了一个基于Java的太阳系小游戏的分析和源码详解。通过对面向对象的知识的学习和实践,作者实现了太阳系各行星绕太阳转的效果。文章详细介绍了游戏的设计思路和源码结构,包括工具类、常量、图片加载、面板等。通过这个小游戏的制作,读者可以巩固和应用所学的知识,如类的继承、方法的重载与重写、多态和封装等。 ... [详细]
  • 本文详细介绍了GetModuleFileName函数的用法,该函数可以用于获取当前模块所在的路径,方便进行文件操作和读取配置信息。文章通过示例代码和详细的解释,帮助读者理解和使用该函数。同时,还提供了相关的API函数声明和说明。 ... [详细]
  • Iamtryingtomakeaclassthatwillreadatextfileofnamesintoanarray,thenreturnthatarra ... [详细]
  • 本文介绍了使用Java实现大数乘法的分治算法,包括输入数据的处理、普通大数乘法的结果和Karatsuba大数乘法的结果。通过改变long类型可以适应不同范围的大数乘法计算。 ... [详细]
  • Mac OS 升级到11.2.2 Eclipse打不开了,报错Failed to create the Java Virtual Machine
    本文介绍了在Mac OS升级到11.2.2版本后,使用Eclipse打开时出现报错Failed to create the Java Virtual Machine的问题,并提供了解决方法。 ... [详细]
  • 本文介绍了一种划分和计数油田地块的方法。根据给定的条件,通过遍历和DFS算法,将符合条件的地块标记为不符合条件的地块,并进行计数。同时,还介绍了如何判断点是否在给定范围内的方法。 ... [详细]
  • 本文介绍了Perl的测试框架Test::Base,它是一个数据驱动的测试框架,可以自动进行单元测试,省去手工编写测试程序的麻烦。与Test::More完全兼容,使用方法简单。以plural函数为例,展示了Test::Base的使用方法。 ... [详细]
author-avatar
Larry_He
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有