【(阶乘的质因数分解)算组合数】【TOJ4111】【Binomial efficient】
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);
}
}
利用阶乘的质因数分解!
比如250!
1*2*3*4*5*6*7*8*9*10*11*12*13*14.....250
中3的质因子个数 除了3后变成(不是倍数的不管)计算3^1次方的为250/3个
又变成 1 2 3 .....250除3
重复上面 知道 3^2 为250/(3^2)
所以阶乘的质因数分解是另外的简单算法
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】的更多相关文章
- P2043 质因子分解(阶乘的质因数分解)
P2043 质因子分解 对$n!$进行质因数分解的一种高效算法 首先,筛出$<=n$的素数 蓝后,对$n$反复除以$prime$,同时$cnt+=n/prime$ $n!$中含有该$prime$ ...
- 济南学习D3T1__线性筛和阶乘质因数分解
[问题描述] 从1− N中找一些数乘起来使得答案是一个完全平方数,求这个完全平方数最大可能是多少. [输入格式] 第一行一个数字N. [输出格式] 一行,一个整数代表答案对100000007取模之后的 ...
- 对于n!的快速质因数分解
N!的阶乘的质因数分解 对于N的阶乘 比如8! 我们要算其中一个质因数出现次数 我们注意到 8!=1 2 3 4 5 6 7 8 1 1 1 1 2的倍数出现的次数8/2=4 1 1 4的倍数出现的次 ...
- 【BZOJ2227】【ZJOI2011】看电影 [组合数][质因数分解]
看电影 Time Limit: 10 Sec Memory Limit: 259 MB[Submit][Status][Discuss] Description 到了难得的假期,小白班上组织大家去看 ...
- poj 3421 X-factor Chains——质因数分解
题目:http://poj.org/problem?id=3421 记忆化搜索竟然水过去了.仔细一想时间可能有点不对,但还是水过去了. #include<iostream> #includ ...
- Poj 1401 Factorial(计算N!尾数0的个数——质因数分解)
一.Description The most important part of a GSM network is so called Base Transceiver Station (BTS). ...
- POj3421 X-factor Chains(质因数分解+排列组合)
POj3421X-factor Chains 一开始没读懂题意,不太明白 Xi | Xi+1 where a | b means a perfectly divides into b的意思,后来才发现 ...
- HDU3988-Harry Potter and the Hide Story(数论-质因数分解)
Harry Potter and the Hide Story Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/65536 ...
- UVA10375 Choose and divide 质因数分解
质因数分解: Choose and divide Time Limit: 3000MS Memory Limit: Unknown 64bit IO Format: %lld & %l ...
随机推荐
- 谋哥:研究App排行榜浮出的神器
昨天发的<App排行榜的秘密>到头条网,阅读量到2万,踩的比顶的多几倍.原因是由于我使用360市场的数据来分析,而且这帮喷子根本不看你分析数据背后的意义,反正看到自己不喜欢的比方" ...
- Hadoop中SequenceFile的使用
1.对于某些应用而言,须要特殊的数据结构来存储自己的数据. 对于基于MapReduce的数据处理.将每一个二进制数据的大对象融入自己的文件里并不能实现非常高的可扩展性,针对上述情况,Hadoop开发了 ...
- MFC多线程内存泄漏问题&解决方法
在用visual studio进行界面编程时(如MFC),前台UI我们能够通过MFC的消息循环机制实现.而对于后台的数据处理.我们可能会用到多线程来处理. 那么对于大多数人(尤其是我这样的菜鸟),一个 ...
- crm使用soap更改下拉框的文本值
//C#代码 //UpdateStateValueRequest updateStateValue = new UpdateStateValueRequest //{ // AttributeL ...
- Node.js实践第一天
实践案例一 表单提交 demo.html <!doctype html> <html> <head> <title>表单提交</title> ...
- DWZ在APS.NET WebForm中的使用(二)
任何框架由于个人理解不到位或者框架自身的局限性,在项目实施中,大家或多或少都会遇到一些问题,下面我就讲述下我在使用DWZ开发过程中的几个问题.如有一点能帮助到你,这篇文章也算有存在的意义了. 1.树菜 ...
- oracle 数据库学习
1.更改数据库用户名密码: alter user '用户名' identified by ’password'. 2.查看数据库有哪些用户:sqlplus system/password ;selec ...
- iOS 事件处理机制与图像渲染过程
Peter在开发公众号功能时触发了一个bug,导致群发错误.对此我们深表歉意,并果断开除了Peter.以下交回给正文时间: iOS 事件处理机制与图像渲染过程 iOS RunLoop都干了什么 iOS ...
- 基本的 html 代码
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...
- 时间的获取和转换,C#和Sql
我们可以通过使用DataTime这个类来获取当前的时间.通过调用类中的各种方法我们可以获取不同的时间:如:日期(2008-09-04).时间(12:12:12).日期+时间(2008-09-04 12 ...