Project Euler 622

題目鏈接：點擊打開鏈接

題意：題目定義完美洗牌方式，s[i]表示i張牌時需要洗多少次回到最初情況，求s[i]=60的所有i的總和。

易知，我們每次洗牌可以轉化爲公式(2*i-1)(mod n-1)(n爲牌數)，也就是s[i]等於 2 mod(i-1)的階數，進而對求s[i]=60的i的總和可以轉化爲求2^60-1的所有因數加一的和(需要特別把2的低次冪-1的因數給排除掉)。

代碼：

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define PI acos(-1.0)
#define mod 1000000007
using namespace std;
typedef long long ll;
map<ll,bool>mp;
ll dfs(ll n,ll pos,ll pp[],ll num[],ll sum){
    if(pos==n){
        if(mp[sum])return 0;
        else {
            mp[sum]=true;
            return sum+1;
        }
    }
    ll ret=1,res=0;
    for(int i=0;i<=num[pos];i++){
        res+=dfs(n,pos+1,pp,num,sum*ret);
        ret*=pp[pos];
    }
    return res;
}
ll update(ll x){
    ll top=0;
    ll pp[100],num[100];
    for(ll i=2;i*i<=x;i++){
        if(x%i==0){
            pp[top++]=i;num[top-1]=0;
            while(x%i==0){
                num[top-1]++;
                x/=i;
            }
        }
    }
    if(x!=1){
        pp[top]=x;
        num[top++]=1;
    }
    return dfs(top,0,pp,num,1);
}
int main(){
    ll n=1;
    for(ll i=1;i<=60;i++){
        n*=2;
        printf("%lld %lld\n",i,update(n-1));
    }
    return 0;
}