题目概述
求 222⋮ mod p (无数个 2 )。
解题报告
直接上降幂大法,直到模数变为
可以证明
示例程序
#include<cstdio>
#include<cctype>
#include<algorithm>
#define fr first
#define sc second
#define mp make_pair
using namespace std;
typedef long long LL;typedef pair<int,bool> data;
const int maxm=10000000;
int te,m;int p[maxm+5],phi[maxm+5];bool pri[maxm+5];
#define Eoln(x) ((x)==10||(x)==13||(x)==EOF)
inline char readc()
{
static char buf[100000],*l=buf,*r=buf;
if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
if (l==r) return EOF;return *l++;
}
inline int readi(int &x)
{
int tot=0,f=1;char ch=readc(),lst='+';
while (!isdigit(ch)) {if (ch==EOF) return EOF;lst=ch;ch=readc();}
if (lst=='-') f=-f;
while (isdigit(ch)) tot=(tot<<3)+(tot<<1)+ch-48,ch=readc();
return x=tot*f,Eoln(ch);
}
void Make(int m)
{
pri[1]=true;phi[1]=1;
for (int i=2;i<=m;i++)
{
if (!pri[i]) p[++p[0]]=i,phi[i]=i-1;
for (int j=1,t;j<=p[0]&&(t=p[j]*i)<=m;j++)
{
if (i%p[j]) phi[t]=phi[i]*phi[p[j]],pri[t]=true; else
{phi[t]=phi[i]*p[j];pri[t]=true;break;}
}
}
}
inline data Pow(LL w,int b,int MOD)
{
LL s=1;bool fl=true;
while (b)
{
if (b&1) {if ((s*=w)>=MOD) fl=false;s%=MOD;}
b>>=1;if (b) {if ((w*=w)>=MOD) fl=false;w%=MOD;}
}
return mp(s,fl);
}
data Solve(int st,int MOD)
{
if (MOD==1) return mp(0,false);
data b=Solve(st+1,phi[MOD]);if (MOD&1) return Pow(2,b.fr,MOD);
if (!b.sc) b.fr+=phi[MOD];return Pow(2,b.fr,MOD);
}
int main()
{
freopen("program.in","r",stdin);
freopen("program.out","w",stdout);
readi(te);Make(maxm);
while (te--) readi(m),printf("%d\n",Solve(1,m).fr);
return 0;
}