自己的第一次容斥原理的实践（还有等差数列求和）终于会做数学题了

zhhx2001

于 2016-05-07 14:20:30 发布

阅读量439

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/zhhx2001/article/details/51337654

版权

数论专栏收录该内容

42 篇文章

订阅专栏

本文展示了一个使用C++编写的程序，该程序通过子集枚举的方法实现了等差数列求和，并运用了容斥原理来解决特定问题。代码中使用了标准库中的各种功能，包括输入输出、算法处理和数学运算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<string>
#define ll long long
#define mod 376544743
using namespace std;

ll sn(ll a1,ll an,ll n)//等差数列求和
{
	return (a1+an)%mod*n%mod*188272372%mod; //当时忘了搞2对于mod的乘法逆元了，然后就wrong了
}
int count(int n)
{
	int ans;
	for (ans=0;n;n=n&(n-1)) ans++;
	return ans;
}
ll n,m,p[50],ans;
int main()
{
	scanf("%lld%lld",&n,&m);
	for (int i=1;i<=n;i++) scanf("%lld",&p[i]);
	for (int i=1;i<(1<<n);i++)//枚举子集
	{
		int j=i;
		ll k=count(j),g=0,t=1;
		for (int l=0;l<n;l++) if (j&(1<<l)) t*=p[l+1]; 
		if (t>m) continue;
		g=sn(t,(m/t)*t,m/t);
		if (k&1) ans=(ans+g)%mod; //容斥原理
		else ans=(ans+mod-g)%mod;	 
	}
	printf("%lld",ans);
	return 0;
}