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));
}
相关推荐
----云烟----17 分钟前
QT中QString类的各种使用
开发语言·qt
lsx20240622 分钟前
SQL SELECT 语句:基础与进阶应用
开发语言
开心工作室_kaic1 小时前
ssm161基于web的资源共享平台的共享与开发+jsp(论文+源码)_kaic
java·开发语言·前端
向宇it1 小时前
【unity小技巧】unity 什么是反射?反射的作用?反射的使用场景?反射的缺点?常用的反射操作?反射常见示例
开发语言·游戏·unity·c#·游戏引擎
武子康1 小时前
Java-06 深入浅出 MyBatis - 一对一模型 SqlMapConfig 与 Mapper 详细讲解测试
java·开发语言·数据仓库·sql·mybatis·springboot·springcloud
九鼎科技-Leo1 小时前
什么是 WPF 中的依赖属性?有什么作用?
windows·c#·.net·wpf
转世成为计算机大神1 小时前
易考八股文之Java中的设计模式?
java·开发语言·设计模式
宅小海2 小时前
scala String
大数据·开发语言·scala
qq_327342732 小时前
Java实现离线身份证号码OCR识别
java·开发语言
锅包肉的九珍2 小时前
Scala的Array数组
开发语言·后端·scala