Upd2019.4.19

yy了一个新的写法
新的东西

CODE:

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
using namespace std;
typedef long long LL;
const int N=3000005;
int MOD;
int f[N];
int add (int x,int y)   {x=x+y;return x>=MOD?x-MOD:x;}
int mul (int x,int y)   {return (LL)x*y%MOD;}
int Pow (int x,int y)
{if (y==0) return 1;if (y==1) return x;int lalal=Pow(x,y>>1);lalal=mul(lalal,lalal);if (y&1) lalal=mul(lalal,x);return lalal;
}
int g[N];
void Init (int n)
{f[0]=0;g[1]=f[1]=1;for (int u=2;u<=n;u++) g[u]=f[u]=add(add(f[u-1],f[u-1]),f[u-2]);for (int u=1;u<=n;u++){  int Inv=Pow(g[u],MOD-2);for (int i=u+u;i<=n;i+=u)  g[i]=mul(g[i],Inv);}
}
int main()
{int T;scanf("%d",&T);while (T--){int n;scanf("%d%d",&n,&MOD);Init(n);int lalal=1,ans=0;for (int u=1;u<=n;u++) {lalal=mul(lalal,g[u]);ans=add(ans,mul(lalal,u));}printf("%d\n",ans);}return 0;
}

bzoj 4833: [Lydsy1704月赛]最小公倍佩尔数相关推荐

BZOJ 4833: [Lydsy1704月赛]最小公倍佩尔数（数论 + 最值反演）
题面令 ${(1+\sqrt 2)}^n=e(n)+f(n)*\sqrt2$ ,其中 $e(n),f(n)$ 都是整数,显然有 \({(1-\sqrt 2)}^n=e(n)-f(n)*\sq ...
【bzoj 4833】[Lydsy1704月赛]最小公倍佩尔数
Description 令 $(1+\sqrt 2)^n=e(n)+\sqrt 2\cdot f(n)$ ,其中 $e(n),f(n)$ 都是整数,显然有 $(1-\sqrt 2)^n=e(n)-\s ...
BZOJ4833: [Lydsy1704月赛]最小公倍佩尔数-数论
传送门题意: 令(1+2–√)n=e(n)+f(n)∗2–√(1+2)n=e(n)+f(n)∗2(1+\sqrt2)^n=e(n)+f(n)*\sqrt2,其中e(n),f(n)e(n),f(n)e ...
BZOJ4833 [Lydsy1704月赛]最小公倍佩尔数
题意已知$e_n+\sqrt2f_n=(1+\sqrt2)^n$,$e_n-\sqrt2f_n=(1-\sqrt2)^n$,$g_n=lcm_{i=1}^nf_i$,求\(\sum_{i ...
BZOJ4833: [Lydsy1704月赛]最小公倍佩尔数
Problem 传送门 Sol 容易得到 fn=en−1+fn−1,en−1=fn−1+en−1,f1=e1=1f_n=e_{n-1}+f_{n-1},e_{n-1}=f_{n-1}+e_{n-1}, ...
【BZOJ4833】最小公倍佩尔数（min-max容斥）
[BZOJ4833]最小公倍佩尔数(min-max容斥) 题面 BZOJ 题解首先考虑怎么求$f(n)$,考虑递推这个东西 \((1+\sqrt 2)(e(n-1)+f(n-1)\sqrt 2) ...
[bzoj 4833]最小公倍佩尔数
传送门 Description Let $(1+\sqrt2)^n=e(n)+f(n)\cdot\sqrt2$ , both $e(n)$ and $f(n)$ are integers ...
[数论反演] BZOJ 4833 最小公倍佩尔数
当时比赛时灵机一动把gigi−1g_i\over g_{i-1}喂给了OEIS 然后就找到了 233 就是这个咯然后就水过去了题解?题解我还没看先挖个坑 UPD:跟这个题是一毛一样的咯 #in ...
[bzoj4833][数论][min-max容斥]最小公倍佩尔数
Description 令(1+sqrt(2))n=e(n)+f(n)*sqrt(2),其中e(n),f(n)都是整数,显然有(1-sqrt(2))n=e(n)-f(n)*sqrt(2).令g( n) ...

bzoj 4833: [Lydsy1704月赛]最小公倍佩尔数

Upd2019.4.19

bzoj 4833: [Lydsy1704月赛]最小公倍佩尔数相关推荐

最新文章

热门文章