[ NOI 2002 ] Robot

\(\\\)

Description


\(\\\)

Solution


垃圾语文题毁我青春

这题其实就是重定义了俩函数....

首先 \(\varphi(1)=0\)

然后 \(2\) 在计算 \(\mu\) 的时候不算做质因子,含 \(2\) 的数 \(\mu\)\(=0\)

\(m\) 的所有因数中 , \(\mu=1,-1,0\)\(\varphi\) 值之和。

\(\\\)

先考虑 \(\mu=1,-1\) 的答案。

根据 \(\mu\) 的定义,此时所有所求的数字都满足,分解质因数每一个质因数的质数至多为 \(1\)

然后 **仅对于这些数 ** 欧拉函数就有了一个新的计算公式:

\[\varphi(n)=n\times {\frac{p_1-1}{p_1}}\times {\frac{p_2-1}{p_2}}\times ...\times {\frac{p_m-1}{p_m}}=(p_1-1)\times (p_2-1)\times ...\times (p_3-1) \]

发现这是一个质数集合选哪些的问题。

注意到每次多选一个质数,所有的数字答案都乘上了相同的数,所以可以直接求答案。

用一种类似动规的写法,依次考虑每一个位置选或不选。

\(ans_1\)表示选奇数个质因子,\(ans_2\) 表示选偶数个质因子的答案。

\[ans_1+=ans_2\times (p_i-1),ans_2+=ans_1\times (p_i-1) \]

可以思考一下加号连接的含义。

\(\\\)

然后就只需要解决 \(\mu=0\) 的部分了。

关于欧拉函数其实有一个等式

\[\sum_{d|m}\varphi(d)=m \]

然后问题解决,答案就是 \(m-ans_1-ans_2-1\),之所以多减掉一个 \(1\) 是因为 \(\varphi(1)=0\)

\(\\\)

Code


#include<cmath>
#include<cstdio>
#include<cctype>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define mod 10000
#define R register
#define gc getchar
using namespace std;
 
inline int rd(){
  int x=0; bool f=0; char c=gc();
  while(!isdigit(c)){if(c=='-')f=1;c=gc();}
  while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=gc();}
  return f?-x:x;
}
 
inline int qpow(int x,int t){
    int res=1;
    while(t){
      if(t&1) (res*=x)%=mod;
      (x*=x)%=mod; t>>=1;
    }
    return res;
}
 
int n,m=1,p,cnt,tmp,ansodd,anseven=1;
 
int main(){
  int n=rd();
  for(R int i=1;i<=n;++i){
    p=rd(); cnt=rd();
    (m*=qpow(p,cnt))%=mod;
    if(p>2){
      tmp=ansodd;
      (ansodd+=anseven*(p-1))%=mod;
      (anseven+=tmp*(p-1))%=mod;
    }
  }
  printf("%d\n%d\n%d\n",(anseven-1+mod)%mod,ansodd,((m-anseven-ansodd)%mod+mod)%mod);
  return 0;
}
posted @ 2018-11-06 18:33  SGCollin  阅读(105)  评论(0编辑  收藏  举报