- 1、本文档共4页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
9999阶乘
using System;
using System.Collections.Generic;
using System.Linq;
using Oyster.Math;
namespace csdnTest
{
class Program
{
static void Main(string[] args)
{
DateTime beginTime = DateTime.Now;
IntX result = Factorial(1000000);
DateTime endTime = DateTime.Now;
Console.WriteLine(endTime - beginTime);
Console.WriteLine(result);
Console.ReadKey();
}
static IntX Factorial(int n)
{
int[] counter = GetPowCounter(n);
SortedDictionaryIntX, bool sDict = new SortedDictionaryIntX, bool();
//计算幂乘并将结果压入优先队列(使用优化过的大数乘法,在计算相等规模的大数乘法时,效率最高)
for (int i = 2; i = n; i++)
{
if (counter[i] 0)
sDict.Add(IntX.Pow(i, (uint)counter[i]), false);
}
IntX valueA = 1, valueB;
//用SortedDictionary模拟优先队列进行最后的计算
while (sDict.Count 1)
{
valueA = sDict.ElementAt(0).Key;
valueB = sDict.ElementAt(1).Key;
sDict.Remove(valueA);
sDict.Remove(valueB);
sDict.Add(valueA * valueB, false);
}
return sDict.ElementAt(0).Key;
}
//做质因数分解,以便使用幂乘进行计算
static int[] GetPowCounter(int n)
{
int[] pList = GetPrime(n);
int[] pCounter = new int[n + 1];
for (int i = 0; i pList.Length; i++)
{
int k = n;
while ((k /= pList[i]) 0)
pCounter[pList[i]] += k;
}
return pCounter;
}
//生成质数列表
static int[] GetPrime(int n)
{
Listint prime = new Listint();
bool[] flags = new bool[n + 1];
for (int i = 2; i = n; i++)
{
if (!flags[i])
prime.Add(i);
for (int j = 0; j prime.
文档评论(0)