HDU 5597 GTW likes function

题意：

现在给出下列两个定义：
	f(x)=f_{0}(x)=\sum_{k=0}^{x}(-1)^{k}2^{2x-2k}C_{2x-k+1}^{k},f_{n}(x)=f(f_{n-1}(x))(n\geq 1)f(x)=f​0​​(x)=∑​k=0​x​​(−1)​k​​2​2x−2k​​C​2x−k+1​k​​,f​n​​(x)=f(f​n−1​​(x))(n≥1)

	$\phi (n)$为欧拉函数。指的是不超过$n$的与$n$互质的正整数个数。

对于每组数据，GTW有两个正整数$n,x$，现在他想知道函数\varphi(f_{n}(x))φ(f​n​​(x))的值。

其实就是喊我们先去找这个函数的规律，然后再带入欧拉函数就可以咯。

但是我没有找到规律，本来我是想乱猜的。但是放弃了。

官方证明：

由打表找规律可得，\sum_{k=0}^{x}(-1)^{k}2^{2x-2k}C_{2x-k+1}^{k}=x+1∑​k=0​x​​(−1)​k​​2​2x−2k​​C​2x−k+1​k​​=x+1，所以显然f_n(x)=n+x+1f​n​​(x)=n+x+1，因此直接求$\phi (n+x+1)$。时间效率O(T\sqrt{n})O(T√​n​​​)

严格证明：

设a_n=\sum_{k=0}^{n}(-1)^k2^{2n-2k}C_{2n-k+1}^ka​n​​=∑​k=0​n​​(−1)​k​​2​2n−2k​​C​2n−k+1​k​​

a_n=2^{2n}+\sum_{k=1}^{n}(-1)^k2^{2n-2k}(C_{2n-k}^k+C_{2n-k}^{k-1})=\sum_{k=0}^{n}(-1)^k2^{2n-2k}C_{2n-k}^k+\sum_{k=0}^{n-1}(-1)^{k+1}2^{2(n-1)-2k}C_{2(n-1)-k+1}^ka​n​​=2​2n​​+∑​k=1​n​​(−1)​k​​2​2n−2k​​(C​2n−k​k​​+C​2n−k​k−1​​)=∑​k=0​n​​(−1)​k​​2​2n−2k​​C​2n−k​k​​+∑​k=0​n−1​​(−1)​k+1​​2​2(n−1)−2k​​C​2(n−1)−k+1​k​​

设b_n=\sum_{k=0}^{n}(-1)^k2^{2n-2k}C_{2n-k}^kb​n​​=∑​k=0​n​​(−1)​k​​2​2n−2k​​C​2n−k​k​​，则b_n=a_n+a_{n-1}b​n​​=a​n​​+a​n−1​​

b_n=2^{2n}+\sum_{k=1}^{n-1}(-1)^k2^{2n-2k}(C_{2n-k-1}^k+C_{2n-k-1}^{k-1})+(-1)^nb​n​​=2​2n​​+∑​k=1​n−1​​(−1)​k​​2​2n−2k​​(C​2n−k−1​k​​+C​2n−k−1​k−1​​)+(−1)​n​​

=4\sum_{k=0}^{n-1}(-1)^k2^{2(n-1)-2k}C_{2(n-1)-k+1}^k+\sum_{k=0}^{n-1}(-1)^{k+1}2^{2(n-1)-2k}C_{2(n-1)-k}^k=4∑​k=0​n−1​​(−1)​k​​2​2(n−1)−2k​​C​2(n−1)−k+1​k​​+∑​k=0​n−1​​(−1)​k+1​​2​2(n−1)−2k​​C​2(n−1)−k​k​​

=4a_n-b_{n-1}=4a​n​​−b​n−1​​

得a_n-a_{n-1}=a_{n-1}-a_{n-2}a​n​​−a​n−1​​=a​n−1​​−a​n−2​​。因为a_0=1,a_1=1a​0​​=1,a​1​​=1，所以a_n=n+1a​n​​=n+1

证明比较费时，打表找规律能很快的得出解，所以本题的关键在于打表找规律。

#include <iostream>
#include <cstdio>
#include <cmath>

using namespace std;
long long f(long long n)
{
    long long cnt=n;
    for(long long i=2; i*i<=n; i++)
    {
        if(n%i==0)
        {
            cnt-=cnt/i;
            while(n%i==0)
                n/=i;
        }
    }
    if(n>1)
        cnt-=cnt/n;
    return cnt;
}
int main()
{
    long long n,x;
    while(~scanf("%lld",&n))
    {
        printf("%lld\n",f(n));
    }
}