二分与倍增
倍增的定义请参考李煜东大佬的书,我们这里给出二分和倍增的模板。
二分:
int l=0,r=v.size()-1;
while(l<r)
{
int mid=l+r>>1;
if(v[mid]>=x) r=mid;
else l=mid+1;
}
倍增:
int r=0,p=1;
while(p)
{
if(r+p<v.size() && v[r+p]<=x) r+=p,p<<=1;
else p>>=1;
}
我可以观察到其实倍增是二分的逆运算,二分与x在数组的位置无关,但倍增却是与x有关的,二分的时间复杂度稳定在( l o g n logn logn),倍增在最优O(1),最差( 2 l o g n 2logn 2logn),普通情况下用二分就够了,一般来说二分比较快,但也有例外,比如测试集中存在大量的比较近的x,我们可以用倍增。
天才ACM
给定一个整数 M,对于任意一个整数集合 S,定义“校验值”如下:
从集合 S 中取出 M 对数(即 2∗M 个数,不能重复使用集合中的数,如果 S 中的整数不够 M 对,则取到不能取为止),使得“每对数的差的平方”之和最大,这个最大值就称为集合 S 的“校验值”。
现在给定一个长度为 N 的数列 A 以及一个整数 T。
我们要把 A 分成若干段,使得每一段的“校验值”都不超过 T。
求最少需要分成几段。
输入格式
第一行输入整数 K,代表有 K 组测试数据。
对于每组测试数据,第一行包含三个整数 N,M,T 。
第二行包含 N 个整数,表示数列A1,A2…AN。
输出格式
对于每组测试数据,输出其答案,每个答案占一行。
数据范围
1≤K≤12,
1≤N,M≤500000,
0≤T≤1018,
0≤Ai≤220
输入样例:
2
5 1 49
8 2 1 7 9
5 1 64
8 2 1 7 9
输出样例:
2
1
这题的思路就是不断归并集合,并用倍增优化。我们可以看到即使我们最坏时间复杂度高达(nlogn)依然能过,是因为倍增优化了许多离原点很近的方案。
#include<iostream>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=500010;
int a[N],t1[N],t2[N];
int n,m;
ll s;
bool check(int l,int mid,int r)
{
for(int i=mid;i<=r;i++) t1[i]=a[i];
sort(t1+mid,t1+r+1);
int i=l,j=mid;
for(int k=l;k<=r;k++)
{
if(j>r||(i<mid&&t1[i]<=t1[j])) t2[k]=t1[i++];
else t2[k]=t1[j++];
}
ll res=0;
int left = l,right=r,cnt=0;
while(cnt<m&&right>left)
{
res+=1ll * (t2[right]-t2[left]) * (t2[right] - t2[left]);
right--,left++,cnt++;
}
if(res<=s)
{
for(int k=l;k<=r;k++) t1[k]=t2[k];
return true;
}
return false;
}
int main(){
//freopen("data.in","r",stdin);
//freopen("data.out","w",stdout);
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d%d%lld",&n,&m,&s);
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
int l=1,r=1,res=0;
while(l<=n)
{
int p=1;
t1[l]=a[l];
res++;
while(p)
{
if(r+p<=n&&check(l,r+1,r+p)) r+=p,p<<=1;
else p>>=1;
}
l=++r;
}
printf("%d\n",res);
}
return 0;
}
ST算法
给大家推荐一下poj 3264。给大家写一下代码,思路请参考李煜东大佬。
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdio>
using namespace std;
const int N=100010;
int a[N];
int fs[N][50],fb[N][50];
int n;
void work()
{
for(int i=1;i<=n;i++){
fs[i][0]=a[i];
fb[i][0]=a[i];
}
int t=log(n)/log(2)+1;
for(int j=1;j<t;j++)
for(int i=1;i<=n;i++)
{
fb[i][j]=max(fb[i][j-1],fb[i+(1<<j-1)][j-1]);
fs[i][j]=min(fs[i][j-1],fs[i+(1<<j-1)][j-1]);
}
}
int query(int l,int r)
{
int k=log(r-l+1)/log(2);
int MAX=max(fb[l][k],fb[r-(1<<k)+1][k]);
int MIN=min(fs[l][k],fs[r-(1<<k)+1][k]);
return MAX-MIN;
}
int main(){
//freopen("data.in","r",stdin);
//freopen("data.out","w",stdout);
int T;
scanf("%d%d",&n,&T);
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
work();
while(T--)
{
int l,r;
scanf("%d%d",&l,&r);
printf("%d\n",query(l,r));
}
return 0;
}