算法介绍:
折半搜索常用于复杂度O(n!)级的搜索问题,当我们发现很显然可以将问题划分为两部分分别搜索枚举,再合二为一求出最终答案时,我们可以选择使用折半搜索。
常见数据规模:
对于答案的值域往往没有要求,只对给出元素个数 n n n 有一定要求:
n ≤ 50 n\leq50 n≤50
例题:
来源:东BOJ:oj.neu.edu.cn
解题思路:
设目标值x为 g o a l goal goal,最终种类数为 a n s ans ans。
如果纯暴力解决,算法复杂度为 O ( 2 n ) O(2^n) O(2n),而 n n n最大可达到 40 40 40,显然是会超时的
所以我选择先枚举出前半部分元素,即 t 1 ∼ t n / 2 {t_1} \sim {t_{n/2}} t1∼tn/2所能生成的所有值 t o t 1 tot_1 tot1所组成的集合 s s s。
再去搜索
t
n
/
2
+
1
∼
t
n
{t_{n/2+1}} \sim {t_{n}}
tn/2+1∼tn能组合产生的所有的值
t
o
t
2
tot_2
tot2,对于每次搜索产生的
t
o
t
2
tot_2
tot2,都去
s
s
s中二分搜索出
x
−
t
o
t
2
x-tot_2
x−tot2的值,如果
x
−
t
o
t
2
x-tot_2
x−tot2存在,就是找到可行方案了,It’s MYGO!, 就把
t
y
p
e
[
x
−
t
o
t
2
]
type[x-tot_2]
type[x−tot2]加到最终答案
a
n
s
ans
ans中。
如果没有找到,就不是可行方案,乐队就要解散了(大悲)。
关于我曲折的debug过程
最开始是
T
L
E
TLE
TLE,这很正常,毕竟
s
s
s规模还是挺大的,硬搜肯定会
T
T
T掉。
显然优化就是
u
n
i
q
u
e
unique
unique去重,然后统计
s
s
s每个元素的出现次数。
结果还是过不了,一堆
W
A
WA
WA。
心态就彻底炸了
这题都做不出来是不是该速速remake了
胃疼头疼的debuff就一起上了。
最后发现,是
u
n
i
q
u
e
unique
unique写错了…
本来我定义的
s
s
s大小是
c
n
t
cnt
cnt,实际我给写成
n
n
n了…
正确代码:
错误代码
无话可说。。。。总之最后是过了,还算不错
完整代码
#include<bits/stdc++.h>
using namespace std;
const int maxn=2e6+60;
long long si[maxn];
int cnt,n;
long long a[50],goal;
long long tot1=0,tot2=0;
long long minn;
long long type[maxn];
int len=0;
void dfs1(int x)
{
if(x==n/2+1)
{
si[++cnt]=tot1;
return;
}
tot1+=a[x];
if(tot1<=goal) dfs1(x+1);
tot1-=a[x];
dfs1(x+1);
}
long long ans=0;
int ans_id=0;
bool jud(long long res)
{
int l=1,r=len;
while(l<r)
{
int mid=(l+r+1)>>1;
if(si[mid]==res)
{
ans_id=mid;
return true;
}
if(si[mid]<=res) l=mid;
else r=mid-1;
}
return false;/*
for(int i=1;i<=cnt;++i)
{
if(si[i]==res) ++ans;
}
return false;*/
}
void dfs2(int x)
{
if(x==n+1)
{
if(tot2==goal) ans+=type[1];
else if(jud(goal-tot2)) ans+=type[ans_id];
return;
}
tot2+=a[x];
if(tot2+minn<=goal) dfs2(x+1);
tot2-=a[x];
dfs2(x+1);
}
int main()
{
scanf("%d%lld",&n,&goal);
for(int i=1;i<=n;++i) scanf("%lld",&a[i]);
dfs1(1);
minn=si[1];
for(int i=2;i<=cnt;++i) minn=min(minn,si[i]);
sort(si+1,si+cnt+1);
long long now=si[1];
long long cnt_now=1;
len=0;
for(int i=2;i<=cnt;++i)
{
if(si[i]==now) ++cnt_now;
else
{
type[++len]=cnt_now;
now=si[i];
cnt_now=1;
}
}
type[++len]=cnt_now;
unique(si+1,si+cnt+1)-si-1;
dfs2(n/2+1);
printf("%lld\n",ans);
// printf("len=%d\n",len);
// for(int i=1;i<=len;++i) printf("%d ",type[i]);
return 0;
}
/*
4 2
1 1 1 1
*/
题都看完了,就来推一首MYGO翻唱吧
【歌ってみた】少女レイ(少女REI) covered by 燈
一首不够?再来一个吧
【歌ってみた】「二息歩行 (Reloaded)」covered by 燈