题目:点击打开链接
题意:在不大于s的所有数中,找出k个gcd不为1的数,求组数。
分析:因为k>=2,所以必定选择两个数不互质,最大的质因数不超过50/2=25,因此枚举不超过25的质因数,对于每个质因数x,与之不互质的数共有s/x,而要从这s/x个数中选出k个数,用到组合数公式,组合数公式可以通过递推预处理。但是这样取会有重复,比如若s=20,k=3,x=2时,不互质的数有2,4,6,8,10,12,14,16,18,20,x=3时,不互质的数有3,6,9,12,15,18,其中6,12,18被重复计算了,显然,被重复计算的组合的质因数6是2和3的最小公倍数,因为这些数是质因数,不互质,所以它们的最小公倍数就是它们之积。运用容斥原理,枚举质因数,将答案减去两个质因数相乘的组数,加上三个质因数相乘的组数...直到枚举所有s/x>=k的质因数,即得答案。
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <stdio.h>
#include <iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
#include <set>
#include <vector>
#include <map>
#define sqr(x) ((x)*(x))
#define PR pair<int,int>
#define MP make_pair
#define fi first
#define se second
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define ll unsigned long long
const ll INF = 1e18;
const int inf=0x3f3f3f3f;
const int M=100010;
const int N=10010;
const ll MOD=1000000007;
const double eps=1e-3;
const double pi=acos(-1.0);
using namespace std;
int k,s;
int prime[9]={2,3,5,7,11,13,17,19,23};
int ans,C[26][26];
void init()
{
C[1][1]=C[1][0]=1;
for(int i=2;i<26;i++)
{
C[i][i]=C[i][0]=1;
for(int j=1;j<i;j++)
C[i][j]=(C[i-1][j-1]+C[i-1][j]);
}
}
int dfs(int idx,int num)
{
int res=0;
for(int i=idx;i<9;i++)
{
if(s/(num*prime[i])>=k)
res+=C[s/(num*prime[i])][k]-dfs(i+1,num*prime[i]);
}
return res;
}
int main()
{
init();
scanf("%d%d",&k,&s);
ans=min(dfs(0,1),10000);
printf("%d\n",ans);
}