【OEIS】Xor Sum AtCoder2272

分析：

OEIS大法好哇。。。

暴力打表然后去OEIS上面搜一下，

get…
但似乎没有给公式啊
嗯。。再看看

似乎是这玩意的前缀和啊，这玩意给了个递推公式。。。
然后用数学老师讲的方法：
$a_1=a_1$
$a_2=a_1$
$a_3=a_1+a_2$
$a_4=\ \ \ \ \ \ \ \ \ a_2$
$a_5=\ \ \ \ \ \ \ \ \ a_2+a_3$
$a_6=\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ a_3$
$a_7=\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ a_3+a_4$
然后就能很轻松地发现：
$S_{2*n}=3\times S_n-a_n$
$S_{2*n+1}=3\times S_n+a_{n+1}$
然后就可以求啦。。。先预处理出$a_1到 a_{1000000}$，然后就可以O(logn)求了。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#define SF scanf
#define PF printf
#define MAXN 1000010
#define INF 0x3FFFFFFF
#define MOD 1000000007
using namespace std;
typedef long long ll;
ll n,n1;
ll a[MAXN];
map<ll,ll> mp;
ll p(ll x){
    if(mp.count(x))
        return mp[x];
    if(x<=1000000)
        return a[x];
    if(x==1)
        return 1;
    if(x%2ll==1ll)
        mp[x]=(p(x/2ll)+p(x/2ll+1ll))%MOD;
    else
        mp[x]=p(x/2ll);
    return mp[x];
}
ll que(ll x){
    if(x==1)
        return 1;
    if(x%2ll==1ll)
        return (que(x/2ll)*3ll+p(x/2ll+1ll))%MOD;
    else
        return (que(x/2ll)*3ll-p(x/2ll)+MOD)%MOD;
}
void prepare(){
    a[1]=1;
    for(int i=1;i<=1000000;i++){
        if(i%2==1)
            a[i]=(a[i/2]+a[i/2+1])%MOD; 
        else
            a[i]=a[i/2];
    }
}   
int main(){
    prepare();
    SF("%lld",&n1);
    n1++;
    PF("%lld\n",que(n1));
}