还看不见绿洲呢,绿洲的存在是来安慰迷失者的吗。
现在的我就像是沙漠中迷失的人,只能奔着我认为对的方向一直前进,向着我认为存在的绿洲前进,虽然过程是困难的,但是只要结果是好的我就已经心满意足了,天多给我点运气吧!
昨天晚上评奖学金我居然评上了有没有,虽然只有600块,但是也是非常开心有没有!出去吃饭!
今天上午在实验室打印东西,发现打印机没有反应,还以为是电脑的问题,搞了好长时间,结果旁边的学长按了下打印机的开关,打印机正常运转了,原来是没有开机,我真的好蠢啊啊啊。
然后就是这几天做题做的比较少了,因为感冒了,很不舒服,更别说在实验室自闭想题了,昨天晚上实验室比较吵,加上自己也没有什么状态,于是一晚上一道题目也没想出来,看到网上也没有这道题目的题解,后来晚上睡不着的时候想了一会,发现我好像证明出来了,今天写篇博客出来。
题目地址
思路的过程比较曲折,我就直接讲做法了。
思路:我们将数字分成两类素数和合数,分别进行处理,对于合数 a a a来说,我们采取这样的策略,选取其最小的质因子 c c c(一定大于1),然后将这个最小的质因子 c c c当成这个合数与某个数的最大公约数,下一步我们只需要找出与 a c \frac{a}{c} ca 互质的数就可以了,我们可以知道 a c \frac{a}{c} ca一定是大于等于2 的,所以我们利用大于1相邻两个数一定是互质的这个结论(很好理解这点),我们可以选出来 a c \frac{a}{c} ca和 a c \frac{a}{c} ca - 1这两个数,然后c, a c \frac{a}{c} ca, a c \frac{a}{c} ca - 1这三个数一定可以组成三角形,至此合数构建结束,还有一个质数,因为质数与某个数的最大公因子d不是其本身就是1,d如果是1 的话一定组不成三角形,这个很好想,所以我们可以按照最大公因数是其本身来考虑,如果最大公因子是其本身,那么组成三角形的某条边就会出现一个1,那么我们剩下的两条边只能选择长度是一样的,所以我们最小最小可以选择这个质数的平方。
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
const int N = 1e6 + 10;
int ans[N];
bool prime[N];
int main ()
{
for (int i = 2; i <= 1e6; i ++)
{
if (prime[i]) continue;
for (int j = i + i; j <= 1e6; j += i) prime[j] = true;
}
ans[1] = 1;
for (int i = 2; i <= 1e6; i ++)
{
if (prime[i]) continue;
ans[i] = 1;
if ((LL) i * i <= 1e6) ans[i * i] -= 1;
}
for (int i = 1; i <= 1e6; i ++) ans[i] += ans[i - 1];
int T; scanf ("%d", &T);
while (T --)
{
int n; scanf ("%d", &n);
printf ("%d\n", ans[n]);
}
return 0;
}