容斥原理 - 分糖（SOJ 747）

最新推荐文章于 2023-06-02 20:23:08 发布

Faithfully__xly

最新推荐文章于 2023-06-02 20:23:08 发布

阅读量389

点赞数

分类专栏： --------数学知识-------- 组合数学文章标签：容斥原理

本文链接：https://blog.csdn.net/weixin_42557561/article/details/100046173

版权

这是一篇关于利用容斥原理解决分糖问题的文章。题目要求在N个糖果和M个小朋友之间分配，确保每个小朋友至少得到一个糖果，且糖果数不能是大于1的正整数的倍数，也不能有剩余。输入包含数据组数和每组的N、M值，输出符合条件的分配方法总数模1000000007。分析中提到，可以通过求解N的质因数并进行容斥计算来得出答案。解决方案涉及到组合计数的问题，即计算非空集合的组合数。

摘要由CSDN通过智能技术生成

分糖

题目描述

有 N 个（相同的）糖果，M 个（不同的）小朋友。M 和 N 满足：1≤M≤N≤100000(105)。
要求：
1.每个小朋友都至少有一个糖果。
2.不存在正整数 X（X>=2），使得每个小朋友的糖果数都是 X 的倍数。
3.糖果不能剩余。
求分糖方法总数。答案模 1000000007（109+7）

输入格式

第一行为数据组数：T<=100000。
接下来 N 行，每行 2 个如上文所示的正整数 N，M。

输出格式

输出 T 行，每行一个整数，为答案。
注意取模！

Analysis

难度不是很大的一道题啦：）
正难则反
我们统计每个小朋友的糖果数都是 X 的倍数的情况
显然X应该是N的因数
那么问题就变成了将N分解质因数
然后容斥一下
方案数的话就是将a个相同物品放入b个不同集合，集合非空的情况
$ans=C_{a-1}^{b-1}$

Code

#include<bits/stdc++.h>
#define int long long
#define re register
#define in read()
using namespace std;
inline int read(){
   
	char ch;int f=1,res=0;
	while((ch=getchar())<'0'||ch>'9') if(ch=='-') f=-1;
	while(ch>