BZOJ 3190: [JLOI2013]赛车

最新推荐文章于 2020-01-13 20:23:43 发布

shjzhqm

最新推荐文章于 2020-01-13 20:23:43 发布

阅读量118

点赞数 2

分类专栏：计算几何

本文链接：https://blog.csdn.net/shjzhqm/article/details/79966526

版权

计算几何专栏收录该内容

2 篇文章 0 订阅

订阅专栏

这题被卡了好久。。。
令我震惊的是暴力可过。。。
这道题可以先把每个赛车看做一个s-t图像，斜率就是v
那么这题就和1007水平可见直线没啥区别了（其实多了三线共点~~~，而且是在第一象限内）

暴力

#include<cmath>
#include<vector> 
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define eps 1e-8
#define inf 0x7fffffff
using namespace std;
struct p
{
    int k,b;
}line[50005];
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct q
{
    double l,r;
};
int n;
int stack[50005],top,tot;
int ans[50005];
int main()
{
    int j,i;
    bool flag;
    n=read();
    for(i=1;i<=n;i++) line[i].b=read();
    for(i=1;i<=n;i++) line[i].k=read();
    for(i=1;i<=n;i++)
    {
        flag=1;
        q x;
        x.l=0,x.r=inf;
        for(j=1;j<=n;j++)
        {
            if(i==j) continue;
            if(line[i].k==line[j].k&&line[j].b>line[i].b) {flag=0;break;}
            if(line[i].k<line[j].k&&line[i].b<line[j].b) {flag=0;break;}
            q y;
            y.l=(double)(line[j].b-line[i].b)/(double)(line[i].k-line[j].k);
            if(line[i].k>line[j].k) x.l=max(x.l,y.l);
            if(line[i].k<line[j].k) x.r=min(x.r,y.l);
            if(x.r-x.l<-eps) {flag=0;break;}
        }
        if(flag) ans[++ans[0]]=i;
    }
    printf("%d\n",ans[0]);
    for(i=1;i<ans[0];i++)
        printf("%d ",ans[i]);
    if(ans[0]) 
        printf("%d",ans[ans[0]]);
    return 0;
}

正解

#include<cmath>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define inf 0x7fffffff
#define eps 1e-10
using namespace std;
struct p
{
    double k,b;
    int id,num;
    bool operator < (const p &a)const
    {
        return k<a.k||(fabs(k-a.k)<eps&&b<a.b);
    }
}line[50005];
struct q
{
    double x,y;
}node[50005],s;
int tot;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n;
int stack[50005],top;
bool cmp(int a,int b)
{
    return line[a].id<line[b].id;
}
double find(int x,int y)
{
    double xx;
    xx=(line[y].b-line[x].b)/(line[x].k-line[y].k);
    if(xx<0) return -inf;
    return xx;
}
bool ans[50005];
double ans1,ans2;
vector<int> bel[50005];
int main()
{
    int i,j,now,cnt=0,sum=0;
    n=read();
    for(i=1;i<=n;i++) cin>>line[i].b,line[i].id=i;
    for(i=1;i<=n;i++) cin>>line[i].k;
    sort(line+1,line+n+1);
    for(i=1;i<=n;i++)
    {
        if(i!=1&&line[i].k==line[i-1].k&&line[i].b==line[i-1].b) line[tot].num++,bel[tot].push_back(line[i].id);
        else line[++tot]=line[i],line[tot].num=1,bel[tot].push_back(line[i].id);
    }
    for(i=1;i<=tot;i++)
    {
        while(top)
        {
            if(fabs(line[stack[top]].k-line[i].k)<eps)top--;
            else if(find(i,stack[top])<-eps)top--;
            else
            {
                ans1=find(i,stack[top-1]);
                ans2=find(stack[top],stack[top-1]);
                if(top&&ans1<ans2-eps)
                    top--;
                else break;
            }
        }
        stack[++top]=i;
    }
    for(i=1;i<=top;i++) 
    {
        for(j=0;j<bel[stack[i]].size();j++)
        {
            ans[bel[stack[i]][j]]=1;
            sum++;
        }
    }
    printf("%d\n",sum);
    for(i=1;i<=n;i++)
    {
        if(ans[i]) 
        {
            cnt++;
            printf("%d",i);
            if(cnt!=sum) printf(" ");
            else return 0;
        }
    }
    return 0;
}