不妨把求φ抽象成把将每个位置上的一个小球左移一格并分裂的过程,那么即求所有球都被移到1号格子的步数. 显然要达到1必须先到达2.可以发现每次分裂一定会分裂出2号位的球,因为2以外的质数一定是奇数.以及,每次移动至多将一个2号位的球移至1号位. 于是我们只要数出每个位置能将几个球分裂至2号位就可以了.注意初始时若2号位没有球答案要+1. 这个数数可以用线性筛搞定.n为质数则有f[n]=f[n-1],否则有f[n]=f[prime]+f[n/prime]. (挤进bzoj前1k了www #inclu…

根绝欧拉函数的计算式,欧拉函数次方每次就是不断把2消掉,并把其他质因子不断变出2,最终弄成1显然我们先预处理每个数会被弄出多少个2出来,次方=弄出来的2的总数-[一开始是否有2](因为一开始没2的话是要先次方一下,无法消掉1个2) ..] of longint; f:..] of boolean; t,i,j,n,x,y:longint; ans:int64; ch:boolean; begin w[]:=; do begin then begin inc(t); p[t]:=i; w[i]:=…

2749: [HAOI2012]外星人 Time Limit: 3 Sec  Memory Limit: 128 MBSubmit: 377  Solved: 199[Submit][Status] Description   Input   Output 输出test行,每行一个整数,表示答案. Sample Input 122 23 1 Sample Output 3 HINT Test<=50 Pi<=10^5,1<=Q1<=10^9 Source 题解: 终于把这题搞掉了.…

外星人 Time Limit: 3 Sec  Memory Limit: 128 MB[Submit][Status][Discuss] Description Input   Output 输出test行,每行一个整数,表示答案. Sample Input 1 2 2 2 3 1 Sample Output 3 HINT Test<=50 Pi<=10^5,1<=Q1<=10^9 Main idea 给定一个数,用Πp[i]^q[i](p<=10^5,q<=10^9)…

2749: [HAOI2012]外星人 Time Limit: 3 Sec  Memory Limit: 128 MBSubmit: 677  Solved: 360[Submit][Status][Discuss] Description Input Output 输出test行,每行一个整数,表示答案. Sample Input 1 2 2 2 3 1 Sample Output 3       [题解]   注意到只有phi[2]=1,所以每个数最终都是要经过phi[2]=1这一步变为1…

1.快速幂 计算a^b的快速算法,例如,3^5,我们把5写成二进制101,3^5=3^1*1+3^2*2+3^4*1 ll fast(ll a,ll b){ll ans=;,a=mul(a,a)))ans=mul(ans,a);return ans;}//一行快速幂 2.快速乘 当模数较大时,直接乘会爆掉long long,需要快速乘法. 即用浮点计算倍数,做差相当于计算余数模2^63的结果,然后再模一下就好了(因为余数不超过long long) typedef long long ll; ll…

NOIP以前可能会持续更新 写在前面 NOIP好像马上就要到了,感觉在校内训练里面经常被虐有一种要滚粗的感觉(雾.不管是普及组还是提高组,我都参加了好几年了,结果一个省一都没有,今年如果还没有的话感觉就真的要滚大粗退役回去念书了QAQ.于是有了压力就来刷(水水水)题.感觉校内OJ的题库还挺多的就开始做校内OJ的题.(本校的其他神犇都在其他各种OJ上屠丧题我感觉好虚啊!)于是把这几年NOIP的原题拿出来做了下. (我蛮立个flag:如果NOIP过了就买BZOJ权限号...) 历年NOIP提高组一句…