2018多校第一场HDU6304

版权声明：欢迎大佬指正！ https://blog.csdn.net/sinat_36215255/article/details/81222821

题目网址：https://cn.vjudge.net/problem/HDU-6304

当时打多校的时候退了好久好久的规律，好久好久，当时我看出来了数目的规律，队友提醒我说用等差数列求，但是我没有听她讲什么，没有交流，gg了，这样不行。

今天重开多校，我和队友一人写了一遍，结果我的一直wa，队友写的一直TE，真是惊讶了，都写自闭了。

最终，Wa的原因是，求每个等差数列数目t的时候，注意（p  - i）/(2*i）+1,这个 i 不能取模，取模之后就错了，超级错。

TE的原因是队友没有把mod设置为常量，不是很明白这样为什么TE。。

关于这个题目的题解，看的一个大神，这里直接把网址给，讲的很明白http://www.cnblogs.com/UnderSilenceee/p/9361371.html

；

include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <string>
#include <cstdio>

using namespace std;
const long long  mod = 1e9+7;
long long  inv2 ;
long long quik_pow(long long x,long long n)
{
    long long res = 1;
    while(n>0)
    {
        if(n&1) res = res*x%mod;
        x = x*x%mod;
        n>>=1;
    }
    return res;
}

long long calc(long long n)
{
    if(n<=1) return n;
    else
        return calc(n/2) + n;
}
long long n;
void  solve()
{
    long long  l =(n/2-30),r = (n/2+30),mid,p = -1;
    while(l<=r)
    {
        mid = (l+r)/2;
        if(calc(mid) > n) r= mid -1;
        else
        {
            l = mid+1;
            p = mid ;
        }
    }
   // cout<<p<<endl;
    long long rest = (( n - calc(p)) %mod + mod)%mod;
    long long ans = 0;
    long long c = 1;
    long long x,y,a1,t,an;
    for(long long i=1; ; i<<=1,c++)
    {
        if(i>p) break;
        x = i%mod;
        y = (i*2)%mod;
        a1 = i%mod;
        t = ((p - i)/(2*i) + 1 )%mod; //注意此处是i 不能取摸
        an = ( y*(t-1)%mod + i )%mod;
        ans = ( ans%mod + c*(a1+an)%mod * t %mod * inv2 % mod) %mod;
    }
    ans = ( ans % mod + (rest *(p+1))% mod +1 )%mod;
    printf("%lld\n",ans%mod);
    //return ans %mod ;
}
int main()
{
     int t;
     inv2 = quik_pow(2,mod-2);
     //cout<<inv2<<endl;
     scanf("%d",&t);
     while(t--)
     {
         scanf("%lld",&n);
         n--;
         solve();
         //printf("%lld\n",ans%mod);
     }
    return 0;
}