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. ACdream OJ 1153 (k-GCD)

    题目链接: http://115.28.76.232/problem?pid=1153 题意: 从给定的n个数中取出k个数,使得他们的最大公约数最大,求这个最大的公约数 分析: 暴力分解不可取,我们能 ...

  2. [Android4.4.3] Nubia Z5S Mokee4.4.3 RC2.0 by syhost

    这个ROM先前在Mokee官网公布过,但一些人測试bug不少,因此已经撤下, 但又有人反馈跟之前RC1.0版的bug差点儿相同, 所以再次在网盘单独公布, 截图以及注意事项见之前的RC1.0的帖子, ...

  3. 力挺8天入门wpf【转载】

    8天入门wpf—— 第八天 最后的补充 摘要: 从这一篇往前看,其实wpf中还有很多东西没有讲到,不过我的原则还是将比较常用的知识点过一遍,如果大家熟悉了这些知识,基本功也就打的差不多了,后续可以等待 ...

  4. Android应用程序窗口（Activity）实现框架简要介绍和学习计划

    文章转载至CSDN社区罗升阳的安卓之旅,原文地址:http://blog.csdn.net/luoshengyang/article/details/8170307 前面我们学习了SurfaceFli ...

  5. webapp 性能优化

    webapp 不像传统页面,它生命周期更长,在手机端上,硬件环境并没有pc上那么好.所以性能的优化尤为重要. webapp的性能优化主要分为两个方面 网络请求优化 和 页面渲染优化 , 我们对于性能优 ...

  6. linux/module.h: No such file or directory 内核模块编译过程

    1.缺少Linux kernel头文件 To install just the headers in Ubuntu: sudo apt-get install linux-headers-$(unam ...

  7. 【模拟】【HDU1443】 Joseph

    Joseph Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Sub ...

  8. NET基础课--代码安全

    1.sqlinjection

  9. jquery中this与$this的区别

    来源:http://www.jb51.net/article/19738.htm jQuery中this与$(this)的区别 $("#textbox").hover( funct ...

  10. canvas、image src、data url、blob file conversion

    //canvas.toDataURL('image/jpeg'), and convert to blob,blob is a File Object. but UC don't support fu ...