URAL 1091 Tmutarakan Exams 【容斥】【组合数】

题目：点击打开链接

题意：在不大于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);
}