二分答案后,贪心做。
每个数对应一个区间,每次取能取到的最小的数,如果取不到了,就无解。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iostream>
#define maxn 50001000
using namespace std;
int a[maxn];
int n,mod,ans;
long long sa,sb,sc,sd;
long long F(long long x)
{
long long x2=x*x%mod,x3=x2*x%mod;
return (sa*x3%mod+sb*x2%mod+sc*x%mod+sd)%mod;
}
bool check(int x)
{
int mx=1;
for (int i=1;i<=n;i++)
{
mx=max(mx,a[i]-x);
if (mx>a[i]+x) return 0;
}
return 1;
}
int main()
{
scanf("%d%d%d%d%d%lld%d",&n,&sa,&sb,&sc,&sd,&a[1],&mod);
a[0]=0;
for (int i=2;i<=n;i++) a[i]=(F(a[i-2])+F(a[i-1]))%mod;
int l=0,r=mod;
while (l<=r)
{
int mid=(l+r)/2;
if (check(mid)) ans=mid,r=mid-1; else l=mid+1;
}
printf("%d\n",ans);
return 0;
}