值域范围内和倍数有关的一种分组方法+分组问题处理出上下界+通过调和级数按顺序枚举因倍数:ARC141D

2 m 2m 2m 个数分成 m m m 组,一种合法方案必然是 m + 1 ∼ 2 M m+1\sim 2M m+12M

所以考虑对每个数进行分组

Trick1 值域范围内和倍数有关的一种分组方法

思考方法:每个数 x = k 2 i x=k2^i x=k2i,放在第 k k k 组里

直观理解:每个奇数不断翻倍在一个组内

通过这样子,就可以正好可以分成 m m m

为什么这样分?

  1. 存在一种方法每一组拿一个数
  2. 每组内最多拿一个

所以已经转化为每组恰好选一个的问题了

Trick 2 分组问题处理出上下界+通过调和级数按顺序枚举因倍数

对于每一组,我们可以求出其上下界

上界我们希望越大越好,下界我们希望越小越好

枚举上界时,我们要满足其因数不被选。所以对应之前的组,选越大越好,正好对应上界。

下界同理。满足其倍数不被选,所以比它大的越小越好。

#include<bits/stdc++.h>
using namespace std;
#define int long long
inline int read(){int x=0,f=1;char ch=getchar(); while(ch<'0'||
ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){
x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}return x*f;}
#define Z(x) (x)*(x)
#define pb push_back
//mt19937 rand(time(0));
//mt19937_64 rand(time(0));
//srand(time(0));
#define N 600010
//#define M
//#define mo
int n, m, i, j, k, T;
int UP[N], DW[N], c[N], a[N]; 
vector<int>v[N]; 

void all_no() {
	for(i=1; i<=n; ++i) printf("No\n"); 
}

int check(int x) {
	for(int y=x; y<=2*m; y+=x) if(c[y]) return 0; 
	return 1; 
}

signed main()
{
//	freopen("in.txt", "r", stdin);
//	freopen("out.txt", "w", stdout);
//	T=read();
//	while(T--) {
//
//	}
	n=read(); m=read(); 
	for(i=1; i<=n; ++i) {
		a[i]=k=read(); 
		while(k%2==0) k/=2; v[k].pb(a[i]); 
//		printf("%lld -> %lld\n", a[i], k); 
		
	}
	for(i=1; i<=2*m; i+=2) {
		if(!v[i].size()) return all_no(), 0; 
		for(j=v[i].size()-1; j>=0 && c[v[i][j]]; --j); 
		if(j==-1) return all_no(), 0; 
		UP[i]=j; 
		for(k=v[i][j]*2; k<=2*m; k+=v[i][j]) c[k]++; 
//		printf("UP[%lld]=%lld\n", i, UP[i]) ; 
	}
	for(i=1; i<=2*m; ++i) c[i]=0; 
	for(i=2*m-1; i>=1; i-=2) {
		if(!v[i].size()) return all_no(), 0; 
		for(j=0; j<v[i].size(); ++j) if(check(v[i][j])) break; 
		if(j==v[i].size()) return all_no(), 0; 
		DW[i]=j; c[v[i][j]]=1; 
	}
//	for(i=1; i<=2*m; i+=2) printf("%lld %lld\n", v[i][DW[i]], v[i][UP[i]]); 
	for(i=1; i<=2*m; i+=2) if(DW[i]>UP[i]) return all_no(), 0; 
	for(i=1; i<=n; ++i) {
		k=a[i]; 
		while(k%2==0) k/=2; 
		if(a[i]>=v[k][DW[k]] && a[i]<=v[k][UP[k]]) printf("Yes\n"); 
		else printf("No\n"); 
	}
	return 0;
}
  • 4
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值