2022CCPC网络赛

Mutiple Set-dfs求因子+组合数

题意：
就是给你一段区间a到b，和权值c，然后问你有多少x，使得区间[a,b]中所有是x的倍数的数组成了一个集合set，然后这个集合的所有子集的价值总和 = c。一个集合的价值就是集合内所有数的和。然后T组测试，T=100,1 ≤ L ≤ R ≤ 1e12，1<=C<=1e14。

思考：

1. 当时这题过的不算太多，我看了遍题意，我以为就是x的所有倍数的总和 = c，没看到子集。然后队友提醒后，那么稍微写了个样例，由于涉及到子集，那么直接看某个数的贡献可以了，比如x的倍数有sum个，那么一个数的贡献就是2^sum-1 组合数定理。
2. 然后既然所有x的倍数组成的和 = c,那么c肯定是x的倍数，那么也就是要枚举c的因子。如果暴力复杂度1e9过不了。但是如果使用dfs求因子，这样就够快了。对于每个因子判断一下是否2^sum-1>c,以及[a,b]这段区间x的倍数的总和是否>c，由于总和可能很大，所以先用int128来判断一下不可能的情况。
3. 如何求一段区间x的倍数的总和呢？首先对于最左边可是x的多少倍，最右边可以是x的多少倍，左边就是a/x上取整，右边就是b/x下取整。然后总和就是x*(一个公差为1的等差数列) = x*(r-l+1)*(l+r)/2。不过由于左边是上取整，有可能右边是0，左边为1，这样就不可以了，要特判掉。然后就判断判断C是不是这两个数的乘积就可以了，不过直接×可能会炸，所以用(C%sum1==0&&C/sum1 == sum2)。
4. 不过要注意题目说的各种情况，当时队友写的差不多了，但是忘记判断No Solution和too many了，所以不要慌，一点一点理清思路。

代码：

#include<bits/stdc++.h>
#define fi first
#define se second
#define pb push_back
#define db double
#define int long long
#define PII pair<int,int >
#define mem(a,b) memset(a,b,sizeof(a))
#define IOS std::ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);

using namespace std;
const int mod = 1e9+7,inf = 1e18;
const int N = 2e5+10,M = 2e6+10;

int T,n,m,k;
PII va[M];
int vb[M];
int pri[M],st[M],pre[M];
int cnt,cnt1,cnt2;

void init(int x)
{
	pre[0] = 1;
	for(int i=1;i<=60;i++) pre[i] = pre[i-1]*2;
	for(int i=2;i<=x;i++)
	{
		if(!st[i]) pri[++cnt] = i;
		for(int j=1;pri[j]*i<=x;j++)
		{
			st[pri[j]*i] = 1;
			if(i%pri[j]==0) break;
		}
	}
}

void dfs(int now,int sum)
{
	if(now>cnt1)
	{
		vb[++cnt2] = sum;
		return ;
	}
	int res = 1;
	for(int i=0;i<=va[now].se;i++)
	{
		dfs(now+1,sum*res);
		res *= va[now].fi;
	}
}

void solve(int a,int b,int c,int cs)
{
	int tmp = c;
	for(int i=1;pri[i]<=tmp/pri[i];i++)
	{
		if(tmp%pri[i]==0)
		{
			int res = 0;
			while(tmp%pri[i]==0) res++,tmp /= pri[i];
			va[++cnt1] = {pri[i],res};
		}
	}
	if(tmp>1) va[++cnt1] = {tmp,1};
	dfs(1,1);
	vector<int > anw;
	for(int i=1;i<=cnt2;i++)
	{
		int x = vb[i];
		int l = (a+x-1)/x,r = b/x;
		__int128 sumt = (r-l+1)*(l+r)/2;
		if(r-l-1>50||sumt>c||sumt*x>c||sumt==0) continue;
		int sum1 = x*(r-l+1)*(l+r)/2;
		int sum2 = pre[r-l+1-1];
		if(c%sum1==0&&c/sum1==sum2) anw.pb(x);
	}
	if(anw.size()==0) cout<<"No Solution";
	else if(anw.size()>1e5) cout<<"Too Many!";
	else
	{
		cout<<anw.size()<<"\n";
		sort(anw.begin(),anw.end());
		for(int i=0;i<anw.size();i++)
		{
			cout<<anw[i];
			if(i!=anw.size()-1) cout<<" ";
		}
	}
	if(cs) cout<<"\n";
}

signed main()
{
	IOS;
	init(2e6+5);
	cin>>T;
	while(T--)
	{
		int a,b,c;
		cin>>a>>b>>c;
		cnt1 = 0,cnt2 = 0;
		solve(a,b,c,T);
	}
	return 0;
}

总结：
多多思考，注重细节，注意输出。