【(阶乘的质因数分解)算组合数】【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 ...
随机推荐
- hdu 5625 Clarke and chemistry
Problem Description Clarke is a patient with multiple personality disorder. One day, Clarke turned i ...
- Python 面向对象基础
面向对象编程——Object Oriented Programming,简称OOP,是一种程序设计思想.OOP把对象作为程序的基本单元,一个对象包含了数据和操作数据的函数. 面向过程的程序设计把计算机 ...
- linux内核中驱动开发常见的相似多态
#include<stdio.h> #include<stdlib.h> struct test { char name[20]; void (*func)(char *); ...
- Flex布局实践
介绍常见布局的Flex写法. 你会看到,不管是什么布局,Flex往往都可以几行命令搞定. 我只列出代码,详细的语法解释请查阅<Flex布局教程:语法篇>.我的主要参考资料是Landon S ...
- H5单页面架构:backbone + requirejs + zepto + underscore
首先,来看看整个项目结构. 跟上一篇angular类似,libs里多了underscore和zepto.三个根目录文件: index.html:唯一的html main.js:requirejs的配置 ...
- Maven 工程下 Spring MVC 站点配置 (二) Mybatis数据操作
详细的Spring MVC框架搭配在这个连接中: Maven 工程下 Spring MVC 站点配置 (一) Maven 工程下 Spring MVC 站点配置 (二) Mybatis数据操作 这篇主 ...
- Mac系统杂项 (持续更新)
一.调整LaunchPad的图标显示列数和行数 .调整每一列显示图标数量.在我的电脑上(1366 * 768),每列7个个人觉得比较不错 defaults write com.apple.dock s ...
- Java IO之序列化
序列化机制是Java语言内建的一种对象持久化方式,可以很容易的在JVM中的活动对象和字节数组之间转换.它的一个重要用途就是远程方法调用的时候,用来对开发人员屏蔽底层实现细节(远端的开发人员不知道这个对 ...
- RecyclerView不同类型Item的展示
代码如下: public class AccessoiresAdapter extends RecyclerView.Adapter { final int VIEW_TYPE_ACCESSORY = ...
- Mongodb 安装和启动
一.首先去官网下载对应的的mongodb ,本人的操作系统是win7 64位 mongodb-win32-x86_64-2.0.6.rar 解压安装:进入到bin目录下,会看到N多的.exe文件 二. ...