码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ3613: [Heoi2014]南园满地堆轻絮

时间:2018-05-07 22:57:31      阅读:204      评论:0      收藏:0      [点我收藏+]

标签:ret   题意   algo   构造   string   时间   mes   int   name   

分析:

构造数据时间有些长,可以用秦九韶优化一下。

二分答案+贪心,即:另每一个b[i]尽可能的小的同时满足题意,在枚举过程中,判断是否存在一个b[i-1]>a[i]+x

如果存在,那么向右找

如果不存在,向左找

附上代码:

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <queue>
//#include <ctime>
using namespace std;
#define N 5000005
#define ll long long
int a[N],f[N],b[N],n,Sa,Sb,Sc,Sd,mod;
int F(int x)
{
	if(f[x]) return f[x];
	return f[x]=((1ll*Sa*a[x]%mod+Sb)*a[x]*a[x]%mod+(1ll*Sc*a[x]%mod+Sd)%mod)%mod;
}
int check(int x)
{
	b[0]=-1<<30;
	for(int i=1;i<=n;i++)
	{
		if(a[i]>b[i-1]+x)
		{
			b[i]=a[i]-x;
		}else
		{
			if(a[i]+x<b[i-1])return 0;
			b[i]=b[i-1];
		}
	}
	return 1;
}
int main()
{
	a[0]=0;
	scanf("%d%d%d%d%d%d%d",&n,&Sa,&Sb,&Sc,&Sd,&a[1],&mod);
	for(int i=2;i<=n;i++)a[i]=(F(i-1)+F(i-2))%mod;
	int l=0,r=1<<30;
	while(l<r)
	{
		int m=(l+r)>>1;
		if(check(m)) r=m;
		else l=m+1;
	}
	printf("%d\n",l);
	return 0;
}

  卡了一下常数(loj上需要卡一下常数,bzoj和洛谷不用,开O2跑的飞快)

BZOJ3613: [Heoi2014]南园满地堆轻絮

标签:ret   题意   algo   构造   string   时间   mes   int   name   

原文地址:https://www.cnblogs.com/Winniechen/p/9005185.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!