当前位置: 代码迷 >> 综合 >> 洛谷 P2822 [NOIP2016 D2T1] 组合数问题 [90分伪题解]
  详细解决方案

洛谷 P2822 [NOIP2016 D2T1] 组合数问题 [90分伪题解]

热度:81   发布时间:2024-01-19 02:35:40.0

题目描述

组合数表示的是从n个物品中选出m个物品的方案数。举个例子,从(1,2,3) 三个物品中选择两个物品可以有(1,2),(1,3),(2,3)这三种选择方法。根据组合数的定 义,我们可以给出计算组合数的一般公式:

其中n! = 1 × 2 × · · · × n

小葱想知道如果给定n,m和k,对于所有的0 <= i <= n,0 <= j <= min(i,m)有多少对 (i,j)满足是k的倍数。

输入输出格式

输入格式:

第一行有两个整数t,k,其中t代表该测试点总共有多少组测试数据,k的意义见 【问题描述】。

接下来t行每行两个整数n,m,其中n,m的意义见【问题描述】。

输出格式:

t行,每行一个整数代表答案。

输入输出样例

输入样例#1:
1 2
3 3
输出样例#1:
1


输入样例#2:
2 5
4 5
6 7
输出样例#2:
0
7









说明

【样例1说明】

在所有可能的情况中,只有是2的倍数。

【子任务】

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

唯一分解定理~

(在老师×威下写了这篇题解,但真的不是正解啊,只能拿90分……)

先分解k,记录下它的每个质因数以及该质因数的次数~

然后对于每个数,计算它关于k的质因数的次数,如果每个次数都不小于k的次数,那么这个数就是k的倍数~

因为组合数公式c[i][j]=i*(i-1)*(i-2)*...*(i-j+2)*(i-j+1)/(1*2*3*...*(j-1)*j),所以如果已知c[i][j-1]的所有次数,就可以直接加i-j+1的次数再减去j的次数,顺便判断是否是倍数,如果是,加1即可,最后输出答案~

但这种解法时间复杂度很高,18和20两个点会T的!


#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;int n,m,t,k,a[1001][2],tota,num[1001],ans;void chang(int u)
{for(int i=2;i<=u;i++){if(!(u%i)){a[++tota][0]=i;while(!(u%i)) a[tota][1]++,u/=i;}}if(u>1) a[++tota][0]=u,a[tota][1]=1;
}bool findd(int u,int kkz)
{bool flag=1;for(int i=1;i<=tota;i++){while(!(u%a[i][0])) num[i]+=kkz,u/=a[i][0];if(num[i]<a[i][1]) flag=0;}return flag;
}bool find2(int u)
{int kkz=u,flag=1;for(int i=1;i<=tota;i++){num[i]=0;while(!(u%a[i][0])) num[i]++,u/=a[i][0];if(num[i]<a[i][1]) flag=0;}return flag;
}int main()
{
//	freopen("problem.in","r",stdin);
//	freopen("problem.out","w",stdout);scanf("%d%d",&t,&k);chang(k);while(t--){scanf("%d%d",&n,&m);ans=0;for(int i=1;i<=n;i++){if(find2(i)) ans++;int kkzvu=min(i,m);for(int j=2;j<=kkzvu;j++){findd(i-j+1,1);if(findd(j,-1)) ans++;}}printf("%d\n",ans);}return 0;
}