function IsPrime(n)
if n == 1 then return false
if n == 2 then return true
for each m from 2if m * m > n then return trueif n mod m = 0 then return falsem := m + 1

template<uint n, uint m>struct TEST{const static uint r = TEST<n, nextM>::r; //nextM为下一个M，暂不实现。在这里用递归结构代替了循环
};
template<uint n>struct ISPRIME{const static uint r = TEST<n, 2>::r; //从2开始，依次判断每一个可能的m取值，判断代码暂未实现。
};
template<>struct ISPRIME<1>{ //对于算法不能计算的特殊值1，判断为0const static uint r = 0;
};
template<>struct ISPRIME<2>{ //对于算法不能计算的特殊值2，判断为1const static uint r = 1;
};

template<uint n, uint m>struct TEST{// if (n % m == 0) n = 0;// if (m * m > n) m = 0; else ++m;const static uint r = TEST<n, m>::r; //上面两行代码不能写在此处，仅说明逻辑。实际的语法下文再做介绍
};
template<uint m>struct TEST<0, m>{ //n为0的情况const static uint r = 0; //即在非特化的模板代码中，n可以被m整除，因此n被赋值为0，故为合数
};
template<uint n>struct TEST<n, 0>{ //m为0的情况const static uint r = 1; //即在非特化的模板代码中，m * m > n，因此n不能被任何比它小的数整除，故为素数
};
template<uint n>struct ISPRIME{const static uint r = TEST<n, 2>::r; //从2开始，依次判断每一个可能的m取值，判断代码暂未实现。
};
template<>struct ISPRIME<1>{ //对于算法不能计算的特殊值1，判断为0const static uint r = 0;
};
template<>struct ISPRIME<2>{ //对于算法不能计算的特殊值2，判断为1const static uint r = 1;
};

#include <iostream>
typedef unsigned int uint;template<uint n, uint m>struct NEXTN{const static uint r = ((n % m != 0) * n);
};template<uint n, uint m>struct NEXTM{const static uint r = (m * m <= n ? (m + 1) : 0);
};template<uint n, uint m>struct TEST{const static uint r = TEST<NEXTN<n, m>::r, NEXTM<n, m>::r>::r;
};template<uint m>struct TEST<0, m>{const static uint r = 0;
};template<uint n>struct TEST<n, 0>{const static uint r = 1;
};template<uint n>struct ISPRIME{const static uint r = TEST<n, 2>::r;
};template<>struct ISPRIME<1>{const static uint r = 0;
};template<>struct ISPRIME<2>{const static uint r = 1;
};int main()
{int primes[] = {ISPRIME<1>::r, ISPRIME<2>::r, ISPRIME<3>::r, ISPRIME<4>::r,ISPRIME<5>::r, ISPRIME<6>::r, ISPRIME<7>::r, ISPRIME<8>::r,ISPRIME<9>::r, ISPRIME<10>::r, ISPRIME<11>::r, ISPRIME<12>::r,ISPRIME<13>::r, ISPRIME<14>::r, ISPRIME<15>::r, ISPRIME<16>::r,};for (int i = 0; i < sizeof(primes) / sizeof(primes[0]); ++i)std::cout << i + 1 << (primes[i] ? " YES" : " NO") <<std::endl;return 0;
}