WC模拟(1.10) T1 Function

标签: NOI 数论 莫比乌斯反演
133人阅读 评论(0) 收藏 举报
分类:

Function

题目背景:

1.10 WC模拟T1

分析:数论 + 莫比乌斯反演

 

将这种丧病的题出到T1,也就这个出题人想的出来······这道题的本质其实是spojdivcnt2,现在考虑直接推式子······


w(d)表示d的质因数个数,这一步的转化我们这么考虑,对于i2的约数,可以被分为两类,是i的约数也是i2的约数,是i的约数,但不是i的约数,那么我们考虑对于i的一个约数x,显然,x的每一个质因数的次幂,一定小于等于i的对应的质因数的次幂,那么我们对于x的某一个质因数的次幂加上i的对应质因数的次幂,即:pb能整除xpa能整除i,那么对于x * pa一定能够对应一个i2的约数,并且这个数一定不会是i的约数,考虑这意味着什么,我们可以对于x任意包含的一个质因数集合,加上对应的在i中的幂次,一定都可以得到一个i2的约数,并且不是i的约数,但是这个质因数x一定要原本包含,那么也就是x的可行的组成集合的元素为w(d)个,那么,可行的选择就有2w(d)种。我们可以发现:


考虑2w(d)的本质,相当于选择一个d的质因数集合,也就是d的不含平方质因子的因数的个数,考虑这类数的相同点就是莫比乌斯函数的值不为0,那么我们直接加莫比乌斯函数值的平方就好了,那么于是就化成了



g为右半边的式子,我们可以发现g(n)相当于1 ~ n的约数个数之和,那么我们可以通过线筛求出,也可以用sqrt(n)的时间直接计算,同样的对于左半边的式子相当于是莫比乌斯函数的平方的前缀和,考虑如何快速求解,我们可以发现原式就是求1 ~ n中不含有平方因子的数的个数,那么我们考虑容斥一下,用所有数先减去一个质数的平方的的倍数,再加上两个质数之积的平方的倍数,减去三个质数之积的平方的倍数······这也就相当于利用莫比乌斯函数容斥一下了所以我们可以得到


所以我们同样可以结合预处理和根号时间计算,考虑类似杜教筛的复杂度分析,我们可以知道预处理n2/3是最优的那么最后的复杂度就是O(n2/3)

 

Source:

/*
	created by scarlyw
*/
#include <cstdio>
#include <string>
#include <algorithm>
#include <cstring>
#include <iostream>
#include <cmath>
#include <cctype>
#include <vector>
#include <set>
#include <queue>
#include <ctime>
#include <bitset>

inline char read() {
	static const int IN_LEN = 1024 * 1024;
	static char buf[IN_LEN], *s, *t;
	if (s == t) {
		t = (s = buf) + fread(buf, 1, IN_LEN, stdin);
		if (s == t) return -1;
	}
	return *s++;
}

///*
template<class T>
inline void R(T &x) {
	static char c;
	static bool iosig;
	for (c = read(), iosig = false; !isdigit(c); c = read()) {
		if (c == -1) return ;
		if (c == '-') iosig = true;	
	}
	for (x = 0; isdigit(c); c = read()) 
		x = ((x << 2) + x << 1) + (c ^ '0');
	if (iosig) x = -x;
}
//*/

const int OUT_LEN = 1024 * 1024;
char obuf[OUT_LEN], *oh = obuf;
inline void write_char(char c) {
	if (oh == obuf + OUT_LEN) fwrite(obuf, 1, OUT_LEN, stdout), oh = obuf;
	*oh++ = c;
}

template<class T>
inline void W(T x) {
	static int buf[30], cnt;
	if (x == 0) write_char('0');
	else {
		if (x < 0) write_char('-'), x = -x;
		for (cnt = 0; x; x /= 10) buf[++cnt] = x % 10 + 48;
		while (cnt) write_char(buf[cnt--]);
	}
}

inline void flush() {
	fwrite(obuf, 1, oh - obuf, stdout);
}

/*
template<class T>
inline void R(T &x) {
	static char c;
	static bool iosig;
	for (c = getchar(), iosig = false; !isdigit(c); c = getchar())
		if (c == '-') iosig = true;	
	for (x = 0; isdigit(c); c = getchar()) 
		x = ((x << 2) + x << 1) + (c ^ '0');
	if (iosig) x = -x;
}
//*/

const int MAXN = 1000000 + 10;

int  prime_cnt, t, cur;
int prime[MAXN], min_cnt[MAXN];
long long miu[MAXN], sum[MAXN], f[MAXN];
bool not_prime[MAXN];

inline void seive() {
	not_prime[1] = true, miu[1] = 1, sum[1] = 1;
	for (int i = 2; i < MAXN; ++i) {
		if (!not_prime[i]) {
			prime[++prime_cnt] = i, miu[i] = -1;
			sum[i] = 2, min_cnt[i] = 1;
		}
		for (int j = 1; j <= prime_cnt && prime[j] * i < MAXN; ++j) {
			int cur = prime[j] * i;
			if (i % prime[j]) {
				miu[cur] = -miu[i], not_prime[cur] = true;
				sum[cur] = sum[i] * 2, min_cnt[cur] = 1;
			} else {
				miu[cur] = 0, not_prime[cur] = true;
				min_cnt[cur] = min_cnt[i] + 1;
				sum[cur] = sum[i] / min_cnt[cur] * (min_cnt[cur] + 1);
				break ;
			}
		}
	}
	for (int i = 1; i < MAXN; ++i) 
		f[i] += f[i - 1] + miu[i] * miu[i], sum[i] += sum[i - 1];
}

inline long long get_miu(int cur) {
	if (cur < MAXN) return f[cur];
	long long ans = 0;
	for (int i = 1, end = sqrt(cur); i <= end; ++i)
		ans += miu[i] * (cur / i / i);
	return ans;
}

inline long long get_sum(long long cur) {
	if (cur < MAXN) return sum[cur];
	long long ans = 0;
	for (int i = 1, last; i <= cur; i = last + 1) {
		last = cur / (cur / i);
		ans += (long long)(last - i + 1) * (long long)(cur / i);
	}
	return ans;
}

inline void solve(long long cur) {
	long long ans = 0;
	for (int i = 1, last; i <= cur; i = last + 1) {
		last = cur / (cur / i);
		ans += (get_miu(last) - get_miu(i - 1)) * get_sum(cur / i);
	}
	W(ans), write_char('\n');
}

int main() {
	freopen("function.in", "r", stdin);
	freopen("function.out", "w", stdout);
	R(t), seive();
	while (t--) R(cur), solve(cur);
	flush();
	return 0;
}


查看评论

ATM取款机模拟实战项目

我们从现实生活中的ATM取款机的广泛使用,来确定我们编写ATM取款机程序的功能。首先,用户的登陆、用户的取款、用户的查询、用户的转帐、用户密码的修改等功能。根据实际功能的需求完成我们的程序。 “ATM取款机项目”主要描述“ATM取款机”的工作原理,与实际生活相结合,编写一个模拟实战项目。
  • 2015年05月12日 14:29

WC模拟(1.6) T1 送你一个DAG

斯特林数 + DP
  • scar_lyw
  • scar_lyw
  • 2018-01-11 15:45:26
  • 185

WC模拟(1.10) T2 Permutation

结论推理
  • scar_lyw
  • scar_lyw
  • 2018-01-17 20:23:25
  • 105

WC模拟(1.10) T3 Subsequence

DP + 复杂度分析
  • scar_lyw
  • scar_lyw
  • 2018-01-17 20:23:48
  • 111

WC2017总结

~~~~~~吾之于WC,如蛙之于井底,仰而望天之高,跃而不及地。Part 1突然开始的WC模拟      ~~~~~~noip2016在11月20日结束,经过一段时间的乱搞,我们于12月开始,为期3个...
  • rzO_KQP_Orz
  • rzO_KQP_Orz
  • 2017-02-09 20:17:14
  • 4112

WC模拟(1.4) T1 序列

迭代加深搜索
  • scar_lyw
  • scar_lyw
  • 2018-01-11 14:56:51
  • 119

WC2018口胡记

作为一个noip考挂选手,我并不能去WC2018。 当天下午一点半左右问@tututu要来题目,一开始只有T1的复述版本,感觉是个数据结构题,但并没有什么思路,而且部分分表也没有,似乎只会O(n2)...
  • zxin__
  • zxin__
  • 2018-02-10 23:13:27
  • 311

WC模拟(1.2) T1 串

结论分析
  • scar_lyw
  • scar_lyw
  • 2018-01-03 14:08:37
  • 113

vb 自动识别拼音

  • 2011年12月18日 15:14
  • 11KB
  • 下载

2017.1.17WC模拟赛总结

boom zero了,很不开心,总结一下问题。 T1boom zero了,原因是常数太大自己又没有去卡。 T2boom zero了,原因是没看好题,不符合它输出要求。 T3boom zero了,...
  • WerKeyTom_FTD
  • WerKeyTom_FTD
  • 2017-01-17 14:47:27
  • 513
    个人资料
    专栏达人 持之以恒
    等级:
    访问量: 2万+
    积分: 2171
    排名: 2万+
    博客专栏
    最新评论