HDU1286+线性筛素数

/*
欧拉函数+求与之互素的数的个数
*/
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>
#include<queue>
#include<vector>
#include<map>
#include<math.h>
typedef long long ll;
//typedef __int64 int64;
const int maxn = 40005;
const int inf = 0x7FFFFFFF;
const double pi = acos(-1.0);
const double eps = 1e-8;
using namespace std;
int prime[ 50001 ],shu[ 50001 ];
int getPrime( int n ){
    for( int i=1;i<n;i+=2 ) shu[ i ] = 1;
    for( int i=0;i<n;i+=2 ) shu[ i ] = 0;
    shu[ 1 ] = 0;
    shu[ 2 ] = 1;
    for( int i=3;i<n;i+=2 ){
        if( shu[i]==1 ){
            int t,delta;
            delta = i*2;
            t = delta+i;
            while( t<n ){
                shu[ t ] = 0;
                t += delta;
            }
        }
    }
    int cnt = 1;
    for( int i=2;i<n;i++ ){
        if( shu[i]==1 ){
            prime[ cnt++ ] = i;
            //printf("%d\n",i);
        }
    }
    return cnt;
}
int main(){
    //freopen( "in2.txt","r",stdin );
    //freopen( "out.txt","w",stdout );
    int ca;
    int cnt = getPrime( 50001 );
    
    scanf("%d",&ca);
    while( ca-- ){
        int n;
        scanf("%d",&n);
        
        if( shu[n]==1 ){
            printf("%d\n",n-1);
            continue;
        }
        
        int ans[ maxn ];
        int tt = 0;
        for( int i=1;i<cnt/*prime[i]*prime[i]<=n*/;i++ ){//在这里    WA了4，5次。。。
            if( n%prime[i]==0 ){
                ans[ tt++ ] = prime[i];
                //printf("%d\n",prime[i]);
            }
        }
        int res = n;
        for( int i=0;i<tt;i++ ){
            res = res*( ans[i]-1 )/ans[i];
        }
        //printf("n:%d ",n);
        printf("%d\n",res);
    }
    return 0;
}

没什么好说的。。。wa了几次。。。因为for循环条件没写好。。

欧拉函数。。。