算法竞赛——进阶指南——NOIP2009 acwing200. Hankson的趣味题

最新推荐文章于 2023-04-25 23:38:45 发布

夕林山寸

最新推荐文章于 2023-04-25 23:38:45 发布

阅读量347

点赞数

分类专栏： # 算法竞赛——进阶指南

本文链接：https://blog.csdn.net/bjfu170203101/article/details/104392565

版权

探讨NOIP2009竞赛中的一道题目，涉及到Hankson的趣味问题。文章指出，解决策略是确定数x是b的倍数且d的约数，注意到1e9内约数最多的数有1536个约数。由于存在2000组数据，简单的枚举方法可能导致超时。为优化，建议预处理平方根内的质数，并找出d的所有质因子，通过DFS枚举以降低复杂度。

摘要由CSDN通过智能技术生成

x一定是b的倍数，d的约数。

由于1e9内的所有数中，约数最多的数的约数为1536个。

所以我们可以枚举d的约数，判是否成立。

但2000组数据会T

我们可以从组成约数的本质来枚举，预处理出sqrt（2e9）的质数。

找出所有d的质因子（有可能会有大于sqrt(d)的质因子，但只可能有一个）判断下。

然后dfs枚举就行最坏情况2000*1536

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define ls (o<<1)
#define rs (o<<1|1)
#define pb push_back
//#define a(i,j) a[(i)*(m+2)+(j)]  //m是矩阵的列数
const int M = 1e5+7;
/*
int head[M],cnt;
void init(){cnt=0,memset(head,0,sizeof(head));}
struct EDGE{int to,nxt,val;}ee[M*2];
void add(int x,int y,int z){ee[++cnt].nxt=head[x],ee[cnt].to=y,ee[cnt].val=z,head[x]=cnt;}
*/
int vs[M],prime[M],m;
int pt[M],sz,ans;
ll a,b,c,d;
void pre()
{
	int n=sqrt(2e9);
	for(int i=2;i<=n;i++)
	{
		if(vs[i])continue;
		prime[++m]