2019河北省大学生程序设计竞赛部分题题解

2019河北省大学生程序设计竞赛

---

A. Battle of Balls

Battle of Balls

即球心到障碍的距离要大于半径。把球变成点，障碍变成圆，找割就行了。
所以是个连通性问题。我喜欢并查集，有人喜欢遍历。
浮点数运算，数都放到不等式左边，不等式右边用eps控制精度。

#include<bits/stdc++.h>
using namespace std;
const int N=1010;
const double eps=1e-8;
int f[N],T,n;
double r,x[N],y[N];

int find(int x)
{
	if(f[x]==x) return f[x];
	else return f[x]=find(f[x]);
}

void merge(int x,int y)
{
	int fx=find(x),fy=find(y);
	if(fx!=fy) f[fy]=fx;
}

double dis(int i,int j)
{
	return sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
}

int main()
{
	cin>>T;
	while(T--)
	{
		scanf("%d%lf",&n,&r);
		for(int i=0;i<=n+1;i++)
			f[i]=i;
		for(int i=1;i<=n;i++)
		{
			scanf("%lf%lf",&x[i],&y[i]);
			if(x[i]-2*r<eps) merge(0,i);//距离小于直径，合并该点和左边界 
			if(100.00-x[i]-2*r<eps) merge(i,n+1);//距离小于直径，合并该点和右边界 
		}
		for(int i=1;i<=n;i++)
			for(int j=i+1;j<=n;j++)
				if(dis(i,j)-2*r<eps) 
					merge(i,j);距离小于直径，合并两点 
		if(find(0)==find(n+1)) printf("No\n") ;
		else printf ("Yes\n") ; 
	}
    return 0;
}

白嫖一份 ^(o o)^

#include<iostream>
#include<cstdio>
using namespace std;
const int N=1e3+3;
const double eps=1e-7;

double x[N],y[N];
int b[N],q[N],hd,tl;

double sqr(double x){return x*x;}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--){
        int n;
        double r;
        bool f=1;
        scanf("%d%lf",&n,&r);
        for(int i=0;i<n;i++)
            scanf("%lf%lf",&x[i],&y[i]),
            b[i]=1;
        hd=tl=0;
        for(int i=0;i<n;i++)
            if(x[i]<2*r+eps){
                b[q[tl++]=i]=0;
                if(100-x[i]<2*r+eps)f=0;
            }
        while(f&&hd<tl){
            int v=q[hd++];
            for(int i=0;i<n;i++)
                if(b[i]&&sqr(x[i]-x[v])+sqr(y[i]-y[v])<4*r*r+eps){
                    b[q[tl++]=i]=0;
                    if(100-x[i]<2*r+eps)f=0;
                }
        }
        printf(f?"Yes\n":"No\n");
    }
    return 0;
}

---

B. Icebound and Sequence

Icebound and Sequence

快速乘&快速幂&逆元

了解逆元

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

ll mod,q,n;

ll ksc(ll x,ll y)
{
	ll ans=0;
	while(y)
	{
		if(y&1) ans=(ans+x)%mod;
		x=2*x%mod;
		y/=2;
	}	
	return ans;
} 

ll ksm(ll x,ll y)
{
	ll ret=1;
	while(y)
	{
		if(y&1) ret=ksc(ret,x)%mod;
		x=ksc(x,x)%mod;
		y/=2;
	}
	return ret;
} 

int main()
{
	int t;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%lld%lld%lld",&q,&n,&mod);	
		mod=(q-1)*mod;
		ll ans=(ksm(q,n+1)-q+mod)%mod/(q-1);
		printf("%lld\n",ans);		
	}	
}

---

C. 分治

分治

区间DP，逆向思考。

#include<bits/stdc++.h>
using namespace std;
int dp[110][110],cost[110],T,n;

int main()
{
	cin>>T;
	while(T--)
	{
		cin>>n;
		memset(dp,0x3f,sizeof(dp));
		for(int i=1;i<=n;i++)
			scanf("%d",&cost[i]),dp[i][i]=0;
		for(int len=2;len<=n;len++)
		{
			for(int l=1;l+len-1<=n;l++)
			{
				int r=l+len-1;
				for(int k=l;k<=r;k++)
				{
					if(k==l) dp[l][r]=min(dp[l][r],dp[k+1][r]+cost[l]*(r-l));
					else if(k==r) dp[l][r]=min(dp[l][r],dp[l][k-1]+cost[r]*(r-l));
					else dp[l][r]=min(dp[l][r],dp[l][k-1]+dp[k+1][r]+cost[k]*(r-l));
				}
			}
		}
		cout<<dp[1][n]<<endl;
	}
    return 0;
}

---

E. Paper Plane Fly Away

Paper Plane Fly Away

树状数组求逆序对。博客

#include<bits/stdc++.h>
using namespace std;
const int N=200010;
int n,pos[N],love[N],seat[N],c[N],ans[N];

int lowbit(int x)
{
	return x&(-x);
}

void add(int p)
{
	for(int i=p;i<=n;i+=lowbit(i))
		c[i]++;
}

int get(int p)
{
	int sum=0;
	for(int i=p;i;i-=lowbit(i))
		sum+=c[i];
	return sum;
}

int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		int x;
		scanf("%d%d",&x,&love[i]);
		pos[x]=i;
	}
	for(int i=1;i<=n;i++)
		seat[i]=pos[love[i]];
	for(int i=n;i;i--)
	{
		ans[i]+=get(seat[i]);
		add(seat[i]);
	}
	memset(c,0,sizeof(c));
	for(int i=1;i<=n;i++)
	{
		ans[i]+=get(n)-get(seat[i]);
		add(seat[i]);
		printf("%d\n",ans[i]);
	}
	return 0;
}

---

F. Take Apples

Take Apples

首先，“一眼能看出”

if(m<=n) cout<<"Alice"<<endl;
if(m>n&&(m-n)%(s+1)==0) cout<<"Alice"<<endl;

考虑当m>n&&(m-n)%(s+1)!=0时，Alice是否有办法用一步变成后手必胜局面。可以斩三堆使得n<s，这样只用考虑m堆的情况（因为从n取不会改变最终到m堆的先后手情况，如果Bob进一步斩三堆也是稳输），假设只剩下m堆了，Alice后手，这在m%(s+1)==0的情况下Alice稳赢。所以Alice可以在第一步的时候就通过调整，使得调整后呈现后手必胜局面。

Bob必胜态大小为$n$的堆Alice取一个Bob取另一个，大小为$m$的堆Bob取$s+1-$Alice取的数量。
否则Alice能一步变成后手必胜态（$n>s$时同时取3个堆，否则取$m$堆）。

然后这个题目有点坑，输入顺序应该是m,n,s…

#include<bits/stdc++.h>
using namespace std;

int main()
{
    int s,m,n;
    while(~scanf("%d%d%d",&m,&n,&s))
    {
        if(n<=s&&m%(s+1)==0) cout<<"Bob"<<endl;
        else cout<<"Alice"<<endl;
    }
}

---

J. 舔狗

舔狗

拓扑排序，从入度小的地方开始（应该算是一种贪心？），一对点一对点地拆。

#include<bits/stdc++.h>
using namespace std;

const int N=1e6+10;
int n,a[N],ind[N],vis[N];

struct node
{
	int x,p;
};
bool operator < (node a,node b){return a.p>b.p;}
priority_queue<node> q;

int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)	
		scanf("%d",&a[i]),ind[a[i]]++;
	for(int i=1;i<=n;i++)	
		q.push({i,ind[i]});
	while(q.size())
	{
		node u=q.top();	
		q.pop();
		if(vis[u.x]||vis[a[u.x]]) continue;
		n-=2;	
		vis[u.x]=vis[a[u.x]]=1;
		int now=a[a[u.x]];	
		ind[now]--;
		q.push({now,ind[now]});
	}
	cout<<n;
	return 0;
}

---

L. smart robot

smart robot

从0开始往上暴力搜索，本以为会TLE却没有，或许答案必然在一个较小的范围内是可证的？

#include<bits/stdc++.h>
using namespace std;
int mp[100][100],n,num[100];
int dx[4]={0,-1,0,1};
int dy[4]={1,0,-1,0};

struct node
{
	int x,y,step;
};

bool bfs(int x)
{
	int len=0;
	if(!x) num[++len]=1;
	while(x)
	{
		num[++len]=x%10;
		x/=10;
	}
	queue<node>q;
	for(int i=1;i<=n;i++)
		for (int j=1;j<=n;j++)
			if(mp[i][j]==num[1])
				q.push(node{i,j,1});
	while(!q.empty())
	{
		node temp=q.front();
		q.pop();
		if(temp.step==len) return true;
		for(int i=0;i<4;i++)
		{
			int xx=temp.x+dx[i];
			int yy=temp.y+dy[i];
			if(xx<1||xx>n||yy<1||yy>n) continue;
			if(mp[xx][yy]!=num[temp.step+1]) continue;
			q.push(node{xx,yy,temp.step+1});	
		}
	}
	return false;			
}
int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			cin>>mp[i][j];
	for(int i=0;;i++)
	{
		if(!bfs(i))
		{
			cout<<i<<endl;
			return 0;
		} 
	}
}