日期:2014-05-20 浏览次数:21299 次
public static uint[] CalculateLargeNumber(int n)
{
if (n < 0) { throw new ArgumentOutOfRangeException("n必须为非负数。"); }
if (n == 0 || n == 1) { return new uint[] { 1 }; }
// 数组的最大长度
const int MaxLength = 100000;
uint[] array = new uint[MaxLength];
// 1! = 1
array[0] = 1;
int i = 0;
int j = 0;
// valid为当前阶乘值的位数(如5! = 120,此时valid = 3)
int valid = 1;
for (i = 2; i <= n; i++)
{
long carry = 0;
for (j = 0; j < valid; j++)
{
long multipleResult = array[j] * i + carry;
// 计算当前位的数值
array[j] = (uint)(multipleResult % 10);
// 计算进到高位的数值
carry = multipleResult / 10;
}
// 为更高位赋值
while (carry != 0)
{
array[valid++] = (uint)(carry % 10);
carry /= 10;
}
}
// 截取有效元素
uint[] result = new uint[valid];
Array.Copy(array, result, valid);
return result;
}
static void Main(string[] args)
{
uint[] arr = CalculateLargeNumber(100);
string str = "";
for (int i = arr.Length-1; i >=0 ; i--)
{
str += arr[i].ToString();
}
Console.WriteLine(str);
Console.ReadKey();
}
------解决方案--------------------
using System;
using System.Text;
using System.Collections.Generic;
namespace Skyiv
{
class Factorial
{
const int n = 100;
static void Main()
{
BigInteger f = 1;
for (int i = 2; i <= n; i++)
{
f *= i;
}
Console.WriteLine("{0}! = {1}", n, f);
}
}
static class Utility
{
public static void Swap<T>(ref T x, ref T y)
{
T z = x;
x = y;
y = z;
}
}
sealed class BigInteger : IEquatable<BigInteger>, IComparable<BigInteger>
{
static readonly int Len = 9; // int.MaxValue = 2,147,483,647
static readonly int Base = (int)Math.Pow(10, Len);
int sign;
List<int> data;
BigInteger(long x)
{
sign = (x == 0) ? 0 : ((x > 0) ? 1 : -1);
data = new List<int>();
for (ulong z = (x < 0) ? (ulong)-x : (ulong)x; z != 0; z /= (ulong)Base) data.Add((int)(z % (ulong)Base));
}
BigInteger(BigInteger x)
{
sign = x.sign; // x != null
data = new List<int>(x.data);
}
public static implicit operator BigInteger(long x)
{
return new BigInteger(x);
}
public static BigInteger Parse(string s)
{
if (s == null) return null;
s = s.Trim().Replace(",", null);
if (s.Length == 0) return 0;
BigInteger z = 0;
z.sign = (s[0] == '-') ? -1 : 1;
if (s[0] == '-' || s[0] == '+') s = s.Substring(1);
int r = s.Length % Len;
z.data = new List<int>(new int[s.Length / Len + ((r != 0) ? 1 : 0)]);
int i = z.data.Count - 1;
if (r != 0) z.data[i--] = int.Parse(s.Substring(0, r));
for (; i >= 0; i--, r += Len) z.data[i] = int.Parse(s.Substring(r, Len));
z.Shrink();
return z;
}
public static BigInteger Abs(BigInteger x)
{
if (x == null) return null;
BigInteger z = new BigInteger(x);
z.sign = Math.Abs(x.sign);
return z;
}
public static BigInteger Pow(BigInteger x, long y)
{
if (x == null) return null;
BigInteger z = 1, n = x;
for (; y > 0; y >>= 1, n *= n) if ((y & 1) != 0) z *= n;
return z;
}
public static BigInteger operator +(BigInteger x)
{
if (x == null) return null;
return new BigInteger(x);
}
public static BigInteger operator -(BigInteger x)
{
if (x == null) return null;
BigInteger z = new BigInteger(x);
z.sign = -x.sign;
return z;
}
public static BigInteger operator ++(BigInteger x)
{
return x + 1;
}
public static BigInteger operator --(BigInteger x)
{
return x - 1;
}
public static BigInteger operator +(BigInteger x, BigInteger y)
{
if (x == null || y == null) return null;
if (x.AbsCompareTo(y) < 0) Utility.Swap(ref x, ref y);
BigInteger z = new BigInteger(x);
if (x.sign * y.sign == -1) z.AbsSubtract(y);
else z.AbsAdd(y);
return z;
}
public static BigInteger operator -(BigInteger x, BigInteger y)
{
if (x == null || y == null) return null;
return x + (-y);
}
public static BigInteger operator *(BigInteger x, BigInteger y)
{
if (x == null || y == null) return null;
BigInteger z = 0;
z.sign = x.sign * y.sign;
z.data = new List<int>(new int[x.data.Count + y.data.Count]);
for (int i = x.data.Count - 1; i >= 0; i--)
for (int j = y.data.Count - 1; j >= 0; j--)
{
long n = Math.BigMul(x.data[i], y.data[j]);
z.data[i + j] += (int)(n % Base);
z.CarryUp(i + j);
z.data[i + j + 1] += (int)(n / Base);
z.CarryUp(i + j + 1);
}
z.Shrink();
return z;
}
public static BigInteger operator /(BigInteger dividend, BigInteger divisor)
{
BigInteger remainder;
return DivRem(dividend, divisor, out remainder);
}
public static BigInteger operator %(BigInteger dividend, BigInteger divisor)
{
BigInteger remainder;
DivRem(dividend, divisor, out remainder);
return remainder;
}
public static BigInteger DivRem(BigInteger dividend, BigInteger divisor, out BigInteger remainder)
{
remainder = null;
if (dividend == null || divisor == null) return null;
if (divisor.sign == 0) throw new DivideByZeroException();
if (dividend.AbsCompareTo(divisor) < 0)
{
remainder = new BigInteger(dividend);
return 0;
}
remainder = 0;
BigInteger quotient = 0;
quotient.sign = dividend.sign * divisor.sign;
quotient.data = new List<int>(new int[dividend.data.Count]);
divisor = Abs(divisor); // NOT: divisor.sign = Math.Abs(divisor.sign);
for (int i = dividend.data.Count - 1; i >= 0; i--)
{
remainder = remainder * Base + dividend.data[i];
int iQuotient = remainder.BinarySearch(divisor, -1);
quotient.data[i] = iQuotient;
remainder -= divisor * iQuotient;
}
quotient.Shrink();
if (remainder.sign != 0) remainder.sign = dividend.sign;
return quotient;
}
public static BigInteger Sqrt(BigInteger x)
{
if (x == null || x.sign < 0) return null;
BigInteger root = 0;
root.sign = 1;
root.data = new List<int>(new int[x.data.Count / 2 + 1]);
for (int i = root.data.Count - 1; i >= 0; i--) root.data[i] = x.BinarySearch(root, i);
root.Shrink();
return root;
}
int BinarySearch(BigInteger divisor, int i)
{
int low = 0, high = Base - 1, mid = 0, cmp = 0;
while (low <= high)
{
mid = (low + high) / 2;
cmp = CompareTo(divisor, mid, i);
if (cmp > 0) low = mid + 1;
else if (cmp < 0) high = mid - 1;
else return mid;
}
return (cmp < 0) ? (mid - 1) : mid;
}
int CompareTo(BigInteger divisor, int mid, int i)
{
if (i >= 0) divisor.data[i] = mid;
return AbsCompareTo(divisor * ((i >= 0) ? divisor : mid));
}
void AbsAdd(BigInteger x)
{
for (int i = 0; i < x.data.Count; i++)
{
data[i] += x.data[i];
CarryUp(i);
}
}
void AbsSubtract(BigInteger x)
{
for (int i = 0; i < x.data.Count; i++)
{
data[i] -= x.data[i];
CarryDown(i);
}
Shrink();
}
void CarryUp(int n)
{
for (; data[n] >= Base; n++)
{
if (n == data.Count - 1) data.Add(data[n] / Base);
else data[n + 1] += data[n] / Base;
data[n] %= Base;
}
}
void CarryDown(int n)
{
for (; data[n] < 0; n++)
{
data[n + 1]--;
data[n] += Base;
}
Shrink();
}
void Shrink()
{
for (int i = data.Count - 1; i >= 0 && data[i] == 0; i--) data.RemoveAt(i);
if (data.Count == 0) sign = 0;
}
public static bool operator ==(BigInteger x, BigInteger y)
{
if (object.ReferenceEquals(x, null)) return object.ReferenceEquals(y, null);
return x.Equals(y);
}
public static bool operator !=(BigInteger x, BigInteger y)
{
if (object.ReferenceEquals(x, null)) return !object.ReferenceEquals(y, null);
return !x.Equals(y);
}
public static bool operator <(BigInteger x, BigInteger y)
{
if (object.ReferenceEquals(x, null)) return !object.ReferenceEquals(y, null);
return x.CompareTo(y) < 0;
}
public static bool operator >(BigInteger x, BigInteger y)
{
if (object.ReferenceEquals(x, null)) return false;
return x.CompareTo(y) > 0;
}
public static bool operator <=(BigInteger x, BigInteger y)
{
if (object.ReferenceEquals(x, null)) return true;
return x.CompareTo(y) <= 0;
}
public static bool operator >=(BigInteger x, BigInteger y)
{
if (object.ReferenceEquals(x, null)) return object.ReferenceEquals(y, null);
return x.CompareTo(y) >= 0;
}
public override string ToString()
{
StringBuilder sb = new StringBuilder();
if (sign < 0) sb.Append('-');
sb.Append((data.Count == 0) ? 0 : data[data.Count - 1]);
for (int i = data.Count - 2; i >= 0; i--) sb.Append(data[i].ToString("D" + Len));
return sb.ToString();
}
public override int GetHashCode()
{
int hash = sign;
foreach (int n in data) hash ^= n;
return hash;
}
public override bool Equals(object other)
{
if (other == null || GetType() != other.GetType()) return false;
return Equals(other as BigInteger);
}
public bool Equals(BigInteger other)
{
return CompareTo(other) == 0;
}
public int CompareTo(BigInteger other)
{
if (object.ReferenceEquals(other, null)) return 1;
if (sign < other.sign) return -1;
if (sign > other.sign) return 1;
if (sign == 0) return 0;
return sign * AbsCompareTo(other);
}
int AbsCompareTo(BigInteger other)
{
if (data.Count < other.data.Count) return -1;
if (data.Count > other.data.Count) return 1;
for (int i = data.Count - 1; i >= 0; i--)
if (data[i] != other.data[i])
return (data[i] < other.data[i]) ? -1 : 1;
return 0;
}
}
}