2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛 题解

最新推荐文章于 2018-05-09 19:08:40 发布
最新推荐文章于 2018-05-09 19:08:40 发布
阅读量857 收藏
点赞数
分类专栏： 比赛题解 计算几何
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/nike0good/article/details/77927857
版权

B Out-out-control cars

给2个三角形，坐标均为整数，问它们分别以 (vx1,vy1),(vx2,vy2) 移动时，是否会相撞（撞一个点就行）

题解：考虑2个三角形相撞拆成线段相撞，进一步转化为其中一条线段的一个端点撞到另一条线上。把运动变为相对运动，于是变成了射线与线段相交的问题。
注意特判初始重合但相对静止的情况。
注意这题卡精度需要用__float128

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define ALL(x) (x).begin(),(x).end()
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
#define double __float128
ll sqr(ll a){return a*a;}
ld sqr(ld a){return a*a;}
double sqr(double a){return a*a;}
const __float128 eps=1e-20;
int dcmp(double x) {
    if (max(x,-x)<eps) return 0; else return x<0 ? -1 : 1; 
}

ld PI = 3.141592653589793238462643383;
class P{
public:
    double x,y;
    P(double x=0,double y=0):x(x),y(y){}
    friend ld dis2(P A,P B){return sqr(A.x-B.x)+sqr(A.y-B.y);   }
    friend ld Dot(P A,P B) {return A.x*B.x+A.y*B.y; }
    friend ld Length(P A) {return sqrt(Dot(A,A)); }
    friend ld Angle(P A,P B) {
        if (dcmp(Dot(A,A))==0||dcmp(Dot(B,B))==0||dcmp(Dot(A-B,A-B))==0) return 0;
        return acos(max((ld)-1.0, min((ld)1.0, Dot(A,B) / Length(A) / Length(B) )) ); 
    }

    friend P operator- (P A,P B) { return P(A.x-B.x,A.y-B.y); }
    friend P operator+ (P A,P B) { return P(A.x+B.x,A.y+B.y); }
    friend P operator* (P A,double p) { return P(A.x*p,A.y*p); }
    friend P operator/ (P A,double p) { return P(A.x/p,A.y/p); }
    friend bool operator< (const P& a,const P& b) {return dcmp(a.x-b.x)<0 ||(dcmp(a.x-b.x)==0&& dcmp(a.y-b.y)<0 );}

}; 
P read_point() {
    P a;
    ll c,d;
    c=read(),d=read();
    a.x=c,a.y=d;
    return a;   
} 
bool operator==(const P& a,const P& b) {
    return dcmp(a.x-b.x)==0 && dcmp(a.y-b.y) == 0;
} 
typedef P V;

double Cross(V A,V B) {return A.x*B.y - A.y*B.x;}
double Area2(P A,P B,P C) {return Cross(B-A,C-A);}
//Cross(v,w)==0(平行)时,不能调这个函数 
P GetLineIntersection(P p,V v,P Q,V w){
    V u = p-Q;
    double t = Cross(w,u)/Cross(v,w);
    return p+v*t;
}
//规范相交-线段相交且交点不在端点 
bool SegmentProperIntersection(P a1,P a2,P b1,P b2) { 
    double  c1 = Cross(a2-a1,b1-a1) , c2 = Cross(a2-a1,b2-a1),
            c3 = Cross(b2-b1,a1-b1) , c4 = Cross(b2-b1,a2-b1);
    return dcmp(c1)*dcmp(c2)<0 && dcmp(c3)*dcmp(c4)<0;
}

//点在线段上（不包含端点） 
bool OnSegment(P p,P a1,P a2) {
    return dcmp(Cross(a1-p,a2-p)) == 0 && dcmp(Dot(a1-p,a2-p))<0;
}


bool is_onshe(P a,V v,V q) {
//  cout<<Dot(q-a,v)<<endl;
//  cout<<dcmp(Dot(q-a,v))<<endl;
    return dcmp(Cross(q-a,v))==0 && dcmp(Dot(q-a,v))>=0;
}
P t1[4],t2[4],v1,v2;
bool check(P A,V b,P C,P D) {
    if (Length(b)<eps) {
        return A==C||A==D||OnSegment(A,C,D);
    }
    if (dcmp(Cross(D-C,b)==0)) {
        return is_onshe(A,b,C)||is_onshe(A,b,D);
    } else {
        if (is_onshe(A,b,C)||is_onshe(A,b,D)) return 1;
        if (SegmentProperIntersection(A,A+b*2000000000LL,C,D)) return 1;
        return 0;
        P q=GetLineIntersection(A,b,C,D-C);
        return is_onshe(A,b,q)&&OnSegment(q,C,D);
    }
}
typedef vector<P> Polygon ;
int isPointInPolygon(P p,Polygon poly) {
    int wn=0;
    int n=poly.size();
    Rep(i,n) {
        if (OnSegment(p,poly[i],poly[(i+1)%n])) return -1; //edge
        int k=dcmp(Cross(poly[(i+1)%n]-poly[i],p-poly[i]));
        int d1 = dcmp(poly[i].y-p.y);
        int d2 = dcmp(poly[(i+1)%n].y-p.y);
        if ( k > 0 && d1 <= 0 && d2 > 0 ) wn++;
        if ( k < 0 && d2 <= 0 && d1 > 0 ) wn--;
    }
    if (wn!=0) return 1; //inside
    return 0; //outside
}

int main()
{
//  freopen("B.in","r",stdin);
//  freopen(".out","w",stdout);
    int T=read();
    For(kcase,T) {
        Rep(i,3) t1[i]=read_point();v1=read_point();
        Rep(i,3) t2[i]=read_point();v2=read_point();
        bool fl=0;
        if (v1==v2) {
            Polygon poly1,poly2;
            Rep(i,3) poly1.pb(t1[i]);
            Rep(j,3) {
                if (isPointInPolygon(t2[j],poly1)) fl=1;
            }
            Rep(i,3) poly2.pb(t2[i]);
            Rep(j,3) {
                if (isPointInPolygon(t1[j],poly2)) fl=1;
            }
        }
        Rep(i,3) {
            P C=t1[i],D=t1[(i+1)%3];
            Rep(j,3) {
                P A=t2[j];
                V b=v2-v1;
                if (check(A,b,C,D)) {
                    fl=1;
                }
            }
        }
        Rep(i,3) {
            P C=t2[i],D=t2[(i+1)%3];
            Rep(j,3) {
                P A=t1[j];
                V b=v1-v2;
                if (check(A,b,C,D)) fl=1;
            }
        }
        printf("Case #%d: ",kcase);
        puts(fl?"YES":"NO");
    }

    return 0;
}

D Hack Portals

有一堆Portal在一条数轴上，每个Portal有一个冷却时间 ci ,位置 xi ，现在你需要hack所有的Portal（一个Portal必须在冷却时间结束后hack），你的起点为0，终点为W，求最小花费时间。

根据多年玩ingress的经验，
先hack一个中间的Portal以后再走到左边hack再走到右边hack，中间Portal最后一次经过的时间会被更新，因此可以假定在任意时刻没被hack的Portal一定是数轴上连续一段区间，区间dp。

#include<bits/stdc++.h> 
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,0x3f,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define MEMx(a,b) memset(a,b,sizeof(a));
#define INF (0x3f3f3f3f)
#define F (1000000007)
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %lld\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
#define ALL(x) (x).begin(),(x).end()
#define gmax(a,b) a=max(a,b);
#define gmin(a,b) a=min(a,b);
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
inline int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
#define MAXN (1010)
int n,L,K;
ll f[MAXN][MAXN][2];
ll Abs(ll x) {
    return max(x,-x);
}
ll x[MAXN],c[MAXN];
pi p[MAXN];
int main()
{
//  freopen("D.in","r",stdin);
//  freopen(".out","w",stdout);
    int T=read();
    For(kcase,T) {
        ll ans=INF;
        n=read(),L=read(),K=read();
        For(i,n) p[i].fi=read(),p[i].se=read();
        sort(p+1,p+1+n);
        For(i,n) x[i]=p[i].fi,c[i]=p[i].se;

        memset(f,127/2,sizeof(f));
        f[1][n][0]=max(x[1],c[1]);
        f[1][n][1]=max(x[n],c[n]);

        ForD(len,n) if (len>1){
            For(i,n-len+1) {
                int j=i+len-1;
                gmin(f[i+1][j][0],max(c[i+1],f[i][j][0]+Abs(x[i]-x[i+1])));
                gmin(f[i+1][j][1],max(c[j],  f[i][j][0]+Abs(x[i]-x[j])));

                gmin(f[i][j-1][0],max(c[i],f[i][j][1]+Abs(x[j]-x[i])));
                gmin(f[i][j-1][1],max(c[j-1],f[i][j][1]+Abs(x[j]-x[j-1])));

            }
        }
        For(i,n) {

(x[i]-K) );
        }
        Pr(kcase,ans);
    }

    return 0;
}

E Half-consecutive Numbers

ci=n(n+1)/2,i=1,2,⋯ ,已知 n ，求最小的i使得i≥n，且 ci 是完全平方数。n<=1e16
考虑 ci 是完全平方数，那么i或i+1中存在一个奇数，这个数也是完全平方数，那么可以在10^8之内暴力枚举满足条件的i，打表。

#include<bits/stdc++.h> 
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,0x3f,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define MEMx(a,b) memset(a,b,sizeof(a));
#define INF (0x3f3f3f3f)
#define F (1000000007)
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %lld\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
#define ALL(x) (x).begin(),(x).end()
#define gmax(a,b) a=max(a,b);
#define gmin(a,b) a=min(a,b);
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
inline int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
ll a[123456];
ll is_sq(ll i) {
    ll t=sqrt(i);t=t*t;
    if (t==i) return 1;return 0;
}
ll check(ll i) {
    ll t=i*(i+1)/2,p=sqrt(t);
    if (p*p==t) return 1;return 0;
}
ll ans[1234]={0,1,8,49,288,1681,9800,57121,332928,1940449,11309768,65918161,384199200,2239277041,13051463048,76069501249,443365544448,2584123765441,15061377048200,87784138523761,511643454094368,2982076586042449,17380816062160328};
int main()
{
//  freopen("E.in","r",stdin);
//  freopen("e.out","w",stdout);
    int T=read();


    For(kcase,T) {
        ll n;
        scanf("%lld",&n);
        int i=1;
        while(ans[i]<n) ++i;
        printf("Case #%d: %lld\n",kcase,ans[i]);
    }
    return 0;
}
nike0good
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
ACM-ICPC 历年竞赛 真题，各大赛区真题详解
04-24
ACM-ICPC分为多个赛区，包括但不限于亚洲、欧洲、北美、南美等。每个赛区的竞赛题目都反映了该地区在计算机科学教育上的特色和优势。例如，亚洲赛区的题目可能更注重算法的创新和实用性，而北美赛区的题目则可能更加...
ACM-ACM-ICPC算法示例之Greedy-题解.zip
最新发布
02-25
本资料包“ACM_ACM-ICPC算法示例之Greedy_题解.zip”主要包含了关于贪心算法的实例解析和题解，旨在帮助参赛者理解和掌握这种算法的应用。 贪心算法的核心思想是在每一步决策时都采取在当前状态下最好或最优的选择...
H. Skiing 2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛
axuhongbo的博客
09-09 709
2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛
计蒜客-2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛-H-Skiing
逐梦者
09-11 353
ACM模版题解我不喜欢计蒜客，我都等了两天了，还没有地方可以提交……现在连题目都看不了，算了，题解就不写了，光发一下代码吧……代码#include <iostream> #include <queue> #include <cstring>using namespace std;const int MAXN = 1e4 + 7; const int MAXM = 2e5 + 7;int n, m,
2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛 G题 Query on a string
t51645的博客
09-09 1356
题意： 给出字符串S和字符串T 然后q次查询 查询分为两种： 1:查询字符串区间[x,y]内T能被匹配几次 2.修改字符串S[x]为c 题解： 树状数组+瞎搞 用树状数组保存前k个位置能匹配几次T 第i个位置能匹配的意思是 S[i~i+len-1] = T[1~len]  len是字符串T的长度 首先预处理出树状数组，然后查询时直接给出答案，修改时先将字符串修改
计蒜客 2017 ACM-ICPC 亚洲区（南宁赛区）网络赛 L题
eletroram的博客
09-26 267
Let SS be a sequence of integers s_{1}s​1​​, s_{2}s​2​​, ......, s_{n}s​n​​ Each integer is is associated with a weight by the following rules: (1) If is is negative, then its weight is 00. (2)
2017 ACM-ICPC 亚洲区（西安赛区）网络赛 C.Sum(找规律)
HJY
09-16 754
Define the function S(x)S(x) for xx is a positive integer. S(x)S(x) equals to the sum of all digit of the decimal expression of xx. Please find a positive integer kk that S(k*x)\%233=0S(k∗x)%233=0.
2017 ACM-ICPC 亚洲区(西安赛区)网络赛 E Maximum Flow
binarycopycode
09-19 636
原题解链接：http://blog.csdn.net/kkkkahlua/article/details/78009087 他用的最小割的来求解最大流。认为只要讨论每一个点到0和n-1那个更小哪条边就一定是最小割中的一条边。感觉很有道理，然而不知道是怎么证的。 然而在考场上我找规律做的。。。
A. Banana 2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛
axuhongbo的博客
09-09 1183
2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛
C. Coconut 2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛
axuhongbo的博客
09-09 253
2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛
2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛 A. Banana（连通性水题）
HJY
09-09 405
Bananas are the favoured food of monkeys. In the forest, there is a Banana Company that provides bananas from different places. The company has two lists. The first list records the types of bananas preferred by different monkeys, and the second one rec
F. Islands 2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛
axuhongbo的博客
09-09 281
2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛
ACM-ACM-ICPC算法示例之Hashing-题解.zip
02-25
在ACM（国际大学生程序设计竞赛）和ICPC（国际编程竞赛）中，算法的熟练掌握是关键。Hashing作为一种高效的数据结构和算法技术，在解决许多问题时起着至关重要的作用。本压缩包中的资源主要提供了Hashing算法的实例...
基于C++的ACM-ICPC模板
02-22
总的来说，这个基于C++的ACM-ICPC模板是参赛者在备赛过程中必不可少的工具，它可以帮助参赛者专注于算法设计，而不需要花费过多精力在代码实现的细节上。通过学习和使用这套模板，选手可以提升自己的编程效率，提高...
求面积 （坐标叉积公式+凹多边形面积-坐标公式）
nike0good |Oier&ACMer | 熟能生巧
10-26 3530
求面积（AREA） 给出一个简单多边形(没有缺口)，它的边要么是垂直的，要么是水平的。要求计算多边形的面积。 多边形被放置在一个X-Y的卡笛尔平面上，它所有的边都平行于两条坐标轴之一。然后按逆时针方向给出各顶点的坐标值。所有的坐标值都是整数(因此多边形的面积也为整数)。   输入 输入文件第一行给出多边形的顶点数n(n≤100)。接下来的几行每行给出多边形一个顶点的坐标值X和Y(都为整数
JSOI2018(战争-Minkowski sum)
nike0good |Oier&ACMer | 熟能生巧
05-09 1989
给你2个点集，分别对这2个点集求凸包A,B，每次询问问凸包B平移向量(dx.dy)(dx.dy)后是否与凸包A有公共点。点集大小n,m，询问次数q的范围为1e5 b⃗ +(dx,dy)=a⃗ \vec{b}+(dx,dy)=\vec{a} (dx,dy)=a⃗ −b⃗ (dx,dy)=\vec{a}-\vec{b} 因此合法的(dx,dy)(dx,dy)范围为：A和-B的Minkowski
POJ 2187(凸包GrahamScan扫描+极角排序+平面最远点对）
nike0good |Oier&ACMer | 熟能生巧
01-18 1676
Language: Default 平面最远点对 Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 22454   Accepted: 6868 Description 求 N (2 的平方.  Input * 第1行:N  接下来每行为
BZOJ 1336([Balkan2002]Alien最小圆覆盖-最小覆盖圆)
nike0good |Oier&ACMer | 熟能生巧
07-22 1670
1336: [Balkan2002]Alien最小圆覆盖 Time Limit: 1 Sec  Memory Limit: 162 MBSec  Special Judge Submit: 733  Solved: 360 [Submit][Status] Description 给出N个点，让你画一个最小的包含所有点的圆。 Input 先给出点的个数N,2 Out
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值