Easy Math（ACM-ICPC 2018 徐州赛区网络预赛）（递归 + 杜教筛）

Easy Math

推式子

∑i=1mμ(in)∑i=1mμ(i×nd×d)，d是n的一个质因子i,d互质项有(−∑i=1mμ(i×nd))，由于减去了多余的非互质项，所以加上，−∑i=1mμ(i×nd)+∑i=1mdμ(i×d×nd)−∑i=1mμ(i×nd)+∑i=1mdμ(i×n)\sum_{i = 1} ^{m} \mu(in)\\ \sum_{i = 1} ^{m} \mu(i \times \frac{n}{d} \times d)，d是n的一个质因子\\ i, d互质项有(-\sum_{i = 1} ^{m} \mu(i \times \frac{n}{d}))，由于减去了多余的非互质项，所以加上，\\ -\sum_{i = 1} ^{m} \mu(i \times \frac{n}{d}) + \sum_{i = 1} ^{\frac{m}{d}} \mu(i \times d \times \frac{n}{d})\\ -\sum_{i = 1} ^{m} \mu(i \times \frac{n}{d}) + \sum_{i = 1} ^{\frac{m}{d}} \mu(i \times n)\\ i=1∑mμ(in)i=1∑mμ(i×dn×d)，d是n的一个质因子i,d互质项有(−i=1∑mμ(i×dn))，由于减去了多余的非互质项，所以加上，−i=1∑mμ(i×dn)+i=1∑dmμ(i×d×dn)−i=1∑mμ(i×dn)+i=1∑dmμ(i×n)

由此我们可以递归加上杜教筛来得到答案：

边界条件n==1n == 1n==1，return∑i=1mμ(i)return\ \sum_{i = 1} ^{m} \mu(i)return ∑i=1mμ(i)。

m==0,return0m == 0, return\ 0m==0,return 0。

代码

/*Author : lifehappy
*/
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include <bits/stdc++.h>using namespace std;typedef long long ll;
const int inf = 0x3f3f3f3f;const int N = 2e6 + 10;int prime[N], mu[N], cnt;bool st[N];void init() {mu[1] = 1;for(int i = 2; i < N; i++) {if(!st[i]) {prime[cnt++] = i;mu[i] = -1;}for(int j = 0; j < cnt && 1ll * i * prime[j] < N; j++) {st[i * prime[j]] = 1;if(i % prime[j] == 0) break;mu[i * prime[j]] = -mu[i];}}for(int i = 1; i < N; i++) {mu[i] += mu[i - 1];}
}map<ll, ll> ans_s;ll S(ll n) {if(n < N) return mu[n];if(ans_s.count(n)) return ans_s[n];ll ans = 1;for(ll l = 2, r; l <= n; l = r + 1) {r = n / (n / l);ans -= (r - l + 1) * S(n / l);}return ans_s[n] = ans;
}ll F(ll n, ll m) {if(n == 1) return S(m);if(m == 0) return 0;for(int i = 0; 1ll * prime[i] * prime[i] <= n; i++) {if(n % prime[i] == 0) {return F(n, m / prime[i]) - F(n / prime[i], m);} }return F(n, m / n) - F(n / n, m);
}int main() {//    freopen("in.txt", "r", stdin);
//    freopen("out.txt", "w", stdout);
//    ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);init();ll n, m;scanf("%lld %lld", &m, &n);ll temp = n; for(int i = 0; i < cnt; i++) {int cnt = 0;while(n % prime[i] == 0) {n /= prime[i];cnt++;}if(cnt >= 2) {puts("0");return 0;}}printf("%lld\n", F(temp, m));return 0;
}

Easy Math（ACM-ICPC 2018 徐州赛区网络预赛）（递归 + 杜教筛）相关推荐

ACM-ICPC 2018 徐州赛区网络预赛 D. Easy Math
Easy Math 问答问题反馈只看题面 16.47% 1000ms 262144K Given a positive integers nn , Mobius function \mu(n)μ(n ...
ACM-ICPC 2018 徐州赛区网络预赛 Features Track（STL二维map）
Morgana is learning computer vision, and he likes cats, too. One day he wants to find the cat moveme ...
ACM-ICPC 2018 徐州赛区网络预赛 D. EasyMath
ACM-ICPC 2018 徐州赛区网络预赛 D. EasyMath 做法: \[f(m,n) = \sum _{i=1}^{m} \mu(in) = \sum_{i=1}^{m}[gcd(i,n)= ...
ACM-ICPC 2018 徐州赛区网络预赛G （单调队列）
传送门题面: There's a beach in the first quadrant. And from time to time, there are sea waves. A wave ( ...
ACM-ICPC 2018 徐州赛区网络预赛 - A. Hard to prepare - （计数递归）
题目链接 After Incident, a feast is usually held in Hakurei Shrine. This time Reimu asked Kokoro to deli ...
ACM-ICPC 2018 徐州赛区网络预赛　I. query　树状数组
I. query 题目链接: Problem Description Given a permutation \(p\) of length \(n\), you are asked to answe ...
ACM-ICPC 2018 徐州赛区网络预赛
BE, GE or NE 题意: 每一轮有三种操作, 加上a 减去b 或者取负当且仅当 a, b, c 不为0时,对应的操作有效; 给出一个上界和一个下界大于等于上界就是 Good Ending ...
ACM-ICPC 2018 徐州赛区网络预赛 Morgana Net
题意:Morgana Net 题解:把a矩阵每一位根据公式推出递推矩阵,然后用矩阵快速幂,比赛没想到啊,, #include<bits/stdc++.h> #define ll long ...
ACM-ICPC 2018 徐州赛区网络预赛 I. Characters with Hash
注意几种情况,前导零是不要的:如果在开头的第一个值为个位数,那么占一位,如果个位数出现在中间,那么占两位,比如6就是06 #include <iostream> #include < ...

Easy Math（ACM-ICPC 2018 徐州赛区网络预赛）（递归 + 杜教筛）

Easy Math

推式子

代码

Easy Math（ACM-ICPC 2018 徐州赛区网络预赛）（递归 + 杜教筛）相关推荐

最新文章

热门文章