前序知识点
:洛谷P1115 最大子段和
三种办法求最大子段和
第一种:二分
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 201000
#define MIN 0xc0c0c0c0
using namespace std;
int n,a[N],ans;
int work(int l,int r)
{
if( l==r ) return a[l];
int mid = (l+r)>>1,sum = 0 ,max1=MIN,max2=MIN;
for(int i=mid;i>=l;i--)
{
sum +=a[i];
max1 = max( max1 , sum );
}
sum = 0;
for(int i=mid+1;i<=r;i++)
{
sum += a[i];
max2 = max( max2 , sum );
}
return max( max( work(l,mid) , work( mid+1,r ) ) , max1+max2 );
}
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
}
ans = work(1,n);
cout<<ans;
return 0;
}
第二种:DP
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 201000
#define MIN 0xc0c0c0c0
using namespace std;
int n,a[N],dp[N],ans = MIN;
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
dp[i] = max( dp[i-1]+a[i] , a[i] );
ans = max( ans , dp[i] );
}
cout<<ans;
return 0;
}
第三种:单调队列
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 201000
#define INF 0x3f3f3f3f3f
#define MIN 0xc0c0c0c0
using namespace std;
int n,a[N],minn=INF,sum[N],ans=MIN;
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
sum[i] = sum[i-1] + a[i];
minn = min( minn , sum[i-1] );
ans = max( ans , sum[i] - minn );
}
cout<<ans;
return 0;
}
正题
idea:首先先乘后加必然不如先加再乘法
好了回归正题:既然要使用二分,那么该问题的单调性是什么呢?
魔法1 :如果使用 n 次可以满足要求,那么 n+1 次的魔法也一定满足要求 。
魔法2 :如果使用 n 次可以满足要求,那么 n+1 次的魔法也一定满足要求 。
但是注意如果枚举加法,TLE送走你,枚举乘法的话是logS级别的,可以接受。
二分的一个小模板是这样的:
void erfen(int now)
{
l = 0,r = MAX;
while(l<r)
{
mid = (l+r)/2;
if( check(mid)==true ) r = mid;
else l = mid+1;
}
}
//check函数是检查魔法1使用次数为now,枚举魔法2使用次数的情况下最大子段和是否大于等于S
ACcode:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MIN 0xc0c0c0c0
#define INF 0x3f3f3f3f
#define N 101000
#define LL long long
using namespace std;
LL lo,n,dp[N],a[N],s,A,B,ans = 1e10;
bool check(LL mid,LL cnt)
{
LL maxx = 0;
memset(dp,0,sizeof(dp));
for(int i=1;i<=n;i++)
{
dp[i] = max( dp[i-1] + a[i] , a[i] ) + mid;
maxx = max( maxx , dp[i] );
}
if( !maxx ) return false;
for(LL i=1;i<=cnt&&maxx<s;i++) maxx*=2;
return maxx>=s;
}
LL work(int cnt)
{
LL l = 0,r = 3e9;
while( l<r )
{
LL mid = (l+r)/2;
if( check(mid,cnt)==true ) r = mid;
else l = mid+1;
}
return l*A + cnt*B;
}
int main()
{
cin>>n>>A>>B>>s;
for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
lo = log2(s)+1;
for(int i=0;i<=lo;i++)
{
ans = min( ans , work(i) );
}
cout<<ans;
return 0;
}