2021ICPC网络赛第二场 :L Euler Function 势能线段树

最新推荐文章于 2024-09-14 11:35:12 发布

原创

最新推荐文章于 2024-09-14 11:35:12 发布 · 290 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#算法

分析

首先考虑欧拉函数的性质

if(p是质数)
	{
   
   
		if(i % p == 0) f[i * p] = f[i] * p;
		else f[i * p] = f[i] * (p - 1);
	}

因为 $100$ 以内的质数很少，所以考虑维护 $b i t s e t$ 判断质因子是否出现过即可，然后用势能线段树维护即可

代码

#pragma GCC optimize(3)
#include <bits/stdc++.h>
#define debug(x) cout<<#x<<":"<<x<<endl;
#define dl(x) printf("%lld\n",x);
#define di(x) printf("%d\n",x);
#define _CRT_SECURE_NO_WARNINGS
#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define SZ(x) ((int)(x).size())
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
typedef vector<int> VI;
const int INF = 0x3f3f3f3f;
const int N = 2e5 + 10;
const ll mod = 998244353;
const double eps = 1e-9;
const double PI = acos(-1);
template<typename T>inline void read(T