UVA1392 DNA Regions

链接

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=4138

题解

这题数据水啊
构建一个$01$序列，如果两个序列的第$i$个位置不同，则是$1$，否则是$0$，这个序列的前$i$项和记为$s_i$
首先把$(i,s_i)$看成是坐标系里面的一些点，这个题就是让你找两个点，使得它们连线的斜率小于等于$\frac{p}{100}$
而且这两个点的横坐标之差尽量大
那我就直接二分法答案L，现在就是判断是否存在横坐标之差大于等于$L$且连线不大于$\frac{p}{100}$的点对
可以先枚举横坐标比较大的点$i$，然后把$1...i-L$这些点里面选出一个与$i$点连线斜率最小的点，比较这个值和$\frac{p}{100}$的大小，我试了一下直接建成凸壳然后暴力扫描是可以$AC$的（所以说数据水），但是如果使用三分复杂度则更稳定些，大概$O(nlog_{1.5}nlog_{2}n)$
这个题得到一个很好的思路就是，题目有时候直接考虑的话，答案不满足二分性质，但是如果把$check$的标准改为：答案如果大于等于这个数，能否行得通，就满足二分性质了

代码

//二分、三分、凸壳 
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cctype>
#define maxn 150010
#define cl(x) memset(x,0,sizeof(x))
#define ll long long
#define eps 1e-8
#define inf 1e60
using namespace std;
int s[maxn], n, i, tot;
double p;
char A[maxn], B[maxn];
struct point{double x, y;}pt[maxn], hull[maxn];
double k(point &p1, point &p2){return (p2.y-p1.y)/(p2.x-p1.x);}
bool init()
{
    scanf("%d%lf",&n,&p);
    if(n==0)return false;
    p/=100;
    scanf("%s%s",A+1,B+1);
    for(i=1;i<=n;i++)s[i]=(A[i]!=B[i])+s[i-1], pt[i].x=i, pt[i].y=s[i];//, printf("(%.2lf,%.2lf)\n",pt[i].x,pt[i].y);
    return true;
}
bool check(int L)
{
    int i, j, l, r, mid1, mid2;
    double t, ans=inf;
    tot=0;
    for(i=L;i<=n;i++)
    {
        for(;tot>=2 and k(pt[i-L],hull[tot])>k(hull[tot],hull[tot-1])-eps;tot--);
        hull[++tot]=pt[i-L];
        for(l=1,r=tot,mid1=l+(r-l)/3,mid2=r-(r-l)/3;r-l>=3;mid1=l+(r-l)/3,mid2=r-(r-l)/3)
        {
            if(k(pt[i],hull[mid1])<k(pt[i],hull[mid2])-eps)r=mid2;
            else l=mid1;
        }
        t=k(pt[i],hull[l]);
        for(j=l+1;j<=r;j++)t=min(t,k(pt[i],hull[j]));
        if(t<p+eps)return true;
    }
    return false;
}
void work()
{
    int l, r, mid;
    l=1, r=n, mid=(l+r+1)>>1;
    if(s[n]==n)
    {
        printf("No solution.\n");
        return;
    }
    while(l^r)
    {
        if(check(mid))l=mid;
        else r=mid-1;
        mid=(l+r+1)>>1;
    }
    printf("%d\n",l);
}
int main()
{
    while(init())work();
    return 0;
}