[luogu1526]智破连环阵

（以下在描述复杂度时，认为$n$和$m$相同，因此一律使用$n$）

称第$i$个炸弹能匹配非空区间$[l,r]$，当且仅当$l$到$r$内所有武器都在$i$攻击范围内，且$r=m$或第$r+1$个武器不在$i$攻击范围内

（更通俗的表示即在第$l$个武器启动时，使用第$i$个炸弹恰好会攻击到第$r$个武器）

关于这个区间的匹配，显然可以在$o(n^{2})$的时间内预处理

下面，将$n$个武器划分为若干个区间，分别表示每一个炸弹所破坏的武器（恰好），爆搜来确定区间划分后，只需判定其是否合法，并用区间数对答案取$min$即可

关于判定，将第$i$个炸弹向划分出的区间中能匹配的区间连边，判定最大匹配是否等于区间数即可

暴搜复杂度为$o(n^{2}2^{n})$（即每一个位置是否划分，以及二分图匹配），下面来考虑剪枝——

剪枝1：求出$[i,n]$这些武器，在允许重复利用炸弹的情况下的最少次数$g_{i}$

若当前划分出$now$个区间，最后一个区间右端点为$r$，最小答案为$ans$，那么若$now+g_{r+1}ge ans$显然不需要搜索下去，直接返回即可

关于$g_{i}$的计算，即$g_{i}=min_{炸弹j能匹配区间[i,k]}g_{k+1}+1$（初始$g_{m+1}=0$），复杂度为$o(n^{3})$

剪枝2：在搜索过程中，每一次划分出一个区间后，即在最终的二分图中新增一个点，根据二分图的过程，直接对其再求一次匹配，当某一次无法匹配直接退出即可

剪枝3：从后往前枚举右端点，在综合剪枝1和剪枝2的基础上，此优化也有一定的用处

#include<bits/stdc++.h>
using namespace std;
#define N 105
struct point{
    int x,y;
}a[N],b[N];
vector<int>v[N<<1];
int n,m,k,ans,g[N],Vis[N][N][N],mm[N][N<<1],vis[N<<1],match[N<<1];
int sqr(int k){
    return k*k;
}
bool check(point x,point y){
    return sqr(x.x-y.x)+sqr(x.y-y.y)<=sqr(k);
}
bool dfs(int k){
    if (vis[k])return 0;
    vis[k]=1;
    for(int i=0;i<v[k].size();i++)
        if ((!match[v[k][i]])||(dfs(match[v[k][i]]))){
            match[k]=v[k][i];
            match[v[k][i]]=k;
            return 1;
        }
    return 0;
} 
void dfs(int k,int s){
    if (s+g[k]>=ans)return;
    if (k>m){
        ans=s;
        return;
    }
    memcpy(mm[s],match,sizeof(mm[s]));
    for(int i=m;i>=k;i--){
        for(int j=1;j<=n;j++)
            if (Vis[j][k][i])v[n+s+1].push_back(j);
        memset(vis,0,sizeof(vis));
        if (dfs(n+s+1)){
            dfs(i+1,s+1);
            memcpy(match,mm[s],sizeof(mm[s]));
        }
        for(int j=1;j<=n;j++)
            if (Vis[j][k][i])v[n+s+1].pop_back();
    }
}
int main(){
    scanf("%d%d%d",&m,&n,&k);
    for(int i=1;i<=m;i++)scanf("%d%d",&a[i].x,&a[i].y);
    for(int i=1;i<=n;i++){
        scanf("%d%d",&b[i].x,&b[i].y);
        for(int j=1;j<=m;j++){
            if (!check(b[i],a[j]))continue;
            for(int k=j;k<=m+1;k++)
                if ((k>m)||(!check(b[i],a[k]))){
                    Vis[i][j][k-1]=1;
                    break;
                }
        }
    }
    g[m+1]=0;
    for(int i=m;i;i--){
        g[i]=n;
        for(int j=1;j<=n;j++)
            for(int k=i;k<=m;k++)
                if (Vis[j][i][k])g[i]=min(g[i],g[k+1]+1);
    }
    ans=n;
    dfs(1,0);
    printf("%d",ans);
}