ACdream 1112 Alice and Bob (sg函数的变形+素数筛)

题意：有N个数，Alice 和 Bob 轮流对这些数进行操作，若一个数 n=a*b且a>1，b>1，可以将该数变成 a 和 b 两个数；

         或者可以减少为a或b，Alice先，问谁能赢

思路：首先单看对每个数进行除法的操作，我们可以知道其实是在除以每个数的素因子或素因子之间的积

         比如 70=2*5*7 我们可以变成 10（2*5）或 14（2*7) 或 35（5*7）或 2 或 5 或 7 或 1 这七种状态

         当我们把他们（2，5，7）当作3个石子也就是一堆时，然而实际上我们是将这堆石子进行nim游戏

         我拿走一个石子 =》 10（2*5） 我拿走了石子7 

                                   14 （2*7） 我拿走了石子5

                                   35 （5*7） 我拿走了石子2

         我拿走两个石子 =》 2    我拿走了石子5 和 石子7 

                                    5    我拿走了石子2 和 石子7 

                                    7    我拿走了石子2 和 石子5

         我拿走三个石子 =》 1     我拿走了石子2 和 石子5 和 石子7

         接下来我们分析把一个数n=a*b变成 a 和 b ，其实这里上面的思想很像，把它当作石子的分堆

         我可以分成             第一种 10（2*5） 和 7

                                    第二种 14（2*7） 和 5

                                    第三种 35（5*7） 和 2

         综上所诉，根据正整数唯一分解定理，任何一个正整数x必然有x=(p1^r1)*(p2^r2)*......*(pn^rn)

         定义sum=r1+r2+...+rn，这个sum的值就是这堆石子的总数，那么sg=sg[sum1]^sg[susm2]^....

         问题又来了？ 这个sum我们应该如何求呢？

         我们可以通过素数筛得到每一个数的最小质因子，我们得到一个类似于递推的公式

         一个正整数的质因子的个数=（这个正整数 / 这个数的最小质因子 所得数） 的质因子个数 + 1（也就是加上这是最小质因子的数量 1）

         接下来代码实现就可以了

代码：

#include <iostream>
#include <cstring>
#include <cstdio>
#define maxn 5000000
using namespace std;

int prime[maxn+5],k=0,samll[maxn],sum[maxn];
bool visit[maxn+5];
int sg[105];

void get_prime()
{
    memset(visit,false,sizeof(visit));
    memset(samll,0,sizeof(samll));
    memset(sum,0,sizeof(sum));
    for(int i=2;i<=maxn;i++)
    {
        if(visit[i]==false)
        {
            prime[k++]=i;
            for(int j=i+i;j<=maxn;j+=i)
            {
                visit[j]=true;
                if(samll[j]==0) samll[j]=i;
            }
            samll[i]=i;
        }
    }
    for(int i=2;i<=maxn;i++)
    sum[i]=sum[i/samll[i]]+1;
}

int get(int n)
{
    if(sg[n]!=-1) return sg[n];
    bool vis[105];
    memset(vis,false,sizeof(vis));
    for(int i=1;i<=n;i++)
    vis[get(n-i)]=true;
    for(int i=1;i<=n/2;i++)
      vis[get(i)^get(n-i)]=true;
    int k;
    for(int i=0;i<105;i++)
    {
        if(vis[i]==false)
        {
           return sg[n]=i;
        }
    }
}

int main()
{
   get_prime();
   memset(sg,-1,sizeof(sg));
   sg[0]=0;
   sg[1]=1;
   for(int i=2;i<=100;i++)
   {
       get(i);
   }
   int n;
   while(cin>>n)
   {
       int ans=0;
       for(int i=1;i<=n;i++)
       {
           int x;
           scanf("%d",&x);
           ans=ans^sg[sum[x]];
       }
       if(ans)
       cout<<"Alice"<<endl;
       else
       cout<<"Bob"<<endl;
   }
   return 0;
}