当前位置: 代码迷 >> 综合 >> 【URAL-1091-Tmutarakan Exams】 容斥原理+组合数
  详细解决方案

【URAL-1091-Tmutarakan Exams】 容斥原理+组合数

热度:34   发布时间:2023-12-29 02:40:42.0

URAL-1091-Tmutarakan Exams
将题意简化一下就是给你一个数s,要用1-s构造一个长度为序列保证n个数的最大公约数为1,且n个数均不相同。
我们直接构造不太好想,于是我们想如何构造出不满足情况的序列,那么就是最大公约数>=2.由于这道题s,k均<=50.所以我们可以枚举gcd为2-S,同样如果一个序列的gcd为4,那么他一定被gcd为2的计算过,所以我们还是只需要枚举所有质数的组合即可,当枚举gcd为2时,一共n个数,一共s/2种选择,那么C(s/2,n)即是最终答案,以此类推,容斥一下就解决了。

#include<stdio.h>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn = 105;
long long c[maxn][maxn];
int pri[20]={
   2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,};
void getC()
{c[0][0]=1;for(int i=1;i<=100;i++){c[i][0]=c[i][i]=1;for(int j=0;j<i;j++)  c[i][j]=(c[i-1][j]+c[i-1][j-1]);}
}
int main()
{int k,s;scanf("%d%d",&k,&s);getC();ll ans=0;for(int i=1;i<(1<<15);i++){ll num=1;int cnt=0;for(int j=0;j<15;j++){if(i&(1<<j)){cnt++;num*=pri[j];}}ll tmp=0;if(s/num>=k) tmp=c[s/num][k];if(cnt&1) ans+=tmp;else ans-=tmp;}printf("%lld\n",min(ans,(ll)10000));return 0;
}