hdu 4279 Number 2012天津网络赛 数论 分析

题目

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4279

题目来源：2012天津网络赛，其实不简单，但是由于是B华丽丽地成了前三题。

简要题意：$f(n)=| \{x\mid \gcd(x,n)\neq1 \wedge x \not\mid n \wedge x\in [1,n]\}|$
$\qquad \quad \ \ \ \ $即$n$内与$n$不互质且不能整除$n$的数的个数。
$\qquad \quad \ \ \ \ $求$[x,y]$内$f(n)$为奇数的数的个数。

数据范围：$T\leqslant 2000;\quad 0<x,y<2^{63}$

题解

在我看来，这是一道很好的数论题目，运用到了不少的知识，但是我第一次看到的时候没能写出来。

定义$\tau(x)$为$x$的正除数个数，$\varphi(x)$为欧拉函数，即小于等于$x$与$x$互质的数的个数。

则易知$f(n)=n-\varphi(n)-\tau(n)+1$，$1$是约数也和$n$互质，要再加上。

首先需要用到的结论是$2\mid \varphi(x)\quad(x>2)$，即大于2的数的欧拉函数值为偶数。

证明：$2\mid \varphi(x)\quad(x>2)$

$$x=p_1^{a_1}p_2^{a_2}\cdots p_k^{a_k};\quad \varphi(x)=\prod_{i=1}^k{(p_i-1)}p_i^{a_1-1}\qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad$$ 若 $k=1,p_1=2$则 $a_1>1,2\mid p_1^{a_1-1},2\mid \varphi(x)$

因为$2$以外正质数都是奇数，所以除此之外的情况必有$2\not\mid p_i$，则$2\mid (p_i-1),2\mid \varphi(x)$

因而$2\mid \varphi(x)\quad(x>2)$

$f(n)\quad(n>2)$为奇数的充要条件变为$n+1-\tau(n)$为奇数，$f(1),f(2)$易知是不满足条件的。

对于奇数$n$我们需要$\tau(n)$为偶数。
对于偶数$n$我们需要$\tau(n)$为奇数。

沿用之前$x,p_i,a_i$我们可以知道有以下公式：

$$\tau(x)=\prod_{i=1}^k(a_i+1)\qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad$$ $\tau(x)$为奇数当且仅当任意 $a_i$为偶数，显然此时的 $x$一定是个完全平方数。

对于奇数$n$我们需要 $n$为完全平方数。
对于偶数$n$我们需要 $n$不为完全平方数。

考虑数的平方和自己奇偶性一致，只知道$\sqrt{n}$内奇数偶数的个数就能知道答案。

结果就是$(n$内偶数的个数$-1)-($$\sqrt{n}$内偶数的个数$)+(\sqrt{n}$内奇数的个数$-1)$需要特判$1,2$。

让结果为$Ans$，可得最终公式为：

$$\begin{equation} Ans= \left\{ \begin{aligned} &0~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~n\in\{1,2\}\\ &\lfloor n/2\rfloor-\lfloor \sqrt{n}/2\rfloor+\lfloor (\sqrt{n}-1)/2\rfloor~~~~~~~else \end{aligned} \right. \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \end{equation}$$

实现

$\quad\ \ $结论出来就很好搞了，但是取根号这个地方有一些精度问题，可以用long double或者是用二分，训练的时候就被精度卡了。复杂度$\Theta(\log n)$就是开方的复杂度。

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <stack>
#include <queue>
#include <string>
#include <vector>
#include <set>
#include <map>

#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define sz(x) ((int)(x).size())
#define fi first
#define se second
using namespace std;
typedef long long LL;
typedef vector<int> VI;
typedef pair<int,int> PII;
LL powmod(LL a,LL b, LL MOD) {LL res=1;a%=MOD;for(;b;b>>=1){if(b&1)res=res*a%MOD;a=a*a%MOD;}return res;}
// head

LL query(LL x) {
    LL sq = sqrt((long double)x);
    LL ans = x > 2 ? x/2-1 : 0;
    ans -= sq/2;
    ans += (sq-1)/2;
    return ans;
}
int main()
{
    int t;
    LL l, r;
    scanf("%d", &t);
    while (t--) {
        scanf("%I64d%I64d", &l, &r);
        printf("%I64d\n", query(r)-query(l-1));
    }
    return 0;
}