CF 900D Unusual Sequence 容斥

题意:问有多少个序列a.长度任意.序列a元素之和为y,序列a的所有元素的gcd为x. x,y<=1e9 


由题意得:x * (b[1]+b[2]+...b[k])=y. 
gcd(b[1],b[2],...b[k])=1. (b[1]+b[2]+..b[k]) =y/x (不能整除则无解.)


令n=y/x f(n)为序列gcd等于1,并且序列元素和为n的方法数.
序列和要想为n,现在序列的个数任意
因为n个1有n-1个间隔,每个间隔可以选择隔离或者不隔离 共2^(n-1)种方案 但是不保证gcd为1.扣掉gcd>1的方案.


f(n)=2^(n-1) - ∑f(ti) ti为n的因子.
f(ti)=2^(ti-1) -∑ f(pi) pi为ti的因子,同时也肯定是n的因子.

处理出n的所有因子,从小到大计算f(ti),按上面式子扣掉重复即可.O(p^2) p为1e9内最多因子个数

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=2e5+5,mod=1e9+7;
ll x,y,n,d[N];
vector<int> v;
ll powmod(ll x,ll n)
{
	ll s=1;
	while(n)
	{
		if(n&1)
			s=(s*x)%mod;
		x=(x*x)%mod;
		n>>=1;
	}
	return s;
}
int main()
{
	cin>>x>>y;
	if(y%x)
	{
		puts("0");
		return 0;
	}
	n=y/x;
	for(int i=1;i*i<=n;i++)
	{
		if(n%i==0)
		{
			v.push_back(i);
			if(i*i!=n)
				v.push_back(n/i);
		}
	}
	sort(v.begin(),v.end());
	for(int i=0;i<v.size();i++)
	{
		d[i]=powmod(2ll,v[i]-1);
		for(int j=0;j<i;j++)
		{
			if(v[i]%v[j]==0)
				d[i]=(d[i]-d[j]+mod)%mod;
		}
	}
	cout<<d[v.size()-1]<<endl;
	return 0;
}