n<=10^6

m<=10^6

p=2^32

用unsigned int 可以避免取模

我写的SB超时 阶乘分解代码

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <cstring>

#include <ctime>

#include <algorithm>

#include <iostream>

#include <sstream>

#include <string>

#define oo 0x13131313

using namespace std;

const unsigned int N=1000000+5;

unsigned int tag[N],p[N],z[N],mm[N];

unsigned int cnt = 0;

unsigned int n,m;

unsigned int quickpow(unsigned int m,unsigned int n)

{

    unsigned int b = 1;

    while (n > 0)

    {

          if (n & 1)

             b = (b*m);

          n = n >> 1 ;

          m = (m*m);

    }

    return b;

}

void get_prime()

{

    tag[1]=1;

    tag[0]=1;

    for (unsigned int i = 2; i < N; i++)

    {

        if (!tag[i])    p[cnt++] = i;

        for (unsigned int j = 0; j < cnt && p[j] * i < N; j++)

        {

            tag[i*p[j]] = 1;

            if (i % p[j] == 0)

                break;

        }

    }

}

unsigned int FIND(unsigned int x)

{

    unsigned int l=0,r=cnt-1;

    while(l<=r)

    {

        unsigned int m=(l+r)/2;

        if(p[m]==x) return m;

        else if(p[m]<x) l=m+1;

        else r=m-1;

    }

}

void fenjie(unsigned int *a,unsigned int d)

{

    for(unsigned int i=0;i<cnt&&p[i]<=sqrt(d);i++)

    {

        while(d%p[i]==0)

        {

            d=d/p[i];

            a[i]++;

        }

    }

    if(tag[d]==0)

    {

        unsigned int t=FIND(d);

        a[t]++;

    }

}

int main()

{

   // freopen("a.in","r",stdin);

    get_prime();

    int T;

    cin>>T;

    while(T--)

    {

        cin>>n>>m;

        memset(z,0,sizeof(z));

        memset(mm,0,sizeof(mm));

        for(unsigned int i=n;i>=n-m+1;i--)

            fenjie(z,i);

        for(unsigned int i=1;i<=m;i++)

            fenjie(mm,i);

        unsigned int ans=1;

        for(unsigned int i=0;i<cnt;i++)

            {

                z[i]=z[i]-mm[i];

                ans=ans*quickpow(p[i],z[i]);

            }

        printf("%u\n",ans);

    }

}

利用阶乘的质因数分解！

比如２５０！

１＊２＊３＊４＊５＊６＊７＊８＊９＊１０＊１１＊１２＊１３＊１４．．．．．２５０

中３的质因子个数　除了３后变成（不是倍数的不管）计算３＾１次方的为２５０／３个

又变成　１　２　３　．．．．．２５０除３　　

重复上面　知道　３＾２　为２５０／（３＾２）　

所以阶乘的质因数分解是另外的简单算法

void getcn(int n)

{

	int ans = 0;

	int i;

	for (i = 1; i <= prime[0] && prime[i] <= n; i++)

	{

		int tmp = n;

		while (tmp)

		{

			num[i] += tmp / prime[i];

			tmp /= prime[i];

		}

	}

	num[0] = i;

}

所以最后的代码是

#include <string.h>

#include <iostream>

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

#include <string>

#include <algorithm>

#include <vector>

#include <string.h>

#include <time.h>

#include <queue>

#include <stack>

#include <map>

#include <set>

#include <sstream>

#define INF 0x3f3f3f3f

#define MAXN 1000005

#define Precision 100005

#define MAX_INT 2147483647

#define Pi acos(-1.0)

#define lowbit(x) ((x)&(-x))

#define Lson root<<1,left,mid

#define Rson root<<1|1,mid+1,right

#define LL long long

#define ULL unsigned long long

#define fresh(x) memset(x,0,sizeof(x))

using namespace std;

int prime[MAXN];

int num[MAXN];

void print()

{

<span style="white-space:pre">	</span>memset(prime, 0, sizeof(prime));

<span style="white-space:pre">	</span>for (int i = 2; i <= 1000000; i++)

<span style="white-space:pre">	</span>{

<span style="white-space:pre">		</span>if (!prime[i]) prime[++prime[0]] = i;

<span style="white-space:pre">		</span>for (int j = 1; j <= prime[0] && prime[j] <= 1000000 / i; j++)

<span style="white-space:pre">		</span>{

<span style="white-space:pre">			</span>prime[prime[j]*i] = 1;

<span style="white-space:pre">			</span>if (i % prime[j] == 0) break;

<span style="white-space:pre">		</span>}

<span style="white-space:pre">	</span>}

}

unsigned qpow(unsigned a, unsigned b)

{

<span style="white-space:pre">	</span>unsigned ans = 1;

<span style="white-space:pre">	</span>while (b)

<span style="white-space:pre">	</span>{

<span style="white-space:pre">		</span>if (b & 1)

<span style="white-space:pre">			</span>ans *= a;

<span style="white-space:pre">		</span>b >>= 1;

<span style="white-space:pre">		</span>a *= a;

<span style="white-space:pre">	</span>}

<span style="white-space:pre">	</span>return ans;

}

void getcn(int n)

{

<span style="white-space:pre">	</span>int ans = 0;

<span style="white-space:pre">	</span>int i;

<span style="white-space:pre">	</span>for (i = 1; i <= prime[0] && prime[i] <= n; i++)

<span style="white-space:pre">	</span>{

<span style="white-space:pre">		</span>int tmp = n;

<span style="white-space:pre">		</span>while (tmp)

<span style="white-space:pre">		</span>{

<span style="white-space:pre">			</span>num[i] += tmp / prime[i];

<span style="white-space:pre">			</span>tmp /= prime[i];

<span style="white-space:pre">		</span>}

<span style="white-space:pre">	</span>}

<span style="white-space:pre">	</span>num[0] = i;

}

void getcm(int m)

{

<span style="white-space:pre">	</span>int ans = 0;

<span style="white-space:pre">	</span>int i;

<span style="white-space:pre">	</span>for (i = 1; i <= prime[0] && prime[i] <= m; i++)

<span style="white-space:pre">	</span>{

<span style="white-space:pre">		</span>int tmp = m;

<span style="white-space:pre">		</span>while (tmp)

<span style="white-space:pre">		</span>{

<span style="white-space:pre">			</span>num[i] -= tmp / prime[i];

<span style="white-space:pre">			</span>tmp /= prime[i];

<span style="white-space:pre">		</span>}

<span style="white-space:pre">	</span>}

}

int main()

{

<span style="white-space:pre">	</span>int n, m, T;

<span style="white-space:pre">	</span>print();

<span style="white-space:pre">	</span>//printf("%d\n", prime[prime[0]]);

<span style="white-space:pre">	</span>//cout << prime[0];

<span style="white-space:pre">	</span>scanf("%d", &T);

<span style="white-space:pre">	</span>while (T--)

<span style="white-space:pre">	</span>{

<span style="white-space:pre">		</span>memset(num, 0, sizeof(num));

<span style="white-space:pre">		</span>scanf("%d%d", &n, &m);

<span style="white-space:pre">		</span>getcn(n);

<span style="white-space:pre">		</span>getcm(m);

<span style="white-space:pre">		</span>getcm(n - m);

<span style="white-space:pre">		</span>unsigned ans = 1;

<span style="white-space:pre">		</span>for (int i = 1; i <= num[0]; i++)

<span style="white-space:pre">		</span>{

<span style="white-space:pre">			</span>ans *= qpow(prime[i], num[i]);

<span style="white-space:pre">		</span>}

<span style="white-space:pre">		</span>printf("%u\n", ans);

<span style="white-space:pre">	</span>}

<span style="white-space:pre">	</span>return 0;

}

【（阶乘的质因数分解）算组合数】【TOJ4111】【Binomial efficient】的更多相关文章

  1. P2043 质因子分解（阶乘的质因数分解）

    P2043 质因子分解 对$n!$进行质因数分解的一种高效算法 首先,筛出$<=n$的素数 蓝后,对$n$反复除以$prime$,同时$cnt+=n/prime$ $n!$中含有该$prime$ ...

  2. 济南学习D3T1__线性筛和阶乘质因数分解

    [问题描述] 从1− N中找一些数乘起来使得答案是一个完全平方数,求这个完全平方数最大可能是多少. [输入格式] 第一行一个数字N. [输出格式] 一行,一个整数代表答案对100000007取模之后的 ...

  3. 对于n!的快速质因数分解

    N!的阶乘的质因数分解 对于N的阶乘 比如8! 我们要算其中一个质因数出现次数 我们注意到 8!=1 2 3 4 5 6 7 8 1 1 1 1 2的倍数出现的次数8/2=4 1 1 4的倍数出现的次 ...

  4. 【BZOJ2227】【ZJOI2011】看电影 [组合数][质因数分解]

    看电影 Time Limit: 10 Sec  Memory Limit: 259 MB[Submit][Status][Discuss] Description 到了难得的假期,小白班上组织大家去看 ...

  5. poj 3421 X-factor Chains——质因数分解

    题目:http://poj.org/problem?id=3421 记忆化搜索竟然水过去了.仔细一想时间可能有点不对,但还是水过去了. #include<iostream> #includ ...

  6. Poj 1401 Factorial(计算N!尾数0的个数——质因数分解）

    一.Description The most important part of a GSM network is so called Base Transceiver Station (BTS). ...

  7. POj3421 X-factor Chains(质因数分解+排列组合)

    POj3421X-factor Chains 一开始没读懂题意,不太明白 Xi | Xi+1 where a | b means a perfectly divides into b的意思,后来才发现 ...

  8. HDU3988-Harry Potter and the Hide Story（数论-质因数分解）

    Harry Potter and the Hide Story Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 ...

  9. UVA10375 Choose and divide 质因数分解

    质因数分解: Choose and divide Time Limit: 3000MS   Memory Limit: Unknown   64bit IO Format: %lld & %l ...

随机推荐

  1. 【MSP是什么】MSP认证之项目管理与项目群管理的区别

    通常所说的项目管理是指运用各种相关知识.技能.方法与工具,为满足或超越项目有关各方对项目的要求与期望,所开展的各种计划.组织.领导.控制等方面的活动.具体包括项目范围管理.项目时间管理.项目成本管理. ...

  2. IOS 掉用系统发短信

    #import <MessageUI/MessageUI.h> MFMessageComposeViewControllerDelegate #pragma mark - 调用ios系统短 ...

  3. DataReader、Table、DataSet和Entity相互转化

    public class CommonService { #region DataReader转化 /// <summary> /// 将DataReader转化为Table /// &l ...

  4. Html、CSS、JavaScript 实时效果在线编辑器 - 学习的好工具，算不算？！

    关于 二维码 与 NFC 之间的出身贫贱说 太阳火神的漂亮人生 (http://blog.csdn.net/opengl_es) 本文遵循"署名-非商业用途-保持一致"创作公用协议 ...

  5. 安装mysql时出错。无法初始化。 libstdc++.so.5

    ./bin/mysqld: error while loading shared libraries: libstdc++.so.5: cannot open shared object file: ...

  6. C#获取本机IP方法,获取本机局域网IP地址方法

    1. private void GetIP() { string hostName = Dns.GetHostName();//本机名 //System.Net.IPAddress[] address ...

  7. Foundation--NSString+NSMutableString

    NSString 字符串创建: 1.NSString *strr = @"0123456789"; 2.NSString *str = [NSString stringWithSt ...

  8. 使用FMDB，libqrencode实现二维码的生成并且保存到数据库

    /**  * 1.首先导入第三方库FMDB和libqrencode,添加库libsqlite3.tbd     

  9. 池化 - Apache Commons Pool

    对于那些创建耗时较长,或者资源占用较多的对象,比如网络连接,线程之类的资源,通常使用池化来管理这些对象,从而达到提高性能的目的.比如数据库连接池(c3p0, dbcp), java的线程池 Execu ...

  10. MySQL server has gone away报错

    1.最近做插入数据库,然后一直报一个错.mysql server has gone away.(如下图) 查了好多资料,终于解决了.. 1.可能是连接超时..进入php.ini,修改wait_time ...