EOJ Monthly 2019.2 ABDEF题解

A. 回收卫星

A题wa自闭了，求哪位过了A的老哥借鉴下代码.........

update：突然发现二分的时候最大范围应该是2e9，之前是1e9，改了这个就过了....

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll inf=2e9;
ll a[5];
void ask(int tp,int m)
{
	for(int i=1;i<=3;i++)a[i]=0;
	a[tp]=m;
	printf("0 %lld %lld %lld\n",a[1],a[2],a[3]);	
}
ll gao(int tp)
{
	ll l=0,r=inf,m,res=0;
	int op;
	while(l<r)
	{
		m=(l+r)/2;
		ask(tp,m);
		cin>>op;
		if(op)l=m+1;
		else r=m;
	}
	ask(tp,l);cin>>op;
	if(!op)l--;
	res+=l;
	l=-inf,r=0;
	while(l<r)
	{
		m=(l+r)/2;
		ask(tp,m);
		cin>>op;
		if(op)r=m-1;
		else l=m;		
	}
	ask(tp,r);cin>>op;
	if(!op)r++;
	res+=r;
	return res/2;
}
int main()
{
	ll x,y,z;
	x=gao(1);
	y=gao(2);
	z=gao(3);
	printf("1 %lld %lld %lld\n",x,y,z);
}

B. 解题

思路：我们先预处理一下，设ai为第 i 位到第 n 为组成的数mod m的值，那么任意区间 [ i , j ] mod m的值可以表示为ai - aj-1，那问题就转化成找 i ，j，使得ai = aj，且这个区间转化的数最小，那么也就是 i 尽量大，j 尽量小，那么我们倒着遍历找到第一个 i ，j，ai=aj 就是答案。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=1e6+10,N=5e7+10;
char s[maxn];
int a[maxn];
void up(int i,int j,int& l,int& r)
{
	if(!l)l=i,r=j;
}
bitset<N>mp;
int main()
{
	int q,m,n;
	scanf("%s",s+1);n=strlen(s+1);
	scanf("%d",&q);
	while(q--)
	{
		scanf("%d",&m);
		int y=1,l=0,r=n+1;
		mp.reset();
		for(int j,i=n;i;i--)
		{
			a[i]=(a[i+1]+y*(s[i]-'0'))%m;
			y=y*10%m;
			if(mp[a[i]])
			{
				for(j=i+1;j<=n;j++)
				if(a[j]==a[i])break;
				up(i,j-1,l,r);
			}
			else if(a[i]%m==0)up(i,n,l,r);
			mp[a[i]]=1;
			if(l)break;
		}
		printf("%d %d\n",l,r);
	}
}

D. 进制转换

签到题，注意别爆 ll

#include<bits/stdc++.h>
#define ll long long
using namespace std;
int get(ll x,int k)
{
	int res=0;
	while(x)
	{
		res++;
		x/=k;
	}
	return res;
}
int main()
{
	int T,k,m;
	scanf("%d",&T);
	while(T--)
	{
		ll l,r;
		scanf("%lld%lld%d%d",&l,&r,&k,&m);
		int tmp=get(r,k);
		if(tmp<=m)
		{
			puts("0");
			continue;
		}
		ll x=pow(1ll*k,m);
		ll ans=r/x-(l-1)/x;
		if(r/x>=k)ans-=(r/(x*k)-(l-1)/(x*k));
		printf("%lld\n",ans);
	}
}

E. 中位数

二分答案ans，把权值大于等于ans的点设为1，小于ans的点设-1，跑起点为1终点为n的最长路，如果最长路权值大于等于0，那么就合法。

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e6+5,inf=1e9;
int n,m,a[maxn],b[maxn],c[maxn],d[maxn],vis[maxn],flag;
int Next[maxn],Laxt[maxn],To[maxn],cnt;
void add(int u,int v)
{
	Next[++cnt]=Laxt[u];Laxt[u]=cnt;To[cnt]=v;
}
void dfs(int u,int fa)
{
	if(vis[u])return;
	if(u==n)
	{
		d[u]=a[u];
		flag=1;
		return;
	}
	d[u]=-inf;
	vis[u]=1;
	for(int i=Laxt[u];i;i=Next[i])
	{
		int v=To[i];
		if(v==fa)continue;
		dfs(v,u);
		d[u]=max(d[u],d[v]);
	}
	d[u]+=a[u];
}
int ok(int k)
{
	for(int i=1;i<=n;i++)
	if(c[i]>=k)a[i]=1;
	else a[i]=-1;
	memset(vis,0,sizeof(vis));
	dfs(1,0);
	return d[1]>=0;
}
int main()
{
	int u,v;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
	scanf("%d",&a[i]),b[i]=c[i]=a[i];
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&u,&v);
		add(u,v);
	}
	dfs(1,0);
	if(!flag)puts("-1"),exit(0);
	sort(b+1,b+1+n);
	int sz=unique(b+1,b+1+n)-b-1;
	int l=1,r=sz,mid;
	while(l<r)
	{
		mid=(l+r)/2;
		if(ok(b[mid]))
		{
			if(l==mid)break;
			l=mid;	
		}
		else r=mid;
	}
	if(ok(b[l+1])&&l<r)l++;
	cout<<b[l];
}

F. 方差

签到题，推一下公式就好了

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=1e6+10;
ll a[maxn],sum[maxn],sum2[maxn],mn=1ll<<60;
ll cal(ll x)
{
	return x*x;
}
int main()
{
	int n,m;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
	sort(a+1,a+1+n);
	for(int i=1;i<=n;i++)
	{
		sum[i]+=sum[i-1]+a[i];
		sum2[i]+=sum2[i-1]+a[i]*a[i];
	}
	for(int i=m;i<=n;i++)
	{
		ll tmp=(sum2[i]-sum2[i-m])*m-cal(sum[i]-sum[i-m]);
		mn=min(mn,tmp);
	}
	cout<<mn;
}