【bzoj3613】[Heoi2014]南园满地堆轻絮 二分答案

二分答案后,贪心做。

每个数对应一个区间,每次取能取到的最小的数,如果取不到了,就无解。


#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;
}


你可能感兴趣的:(【bzoj3613】[Heoi2014]南园满地堆轻絮 二分答案)