BZOJ 3884 上帝与集合的正确用法 (欧拉定理)

最新推荐文章于 2019-09-18 10:41:05 发布

AbEver

最新推荐文章于 2019-09-18 10:41:05 发布

阅读量798

点赞数 3

分类专栏：数论 & 数学 BZOJ

本文链接：https://blog.csdn.net/Ab_Ever/article/details/70208793

版权

BZOJ 同时被 2 个专栏收录

70 篇文章 0 订阅

订阅专栏

数论 & 数学

28 篇文章 1 订阅

订阅专栏

Description

根据一些书上的记载，上帝的一次失败的创世经历是这样的：
第一天，上帝创造了一个世界的基本元素，称做“元”。
第二天，上帝创造了一个新的元素，称作“α”。“α”被定义为“元”构成的集合。容易发现，一共有两种不同的“α”。
第三天，上帝又创造了一个新的元素，称作“β”。“β”被定义为“α”构成的集合。容易发现，一共有四种不同的“β”。
第四天，上帝创造了新的元素“γ”，“γ”被定义为“β”的集合。显然，一共会有16种不同的“γ”。
如果按照这样下去，上帝创造的第四种元素将会有65536种，第五种元素将会有2^65536种。这将会是一个天文数字。
然而，上帝并没有预料到元素种类数的增长是如此的迅速。他想要让世界的元素丰富起来，因此，日复一日，年复一年，他重复地创造着新的元素……
然而不久，当上帝创造出最后一种元素“θ”时，他发现这世界的元素实在是太多了，以致于世界的容量不足，无法承受。因此在这一天，上帝毁灭了世界。
至今，上帝仍记得那次失败的创世经历，现在他想问问你，他最后一次创造的元素“θ”一共有多少种？
上帝觉得这个数字可能过于巨大而无法表示出来，因此你只需要回答这个数对p取模后的值即可。
你可以认为上帝从“α”到“θ”一共创造了10^9次元素，或10^18次，或者干脆∞次。
一句话题意：
这里写图片描述

Input

接下来T行，每行一个正整数p，代表你需要取模的值
对于100%的数据，T<=1000,p<=10^7

Output

T行，每行一个正整数，为答案对p取模后的值

Sample Input

3
2
3
6

Sample Output

0
1
4

Solution

这题真心是一道好题，题目很有趣，而且解题方法十分巧妙。不得不说出题人真心强无敌。
求2的无穷次幂Mod p的值，关于这种与无穷次计算相关的题目，我们可以类比当幂数很大时，我们该怎么做。

当计算步骤有限但很多时，我们可以想到用快速幂之类的方法快速计算。或者有些题目可以根据给出的数据的性质找出规律以推出公式(打表大法好。。。)

但在处理无穷多次运算时，我们必须找到门路，将其转化成有限次数的运算，方能用程序求解。有时候可以找到一个值q，发现计算q次后值不变了，或者说当次数趋于q时，值趋于一个定值f(q)。这时我们可以转化为有穷的情况来求解。但有时候无穷次计算和有限次计算迥然不同，举个栗子，就像无穷级数求和不能用求有限和的方法一样。这时就需要我们进行分析，将问题转化，最终的目标一定是可以用有限次计算得出的。

———————————————分割线———————————————————-

好像走远了。。
不要紧，回到此题。诚然，我们直接做是不行的。我们需要欧拉降幂(就是请欧拉出来帮我们降一下幂(๑•ᴗ•๑))，先讲以下有关欧拉定理的知识。

欧拉函数

$\phi(n) = n \times\prod (pi-1)/pi$

其中 $n=p1^{e1}p2^{e2}...pk^{ek}$
根据质数唯一分解定理，我们可以O( $\sqrt{n}$ )求一个数的欧拉函数。相信每个人都知道 $\phi(n)$ 还代表与不超过n且与n互质的数的个数，我们可以线性筛筛出欧拉函数表。

然而这不是重点。

欧拉定理

对于与x互质的数m,有 $x^{\phi(m)}\equiv1(mod\ m)$
欧拉定理可看作
费马小定理的推广，这里就不说费马小定理了。

回归题目，我们令 $p=2^k*q$ , 其中 $q$ 为奇数。明显这分解唯一。
然后就得到
$2^{2^{2^{2^{...}}}} mod\ p= 2^k(2^{2^{2^{2^{...}}}-k} \ mod \ q)$
我们设 $x=2^{2^{2^{...}}}-k$ ,
同时令 $w=x/\phi(q)$ , 得到 $2^x\ mod\ q=2^{x\ mod\ \phi(q)}*(2^{\phi(q)})^w\ mod\ q$
这数学公式可真心难打啊。。
接着根据欧拉定理(q与2互质)，有 $2^{\phi(q)}\equiv1(mod\ q)$ ,接着就得到了我们想要的
$2^k(2^{2^{2^{2^{...}}}-k} \ mod \ q)=2^k(2^{(2^{2^{2^{...}}}-k)\ mod \ \phi(q)} \ mod \ q)$
简直完美，因为我们发现
由于幂的无穷性，指数与原来一毛一样！好了我们递归吧。
引用出题人的话，我们可以发现

除第一次外模数都是偶数，故每次递归模数都会至少除掉2。因此在不超过O( $\log2(p)$ )次递归之后，模数就会变成1。由于任何数mod 1的结果都是0，故此时递归结束，回溯并计算答案。

欧拉函数处理 $\phi(1),\phi(2)$ ，值都是偶数！！感性一点，如果n(n>2)与x(x<=n)互质，那n-x也与n互质嘛。
我用线性筛来求欧拉函数，时间复杂度O(p+Tlog2p)，如果不考虑快速幂的话(考虑也是log级别的)。
到这一步，我们已经成功从嘴巴处AC此题，具体代码很简单。

代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#define N 10000000
using namespace std;


int T, p, q;
int cnt, prime[N+10], phi[N+10];
bool vis[N+10];

void Get_phi(){
    phi[1] = 1;
    vis[1] = true;
    for(int i = 2; i <= N; i++){
      if(!vis[i]){
        phi[i] = i - 1;
        prime[++cnt] = i;
      }
      for(int j = 1; j <= cnt && i * prime[j] <= N; j++){
        vis[i * prime[j]] = true;
        if(i % prime[j] == 0){
          phi[i * prime[j]] = phi[i] * prime[j];
          break;
        }
        else  phi[i * prime[j]] = phi[i] * (prime[j] - 1);
      }
    }
}

int Pow(long long x, int n, int Mod){
    if(n == 0)  return 1 % Mod;
    long long res = Pow(x*x%Mod, n>>1, Mod);
    if(n & 1)  res = res * x % Mod;
    return res;
}

int Work(int p){
    if(p == 1)  return 0;
    int k = 0;
    while(!(p & 1))  p >>= 1, k ++;
    int tmp = Work(phi[p]);
    tmp = (tmp%phi[p] - k%phi[p] + phi[p]) % phi[p];
    tmp = Pow(2, tmp, p);
    return tmp << k;
}

int main(){    

    scanf("%d", &T);

    Get_phi();

    while(T --){
      scanf("%d", &p);
      printf("%d\n", Work(p));
    }

    return 0;
}