其实就是求三维偏序最长链。
我们考虑cdq分治+dp。
注意要求严格小于,因此cdq分治时要注意一些细节的调整。
别忘了初始化!【哭泣】
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <set>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 50010
inline char gc(){
static char buf[1<<16],*S,*T;
if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
return *S++;
}
inline int read(){
int x=0,f=1;char ch=gc();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
return x*f;
}
int aa[N],n,m,ans=0,c[N];
struct data{
int x,y,z,f;
friend bool operator<(data a,data b){return a.y<b.y;}
}e[N];
inline bool cmp(data a,data b){return a.x==b.x?(a.y==b.y?a.z<b.z:a.y<b.y):a.x<b.x;}
inline void add(int x,int val){for(;x<=m;x+=x&-x) c[x]=max(c[x],val);}
inline void clr(int x){for(;x<=m;x+=x&-x) c[x]=0;}
inline int ask(int x){int res=0;for(;x;x-=x&-x) res=max(res,c[x]);return res;}
inline void cdq(int l,int r){
if(l==r) return;int mid=l+r>>1;
while(mid>=l&&e[mid].x==e[mid-1].x) --mid;//严格小于,调整一下
if(mid<l) return;cdq(l,mid);sort(e+l,e+mid+1);
sort(e+mid+1,e+r+1);int p1=l,p2=mid+1;
while(p1<=mid&&p2<=r){
if(e[p1]<e[p2]) add(e[p1].z,e[p1].f),++p1;
else e[p2].f=max(e[p2].f,ask(e[p2].z-1)+1),++p2;
}while(p2<=r) e[p2].f=max(e[p2].f,ask(e[p2].z-1)+1),++p2;
for(int i=l;i<p1;++i) clr(e[i].z);sort(e+mid+1,e+r+1,cmp);cdq(mid+1,r);
}
int main(){
// freopen("a.in","r",stdin);
int x=read(),mod=read(),y=1;n=read();e[0].x=2000000000;
for(int i=1;i<=n;++i){
y=(ll)y*x%mod;e[i].x=y;e[i].f=1;
y=(ll)y*x%mod;e[i].y=y;
y=(ll)y*x%mod;e[i].z=y;
if(e[i].x>e[i].y) swap(e[i].x,e[i].y);
if(e[i].x>e[i].z) swap(e[i].x,e[i].z);
if(e[i].y>e[i].z) swap(e[i].y,e[i].z);aa[i]=e[i].z;
}sort(aa+1,aa+n+1);m=unique(aa+1,aa+n+1)-aa-1;
for(int i=1;i<=n;++i) e[i].z=lower_bound(aa+1,aa+m+1,e[i].z)-aa;
sort(e+1,e+n+1,cmp);cdq(1,n);
for(int i=1;i<=n;++i) ans=max(ans,e[i].f);
printf("%d\n",ans);
return 0;
}