• CF1143F/1142C U2


    CF1143F/1142C U2

    • 巧妙的思维题.注意到这里只用两个点就可以确定一根抛物线,联想到两点确定一条直线,尝试转化.
    • (y=x^2+bx+c) 就可以写成 (y-x^2=bx+c) ,可以发现,将点 ((x_i,y_i)) 变为 ((x_i,y_i-x_i^2)) 后,就变成了对每对点连一条直线,答案就是上凸壳的边数.
    • 注意在统计答案时,不能计算两点 (x) 相同的直线.
    #include<bits/stdc++.h>
    using namespace std;
    #define ll long long
    #define mp make_pair
    #define pii pair<int,int>
    inline int read()
    {
    	int x=0;
    	bool pos=1;
    	char ch=getchar();
    	for(;!isdigit(ch);ch=getchar())
    		if(ch=='-')
    			pos=0;
    	for(;isdigit(ch);ch=getchar())
    		x=x*10+ch-'0';
    	return pos?x:-x;
    }
    const int MAXN=1e5+10;
    int n;
    struct v2
    {
    	double x,y;
    	v2(double x=0,double y=0):x(x),y(y) {}
    	friend double operator * (const v2 &a,const v2 &b)
    		{
    			return a.x*b.y-a.y*b.x;
    		}
    	friend v2 operator + (const v2 &a,const v2 &b)
    		{
    			return v2(a.x+b.x,a.y+b.y);
    		}
    	friend v2 operator - (const v2 &a,const v2 &b)
    		{
    			return v2(a.x-b.x,a.y-b.y);
    		}
    	bool operator < (const v2 &rhs) const
    		{
    			return x==rhs.x?y<rhs.y:x<rhs.x;
    		}
    	double modulus()
    		{
    			return sqrt(x*x+y*y);
    		}
    	double angle()
    		{
    			return atan2(y,x);
    		}
    }p[MAXN];
    v2 stk[MAXN];
    int tp=0;
    v2 origin;
    bool cmp(const v2 &a,const v2 &b)
    {
    	double a1=(a-origin).angle();
    	double a2=(b-origin).angle();
    	return a1==a2?a.x<b.x:a1<a2;
    }
    void ConvexHull()
    {
    	for(int i=2;i<=n;++i)
    		if(p[i]<p[1])
    			swap(p[i],p[1]);
    	origin=p[1];
    	sort(p+2,p+1+n,cmp);
    	for(int i=1;i<=n;++i)
    	{
    		while(tp>=2 && (stk[tp]-stk[tp-1])*(p[i]-stk[tp])<=0)
    			--tp;
    		stk[++tp]=p[i];
    	}
    	stk[++tp]=p[1];
    }
    #define eps 1e-12
    int main()
    {
    	n=read();
    	for(int i=1;i<=n;++i)
    	{
    		int x=read(),y=read();
    		p[i]=v2((double)x,(double)(y-1LL*x*x));
    	}
    	ConvexHull();
    	if(tp==3)
    	{
    		if(fabs(stk[1].x-stk[2].x)<eps)
    			cout<<0<<endl;
    		else
    			cout<<1<<endl;
    		return 0;
    	}
    	int ans=0;
    	for(int i=1;i<tp;++i)
    	{
    		int j=(i==1)?3:i-1;
    		if(fabs(stk[i+1].x-stk[i].x)<eps)
    			continue;
    		double y=stk[i].y+(stk[j].x-stk[i].x)/(stk[i+1].x-stk[i].x)*(stk[i+1].y-stk[i].y);
    		if(y-stk[j].y>eps)
    			++ans;
    	}
    	cout<<ans<<endl;
    	return 0;
    }
    
  • 相关阅读:
    ZooKeeper 授权访问
    jQuery代码优化:事件委托篇
    设计一个学生作业完成情况的管理程序
    UVa 10192 Vacation(LCS水题)
    BBS快照源代码代写
    ftp服务器系统编程编码程序下载代写代开发
    文件夹病毒专杀工具
    C#生成软件注册码
    《恐怖丛林生存》肉搏攻略,解开星星系统之谜
    《恐怖丛林生存》各种bug 各种解密
  • 原文地址:https://www.cnblogs.com/jklover/p/10660190.html
Copyright © 2020-2023  润新知