C#,《小白学程序》第二十五课:大数乘法(BigInteger Multiply)的Karatsuba算法及源代码

1 文本格式

/// <summary>

/// 《小白学程序》第二十五课:大数(BigInteger)的Karatsuba乘法

/// Multiplies two bit strings X and Y and returns result as long integer

/// </summary>

/// <param name="a"></param>

/// <param name="b"></param>

/// <returns></returns>

public static string karatsuba_multiply(string a, string b)

{

// Find the maximum of lengths of x and Y and make

// length of smaller string same as that of larger string

int n = Math.Max(a.Length, b.Length);

if (a.Length != n) a = a.PadLeft(n, '0');

if (b.Length != n) b = b.PadLeft(n, '0');

// Base cases

if (n == 0)

{

return "0";

}

else if (n == 1)

{

return ((a[0] - '0') * (b[0] - '0')).ToString();

}

else if (n <= 9)

{

// int max 21474 83647

// long max 9223 37203 68547 75807

return (ulong.Parse(a) * ulong.Parse(b)).ToString();

}

int fh = n / 2; // First half of string

int sh = n - fh; // Second half of string

// Find the first half and second half of first string.

string x1 = a.Substring(0, fh);

string x2 = a.Substring(fh);

// Find the first half and second half of second string

string y1 = b.Substring(0, fh);

string y2 = b.Substring(fh);

// Recursively calculate the three products of

// inputs of size n/2

string p1 = karatsuba_multiply(x1, y1);

string p2 = karatsuba_multiply(x2, y2);

//string P3 = Karatsuba(AddBitStrings(Xl, Xr), AddBitStrings(Yl, Yr));

string p3 = karatsuba_multiply(big_integer_plus(x1, x2), big_integer_plus(y1, y2));

// Combine the three products to get the final result.

//return P1 * (1L << (2 * sh)) + (P3 - P1 - P2) * (1L << sh) + P2;

int[] t1 = new int[sh];

string w1 = String.Join("", t1);

string v1 = p1 + w1 + w1;

string v2 = big_integer_subtract(p3, big_integer_plus(p1, p2)) + w1;

return big_integer_plus(v1, big_integer_plus(v2, p2));

}

2 代码格式

cs 复制代码
/// <summary>
/// 《小白学程序》第二十五课:大数(BigInteger)的Karatsuba乘法
/// Multiplies two bit strings X and Y and returns result as long integer
/// </summary>
/// <param name="a"></param>
/// <param name="b"></param>
/// <returns></returns>
public static string karatsuba_multiply(string a, string b)
{
    // Find the maximum of lengths of x and Y and make
    // length of smaller string same as that of larger string
    int n = Math.Max(a.Length, b.Length);
    if (a.Length != n) a = a.PadLeft(n, '0');
    if (b.Length != n) b = b.PadLeft(n, '0');

    // Base cases
    if (n == 0)
    {
        return "0";
    }
    else if (n == 1)
    {
        return ((a[0] - '0') * (b[0] - '0')).ToString();
    }
    else if (n <= 9)
    {
        // int max 21474 83647
        // long max 9223 37203 68547 75807
        return (ulong.Parse(a) * ulong.Parse(b)).ToString();
    }

    int fh = n / 2;  // First half of string
    int sh = n - fh; // Second half of string

    // Find the first half and second half of first string.
    string x1 = a.Substring(0, fh);
    string x2 = a.Substring(fh);

    // Find the first half and second half of second string
    string y1 = b.Substring(0, fh);
    string y2 = b.Substring(fh);

    // Recursively calculate the three products of
    // inputs of size n/2
    string p1 = karatsuba_multiply(x1, y1);
    string p2 = karatsuba_multiply(x2, y2);
    //string P3 = Karatsuba(AddBitStrings(Xl, Xr), AddBitStrings(Yl, Yr));
    string p3 = karatsuba_multiply(big_integer_plus(x1, x2), big_integer_plus(y1, y2));

    // Combine the three products to get the final result.
    //return P1 * (1L << (2 * sh)) + (P3 - P1 - P2) * (1L << sh) + P2;
    int[] t1 = new int[sh];
    string w1 = String.Join("", t1);
    string v1 = p1 + w1 + w1;
    string v2 = big_integer_subtract(p3, big_integer_plus(p1, p2)) + w1;
    return big_integer_plus(v1, big_integer_plus(v2, p2));
}
相关推荐
前端世界15 小时前
HarmonyOS 数据处理性能优化:算法 + 异步 + 分布式实战
算法·性能优化·harmonyos
楼田莉子15 小时前
C++算法专题学习:栈相关的算法
开发语言·c++·算法·leetcode
晨非辰16 小时前
#C语言——刷题攻略:牛客编程入门训练(九):攻克 分支控制(三)、循环控制(一),轻松拿捏!
c语言·开发语言·经验分享·学习方法·visual studio
_oP_i16 小时前
Java 服务接口中解决跨域(CORS,Cross-Origin Resource Sharing)问题
java·开发语言
陈序猿(代码自用版)16 小时前
【考研C语言编程题】数组元素批量插入实现(含图示+三部曲拆解)
c语言·开发语言·考研
唐•苏凯16 小时前
ArcGIS Pro 遇到严重的应用程序错误而无法启动
开发语言·javascript·ecmascript
kyle~16 小时前
排序---冒泡排序(Bubble Sort)
c语言·c++·算法
l1t16 小时前
我改写的二分法XML转CSV文件程序速度追上了张泽鹏先生的
xml·c语言·人工智能·算法·expat
一碗白开水一16 小时前
【论文阅读】Far3D: Expanding the Horizon for Surround-view 3D Object Detection
论文阅读·人工智能·深度学习·算法·目标检测·计算机视觉·3d
轮到我狗叫了16 小时前
力扣.1054距离相等的条形码力扣767.重构字符串力扣47.全排列II力扣980.不同路径III力扣509.斐波那契数列(记忆化搜索)
java·算法·leetcode