AOJ7Redraiment猜想

redraiment在家极度无聊，于是找了张纸开始统计素数的个数。
设函数f(n)返回从1->n之间素数的个数。
redraiment发现:

f(1) = 0
f(10) = 4
f(100) = 25
…

满足g(m) = 17 * m2 / 3 - 22 * m / 3 + 5 / 3
其中m为n的位数。
他很激动，是不是自己发现了素数分布的规律了！
请你设计一个程序，求出1->n范围内素数的个数，来验证redraiment是不是正确的，也许还可以得诺贝尔奖呢。^_^

Input
输入包括多组数据。
每组数据仅有一个整数n (1≤n≤10,000,000)。
输入以０结束

Output
对于每组数据输入，输出一行，为1->n(包括n)之间的素数的个数

Sample Input
Original Transformed

1
10
65
100
0
1[EOL] 10[EOL] 65[EOL] 100[EOL] 0[EOL] [EOF]

Sample Output
Original Transformed

0
4
18
25

代码：

    #include <stdio.h>
    int a[10000001];
    int main()
    {
    int i, j;
    int n;
    while(scanf("%ld", &n) != EOF&&n)
    {
         int k = 1;
         if(n == 1)
            printf("0\n");
         else if(n == 2)
            printf("1\n");
         else
         {
            int s = 1;
            for(i = 1; i <= n; i++)
                a[i] = i;
            a[1] = 0;
            for(i = 3; i <= n; i += 2)
            {
                for(j = 3; j <= i&&i * j <= n; j += 2)
                    a[i*j] = 1;
                if(a[i] != 1)
                    s++;
            }
            printf("%ld\n", s);
         }
    }
     return 0;
}