codeforces 900D 数论+组合+容斥原理

问有多少个这样的数字序列

所有数的GCD等于x 并且 所有数的和等于y

题解：

非常难有思路啊 看题解后过的。

考虑序列GCD为x的倍数 即GCD = n*x 和当然都为y 这个条件不要忘了

这样我们可以用  容斥原理来递推的计算GCD为n*x的序列个数是多少

怎么计算呢

以样例为例子 3 9

当GCD = 3 的时候 可以有9 / 3 = 3 个3 序列是这样的 3 3 3

那么有三个空 用插板法 可以计算可以插板的方式数位2**(3-1) = 2**2 = 4种

这里解释插板的意义 3|3 3插一个板就表示相邻的数求和 那么3|3 3 就是 6 3

同理 3 3|3 -> 3 6; 3|3|3 -> 9； 但是这样插板出现了问题 就是出现了GCD 并不为3的序列 即 9

这个时候就需要用容斥原理来 递推

设a[i] 表示GCD为i的序列个数， a[j] 表示GCD为j个数

不妨设i > j

if (i % j == 0) a[j] = (a[j]+MOD-a[i]) % MOD; 因为 GCD 为i的情况是一定可以通过插板 得到GCD 为j的情况

而通过从大到小的递推 a[i]已经是容斥后的结果

这样最终得到GCD最小为x 的结果就是答案 网上题解写的有点不太清楚 这里自己补充点自己的理解。轻喷。。

代码君：

 #include <bits/stdc++.h>
 #include <string.h>
 #include <iostream>
 #include <stdio.h>
 #define pb push_back

 const int MOD = 1e9+;
 const int MAXN = 1e5+;
 typedef long long ll;

 using namespace std;

 ll x, y;

 ll mypow(ll base, ll p)
 {
     if (p == ) return ;
     ll tmp = mypow(base, p/);
     if (p & ) tmp = (tmp*tmp*base) % MOD;
     else tmp = (tmp*tmp) % MOD;
     return tmp;
 }
 vector<ll> a;
 ll dp[MAXN];
 int main()
 {
     //freopen("in.txt", "r", stdin);
     while (cin >> x >> y)
     {
         a.clear();
         if (y % x != )
         {
             cout <<  << endl;
             continue;
         }
         for (ll i = ; i*i <= y; i++)
         {
             if (i % x ==  && y % i == ) a.pb(i);
             if (i*i != y && y % i ==  && (y/i)%x== ) a.pb(y/i);
         }
         sort(a.begin(), a.end());
         for (int i = ; i < a.size(); i++) dp[i] = mypow(, (y/a[i]-));
         for (int i = (int)a.size()-; i >= ; i--)
             for (int j = i+; j < a.size(); j++)
                 if (a[j] % a[i] == )
                 {
                     dp[i] -= dp[j];
                     dp[i] = (dp[i] + MOD) % MOD;
                 }
         cout << dp[] << endl;
     }
     return ;
 }