Gym-101673 :East Central North America Regional Contest (ECNA 2017)(寒假自训第8场)

A .Abstract Art

题意：求多个多边形的面积并。

思路：模板题。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const double inf=1e200;
const double eps=1e-;
const double pi=*atan(1.0);
int dcmp(double x){ return fabs(x)<eps?:(x<?-:);}
struct point{
    double x,y;
    point(double a=,double b=):x(a),y(b){}
};
point operator +(point A,point B) { return point(A.x+B.x,A.y+B.y);}
point operator -(point A,point B) { return point(A.x-B.x,A.y-B.y);}
point operator *(point A,double p){ return point(A.x*p,A.y*p);}
point operator /(point A,double p){ return point(A.x/p,A.y/p);}
bool operator ==(const point& a,const point& b){
    return fabs(a.x-b.x)<eps&&fabs(a.y-b.y)<eps;
}
double dot(point A,point B){ return A.x*B.x+A.y*B.y;}
double det(point A,point B){ return A.x*B.y-A.y*B.x;}
double det(point O,point A,point B){ return det(A-O,B-O);}
double length(point A){ return sqrt(dot(A,A));}
double area(vector<point>p){
    double ans=; int sz=p.size();
    for(int i=;i<sz-;i++) ans+=det(p[i]-p[],p[i+]-p[]);
    return ans/2.0;
}
double seg(point O,point A,point B){
    if(dcmp(B.x-A.x)==) return (O.y-A.y)/(B.y-A.y);
    return (O.x-A.x)/(B.x-A.x);
}
vector<point>pp[];
pair<double,int>s[*];
double polyunion(vector<point>*p,int N){
    double res=;
    for(int i=;i<N;i++){
        int sz=p[i].size();
        for(int j=;j<sz;j++){
            int m=;
            s[m++]=make_pair(,);
            s[m++]=make_pair(,);
            point a=p[i][j],b=p[i][(j+)%sz];
            for(int k=;k<N;k++){
                if(i!=k){
                    int sz2=p[k].size();
                    for(int ii=;ii<sz2;ii++){
                        point c=p[k][ii],d=p[k][(ii+)%sz2];
                        int c1=dcmp(det(b-a,c-a));
                        int c2=dcmp(det(b-a,d-a));
                        if(c1==&&c2==){
                            if(dcmp(dot(b-a,d-c))){
                                s[m++]=make_pair(seg(c,a,b),);
                                s[m++]=make_pair(seg(c,a,b),-);
                            }
                        }
                        else{
                            double s1=det(d-c,a-c);
                            double s2=det(d-c,b-c);
                            if(c1>=&&c2<) s[m++]=make_pair(s1/(s1-s2),);
                            else if(c1<&&c2>=) s[m++]=make_pair(s1/(s1-s2),-);
                        }
                    }
                }
            }
            sort(s,s+m);
            double pre=min(max(s[].first,0.0),1.0),now,sum=;
            int cov=s[].second;
            for(int j=;j<m;j++){
                now=min(max(s[j].first,0.0),1.0);
                if(!cov) sum+=now-pre;
                cov+=s[j].second;
                pre=now;
            }
            res+=det(a,b)*sum;
        }
    }
    return res/;
}
int main()
{
    int N,M,i,j; point tp;
    scanf("%d",&N);
    for(i=;i<N;i++){
        scanf("%d",&M);
        for(j=;j<M;j++){
            scanf("%lf%lf",&tp.x,&tp.y);
            pp[i].push_back(tp);
        }
    }
    double t1=,t2=polyunion(pp,N);
    for(i=;i<N;i++) t1+=area(pp[i]);
    printf("%f %f\n",-t1,-t2);
    return ;
}

B .Craters

题意：给定N个圆，让你用最小的周长把这些圆包起来，且满足到圆的最近距离不小于10.

思路：把每个圆的半径增加10，然后等分1000份，然后求凸包即可。

#include<bits/stdc++.h>
#define mp make_pair
using namespace std;
typedef long long ll;
const double inf=1e200;
const double eps=1e-;
const double pi=*atan(1.0);
int dcmp(double x){ return fabs(x)<eps?:(x<?-:);}
struct point{
    double x,y;
    point(double a=,double b=):x(a),y(b){}
};
point operator +(point A,point B) { return point(A.x+B.x,A.y+B.y);}
point operator -(point A,point B) { return point(A.x-B.x,A.y-B.y);}
point operator *(point A,double p){ return point(A.x*p,A.y*p);}
point operator /(point A,double p){ return point(A.x/p,A.y/p);}
point rotate(point A,double rad){
    return point(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));
}
bool operator ==(const point& a,const point& b) {
     return dcmp(a.x-b.x)==&&dcmp(a.y-b.y)==;
}
double dot(point A,point B){ return A.x*B.x+A.y*B.y;}
double det(point A,point B){ return A.x*B.y-A.y*B.x;}
double det(point O,point A,point B){ return det(A-O,B-O);}
double length(point A){ return sqrt(dot(A,A));}
double angle(point A,point B){ return acos(dot(A,B)/length(A)/length(B));}
double area(vector<point>p){
    double ans=; int sz=p.size();
    for(int i=;i<sz-;i++) ans+=det(p[i]-p[],p[i+]-p[]);
    return ans/2.0;
}
double seg(point O,point A,point B){
    if(dcmp(B.x-A.x)==) return (O.y-A.y)/(B.y-A.y);
    return (O.x-A.x)/(B.x-A.x);
}
bool cmp(point a,point b){ return a.x==b.x?a.y<b.y:a.x<b.x; }
void convexhull(point *a,int n,point *ch,int &top)
{
    sort(a+,a+n+,cmp);
    top=;
    for(int i=;i<=n;i++){
        while(top>=&&det(ch[top-],ch[top],a[i])<=) top--;
        ch[++top]=a[i];
    }
    int ttop=top;
    for(int i=n-;i>=;i--){
        while(top>ttop&&det(ch[top-],ch[top],a[i])<=) top--;
        ch[++top]=a[i];
    }
}
point ch[],p[];
int main()
{
    int N,i,j,tot=,top; double ans,one=pi/,x,y,r;
    scanf("%d",&N);
    for(i=;i<=N;i++){
        scanf("%lf%lf%lf",&x,&y,&r); r+=;
        for(j=;j<=;j++){
            tot++; p[tot].x=x+r*cos(one*j); p[tot].y=y+r*sin(one*j);
        }
    }
    convexhull(p,tot,ch,top);
    for(i=;i<top;i++) ans+=length(ch[i]-ch[i+]);
    printf("%.10lf\n",ans);
    return ;
}

C .DRM Messages

题意：字符串操作模拟。

思路：模拟。

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
char c[maxn];
void solve(int L,int R)
{
    int sum=;
    rep(i,L,R) sum+=c[i]-'A';
    sum%=;
    rep(i,L,R) c[i]=(c[i]-'A'+sum+)%+'A';
}
int main()
{
    int N; scanf("%s",c+);
    N=strlen(c+);
    solve(,N/);
    solve(N/+,N);
    rep(i,,N/){
        c[i]=(c[i]-'A'+(c[i+N/]-'A')+)%+'A';
    }
    rep(i,,N/) putchar(c[i]);
    return ;
}

D .Game of Throwns

题意：模拟题。

思路：模拟。

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
int q[maxn],tot;
int read()
{
    int x=,F=; char c=getchar();
    while(c!='-'&&c!='u'&&!(c>=''&&c<='')) c=getchar();
    if(c=='-') F=-,c=getchar();
    else if(c=='u') {
        rep(i,,) c=getchar(); return maxn;
    }
    while(c>=''&&c<='') x=x*+c-'',c=getchar();
    return F*x;
}
int main()
{
    int N,K,x,P=; scanf("%d%d",&N,&K);
    rep(i,,K) {
        x=read();
        if(x==maxn)  {
            x=read();
            tot-=x;
        }
        else q[++tot]=x;
    }
    rep(i,,tot) P=((P+q[i])%N+N)%N;
    printf("%d\n",P);
    return ;
}

E .Is-A? Has-A? Who Knowz-A?

题意：题意可能没有解释清楚。N个名字，M个传递关系，Q次询问，反正就是有两种传递关系，is和has，A（has or is）->B; B->C; C->D;像这样的传递关系，假如全部都是is，则A可以推出isD；否则可以推出A has D。N<=500。

思路：因为N<500，直接DFS即可，复杂度O(N^3)；注意自己和自己有is关系，但是没有has关系。

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
int dis[maxn][maxn][],tot;
struct in{
    int x,y,opt;
    in(){}
    in(int xx,int yy,int oo):x(xx),y(yy),opt(oo){}
};
queue<in>q; map<string,int>mp;
map<int,string>F;
void name(string s){
    if(mp.find(s)==mp.end()) {
            mp[s]=++tot; F[tot]=s;
    }
}
void DFS()
{
    while(!q.empty()){
        int x=q.front().x,y=q.front().y,op=q.front().opt; q.pop();
        //cout<<":"<<F[x]<<" "<<F[y]<<" "<<op<<endl;
        dis[x][y][op]=;
        for(int i=;i<=tot;i++){
            if(!dis[x][i][]&&dis[y][i][]){
                dis[x][i][]=; q.push(in(x,i,));
            }
            if(!dis[x][i][]&&dis[y][i][]&&op==){
                dis[x][i][]=; q.push(in(x,i,));
            }
            if(!dis[x][i][]&&dis[y][i][]&&op==){
                dis[x][i][]=; q.push(in(x,i,));
            }
        }
    }
}
int main()
{
    int N,M,u,v,Ca=; string A,B,C;
    scanf("%d%d",&N,&M);
    rep(i,,N){
        cin>>A>>B>>C;
        name(A); name(C);
        if(B[]=='i') q.push(in(mp[A],mp[C],)),dis[mp[A]][mp[C]][]=;
        else q.push(in(mp[A],mp[C],)),dis[mp[A]][mp[C]][]=;
    }
    DFS();
    rep(i,,M){
        cin>>A>>B>>C;
        printf("Query %d: ",++Ca);
        if(A==C&&B[]=='i') {puts("true"); continue;}
        u=mp[A]; v=mp[C];
        if(B[]=='i') {
            if(dis[u][v][]) puts("true");
            else puts("false");
        }
        else {
            if(dis[u][v][]) puts("true");
            else puts("false");
        }
    }
    return ;
}

F .Keeping On Track

题意：给定一棵树，问删去一个点后，不连通的点最对有多少，在此基础上加一条边可以最多让多少对点对恢复连通。

思路：就是dfs一次，然后看所有儿子和父亲之上的连通块的大小； 第二问取最大的两个连通块连通即可。

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
int Laxt[maxn],Next[maxn],To[maxn];
int sz[maxn],q[maxn],cnt,ans,fcy,tot,N;
void add(int u,int v)
{
    Next[++cnt]=Laxt[u];Laxt[u]=cnt;To[cnt]=v;
}
void dfs(int u,int f)
{
    sz[u]=;
    for(int i=Laxt[u];i;i=Next[i]){
        if(To[i]!=f){
            dfs(To[i],u);
            sz[u]+=sz[To[i]];
        }
    }
    tot=; q[++tot]=N+-sz[u];
    for(int i=Laxt[u];i;i=Next[i])
      if(To[i]!=f) q[++tot]=sz[To[i]];
    sort(q+,q+tot+);
    int sum=,tmp=;
    rep(i,,tot) sum+=q[i];
    rep(i,,tot) tmp+=(sum-q[i])*q[i];
    tmp/=;
    if(tmp>ans) ans=tmp,fcy=q[tot]*q[tot-];
    else if(tmp==ans) fcy=max(fcy,q[tot]*q[tot-]);
}
int main()
{
    int u,v;
    scanf("%d",&N);
    rep(i,,N) {
        scanf("%d%d",&u,&v);
        add(u,v); add(v,u);
    }
    dfs(,-);
    printf("%d %d\n",ans,ans-fcy);
    return ;
}

G .A Question of Ingestion

题意：一个人开始的饭量是M，如果今天吃饭了，第二天的饭量变为今天的2/3；如果今天不吃，第二天的饭量和昨天饭量一样（即恢复）； 如果连续两天不吃，则第三天的饭量恢复到M。  现在给出N天的食物量， 问怎么吃可以迟到最多的食物。

思路：DP即可，dp[i][j]表示第i天的时候是连续第j天吃饭的最大食物量。

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
int a[maxn],f[maxn],dp[maxn][maxn],ans;
int main()
{
    int N,M;
    scanf("%d%d",&N,&M); f[]=M;
    rep(i,,N) scanf("%d",&a[i]);
    rep(i,,N) f[i]=f[i-]*/;
    rep(i,,N){
        rep(j,,i){
            dp[i][j]=max(dp[i][j],dp[i-][j-]+min(f[j],a[i]));
            if(i>)dp[i][j]=max(dp[i][j],dp[i-][j]+min(f[j],a[i]));
            if(i>)dp[i][]=max(dp[i][],dp[i-][j]+min(f[],a[i]));
        }
    }
    rep(i,,N) ans=max(ans,dp[N][i]);
    printf("%d\n",ans);
    return ;
}

H .Sheba's Amoebas

题意：问有几个连通块。

思路：DFS或者并查集即可。

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
char c[maxn][maxn];
int fa[maxn*maxn],ans,N,M;
int x[]={,,,-,-,-,,};
int y[]={-,,,-,,,-,};
int find(int x){
    if(x==fa[x]) return x;
    return fa[x]=find(fa[x]);
}
void merge(int u,int v)
{
    int fu=find(u),fv=find(v);
    if(fu!=fv) fa[fu]=fv;
}
int main()
{
    scanf("%d%d",&N,&M);
    rep(i,,N*M) fa[i]=i;
    rep(i,,N) scanf("%s",c[i]+);
    rep(i,,N)
     rep(j,,M){
        if(c[i][j]=='#'){
            rep(k,,){
                if(i+x[k]>=&&i+x[k]<=N&&j+y[k]>=&&j+y[k]<=M&&c[i+x[k]][j+y[k]]=='#'){
                    merge((i-)*M+j,(i+x[k]-)*M+j+y[k]);
                }
            }
        }
    }
    rep(i,,N)
     rep(j,,M){
       if(c[i][j]!='#') continue;
       int t=(i-)*M+j;
       if(find(t)==t) ans++;
    }
    printf("%d\n",ans);
    return ;
}

太累了，后面了两个题先鸽了