URAL 1091

TeJoy

于 2020-07-01 11:14:01 发布

阅读量108

点赞数

分类专栏：容斥原理

本文链接：https://blog.csdn.net/weixin_45539557/article/details/107059331

版权

容斥原理专栏收录该内容

9 篇文章 0 订阅

订阅专栏

题意：
从 s 个数里挑出 k 个, 这个 k 个数的共因子大于 1 一共有多少种方法

思路：
很裸的容斥，首先只选择小于k的质因子i，个数显然是C（s/i,k)，由k>=2知道，i最多枚举到s/2,所以预处理组合数，跑一次容斥即可

#include<iostream>
#include<cstdio>
#include<cstring>

using namespace std;
typedef long long ll;
int c[30][30],ans,cnt;
int k,s;
int a[]={ 0,2,3,5,7,11,13,17,19,23};
void dfs(int pos,int sum,int li)
{
    sum=a[pos]*sum;
    if(s/sum<k)return;
    if(li&1)
        ans+=(c[s/sum][k]);
    else
        ans-=(c[s/sum][k]);
    if(ans>10000)return;
    for(int i=pos+1;i<=cnt;++i)
        dfs(i,sum,li+1);
}
int main()
{
   c[0][0]=1;
   for(int i=1;i<=26;++i)
   {
      c[i][0]=c[i][i]=1;
      for(int j=1;j<i;++j)	
         c[i][j]=c[i-1][j-1]+c[i-1][j]; 
   } 
   scanf("%d%d",&k,&s);
   for(int i=9;i>=1;--i)
   {
       if(s/a[i]>=k){ 
            cnt=i;
            break;
       }
   } 
   for(int i=1;i<=cnt;++i)
       dfs(i,1,1);
   if(ans>10000)
       puts("10000");
   else 
       printf("%d\n",ans);
   return 0; 
}