[HEOI2014]南园满地堆轻絮

原题链接

思路一 (二分+贪心)：

二分答案 ans 每个数就变成一个区间 a[i]-ans,a[i]+ans ,问题转化为 判断n个区间各任选一点是否可以形成一个非严格递增的序列 ，解决这个问题贪心即可 。 具体看代码

参考代码：

#include <iostream>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <algorithm>
#include <vector>
#include <string>
#include <iomanip>
#include <cmath>
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <climits>
//#include <unordered_map>
#define guo312 std::ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)
#define ll long long
#define Inf LONG_LONG_MAX
#define inf INT_MAX
#define endl "\n"
#define PI 3.1415926535898
using namespace std;
const int N=5e6+10;
ll sa,sb,sc,sd,mod,n;
ll a[N],b[N];
ll f(ll x){     // 生成数据
    return (sa * x % mod * x % mod * x % mod + sb * x % mod * x % mod + sc * x % mod + sd) % mod;
}
bool check(ll x){   // 贪心选点 （每个数都是越小越好）
	b[1]=max((ll)1,a[1]-x);  // b[1] 越小越好
	for(int i=2;i<=n;i++){
		if(a[i]>=b[i-1]){   // 本身就大于 b[i-1]
			b[i]=max(b[i-1],a[i]-x);  
		}
		else{
			if(a[i]+x<b[i-1]) return 0;  // 不满足单调递增
			else b[i]=b[i-1];   // 保证 b[i] 尽量小
		}
	}
	return 1;
}
int main(){
guo312;
	cin>>n>>sa>>sb>>sc>>sd>>a[1]>>mod;
	for(int i=2;i<=n;i++){  // 生成数据
		a[i]=(f(a[i-1])+f(a[i-2]))%mod;
	}
	ll l=0,r=mod;  // 二分答案 b[i] 是非负数 
	while(l<r){
		ll mid=l+r>>1;
		if(check(mid)) r=mid;
		else l=mid+1;
	}
	cout<<l;
	return 0;
}

思路二（逆序对）：

仔细观察会发现 影响单调递增的是逆序对

这里需要知道两个数 x y （x>y) 要使 x<=y 对本题最好就是 x-(x-y+1)/2 y+(x-y+1)/2 （加1是存在差值是奇数的情况)

寻找差值最大的逆序对即可

参考代码：

#include <iostream>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <algorithm>
#include <vector>
#include <string>
#include <iomanip>
#include <cmath>
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <climits>
//#include <unordered_map>
#define guo312 std::ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)
#define ll long long
#define Inf LONG_LONG_MAX
#define inf INT_MAX
#define endl "\n"
#define PI 3.1415926535898
using namespace std;
const int N=5e6+10;
ll sa,sb,sc,sd,mod,n;
ll a[N],b[N];
ll f(ll x){
    return (sa * x % mod * x % mod * x % mod + sb * x % mod * x % mod + sc * x % mod + sd) % mod;
}
int main(){
guo312;
	cin>>n>>sa>>sb>>sc>>sd>>a[1]>>mod;
	ll re=0,maxn=a[1]; // maxn记录前面的最大值
	for(int i=2;i<=n;i++){   // 生成数
		a[i]=(f(a[i-1])+f(a[i-2]))%mod;
		re=max(re,maxn-a[i]);  // 更新最大差值
		maxn=max(maxn,a[i]);  // 更新最大值
	}
	cout<<(re+1)/2;
	return 0;
}