zoj 2286 Sum of Divisors

姐给我找的数论水题。

 

题目是，给出一个数字，求出在100W内因子和不超过这个数字的数的个数。因子和不包括它本身。

 

很显然的方法，打表，然后判断。

 

打表是用类似筛法打的，初始值都赋为1（都是1的倍数），就是2加到2的倍数上，3加到3的倍数上就这么下去。。。。。因为只要求个数，然后排下序就好了。

 

开始用二分，二分没写对，一直出不来结果，后来把答案全打出来了，300+W个数。。。悲剧，机子卡死机了。

 

大于等于3392928的直接输出100W。今早起来把二分改改，顺利AC。。。

 

#include <cstdio> #include <cstdlib> #include <iostream> #include <string.h> #include <queue> #include <stack> #include <math.h> #include <algorithm> #include <limits.h> #define MAX 1000001 using namespace std; int a[MAX]; void init() { int i,k; for(i=1; i<MAX; i++) a[i] = 1; for(i=2; i<MAX; i++) for(k=2; k*i<MAX; k++) a[k*i] += i; } int binsearch(int x) { int beg = 0,end = MAX; int mid; while( beg <= end ) { mid = (beg+end)/2; if( x >= a[mid] && x < a[mid+1] ) return mid; if( x >= a[mid] ) beg = mid; if( x < a[mid] ) end = mid; } return -1; } int main() { int n,ans,i,k; init(); sort(a+1,a+MAX); while( ~scanf("%d",&n) ) { if( n == 1 ) { printf("78499/n"); continue; } if( n >= 3392928 ) { printf("1000000/n"); continue; } if( n == 3392927 ) { printf("999999/n"); continue; } ans = binsearch(n); printf("%d/n",ans); } return 0; }