HDU 4342 History repeat itself 模拟

来源：http://acm.hdu.edu.cn/showproblem.php?pid=4342

题意：首先让求第几个非平方数，然后求从1到该数之间的每个sqrt(i)的下取整的和。

思路：一个简单的模拟题目，但是由于数据范围大，需要用__int64。我们可以首先把平方数筛选出来，假如让求第n个非平方数的话，看n前面有多少个平方数，假设有x个，则第n个非平方数就是n+x。注意两种特殊情况，即n + x是完全平方数，则加1即可。另一种是n是完全平方数。之后求和的时候，注意到一个平方数到另一个平方数之间的数开根号下取整的结果是一样的。因此不需要一个一个循环。

代码：

#include <iostream>
#include <cstdio>
#include <string.h>
#include <cmath>
using namespace std;

typedef long long LL;
#define CLR(arr,val) memset(arr,val,sizeof(arr))
LL N = 3000000000;
const int M = 1000010;
LL squre[M];
void init(){
	CLR(squre,0);
	int cnt = 0;
	for(long long i = 1;;++i){
	  if(i * i > N)
		  break;
	  //printf("ss");
	  squre[cnt++] = i * i;
	}
	//for(int i = 0;i < 10; ++i)
	//	printf("%d ",squre[i]);
	//printf("\n");
}
int find(int x){
	int ans = 0;
	for(int i = 0;;++i){
		if(squre[i] <= x && squre[i+1] > x){
		  ans = i;
		  break;
		}
	}
	return ans + 1;
}
int min(int a,int b){
	return a<b?a:b;
}
int main(){
	init();
	int numcase;
	scanf("%d",&numcase);
	while(numcase--){
	   __int64 n;
	   scanf("%I64d",&n);
	   int x = find(n);
	   n = n + x;
	   __int64 y = sqrt(n*1.0);
	   __int64 z = n / y;
	   if(y == z)
		   n++;
	   __int64 sum = 1;
	   __int64 p = 2,q = 1;
	   while(p <= n){
		     __int64 s = sqrt(squre[q] + 0.5) - 1;
			 __int64 mmin = min(n,squre[q] - 1);
			 __int64 num = mmin - p + 1 ;
			 sum += (s * num);
			 p = mmin + 1;
			 q++;
	   }
	   printf("%I64d %I64d\n",n,sum);
	}
	return 0;
}