版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/qq_34454069/article/details/83116679
题意:
分析:
很显然,以第一种方式的使用次数为自变量,以第二种方式的使用次数为因变量,一定是一个不上升函数。
换言之,已知f(x)为不下降函数,通过神奇的证明发现:f(x)+x为单谷函数。
所以。。。三分就行了第一种方式的使用次数就行了。。
求总方案可以用贪心的思想:从左往右依次处理每个位置,求出覆盖当前位置的最远点。如果要在当前位置加值,则每次都加到最远点为止。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define SF scanf
#define PF printf
#define MAXN 100010
#define INF 1000000000000000ll
#define MAXR 10000000ll
using namespace std;
typedef long long ll;
int n,h,t;
int p[MAXN],l[MAXN],r[MAXN];
void init(vector<int> val0,vector<int> a,vector<int> b,vector<int> m){
/*val0.resize(3);
a.resize(3);
b.resize(3);
m.resize(3);
for(int i=0;i<3;i++)SF("%d",&val0[i]);
for(int i=0;i<3;i++)SF("%d",&a[i]);
for(int i=0;i<3;i++)SF("%d",&b[i]);
for(int i=0;i<3;i++)SF("%d",&m[i]);*/
p[0]=val0[0];
for(int i=1;i<n;i++)
p[i]=(1ll*a[0]*p[i-1]+b[0])%m[0];
l[0]=val0[1];
r[0]=val0[2];
for(int i=1;i<h;i++){
l[i]=min(1ll*n-1,(1ll*a[1]*l[i-1]+b[1])%m[1]);
int dist=r[i-1]-l[i-1];
r[i]=min(1ll*n-1,l[i]+(1ll*a[2]*dist+b[2])%m[2]);
}
}
vector<int> v[MAXN];
int pre[MAXN];
ll calc(int x){
ll res=1ll*x*t;
ll sum=x;
int far=-1;
memset(pre,0,sizeof pre);
for(int i=0;i<n;i++){
sum+=pre[i];
for(int j=0;j<int(v[i].size());j++)
far=max(far,v[i][j]);
ll add=max(0ll,p[i]-sum);
if(add!=0&&far<i)
return INF;
sum+=add;
res+=add;
pre[far+1]-=add;
}
return res;
}
class SettingShield{
public:
ll getOptimalCost(int n1,int h1,int t1,vector<int> val0,vector<int> a,vector<int> b,vector<int> m){
n=n1;
h=h1;
t=t1;
//SF("%d%d%d",&n,&h,&t);
init(val0,a,b,m);
for(int i=0;i<h;i++)
v[l[i]].push_back(r[i]);
/*for(int i=0;i<n;i++)
PF("%d ",p[i]);
for(int i=0;i<h;i++)
PF("\n%d %d",l[i],r[i]);*/
ll lx=0,rx=MAXR,ans=MAXR;
while(lx<=rx){
ll midl=lx+(rx-lx)/3;
ll midr=rx-(rx-lx)/3;
ll v1=calc(midl),v2=calc(midr);
//PF("<%lld %lld %lld:%lld %lld:%lld>\n",lx,rx,midl,v1,midr,v2);
if(v1==INF||v1>v2){
lx=midl+1;
ans=midr;
}
else{
rx=midr-1;
ans=midl;
}
}
return calc(ans);
}
};