版权声明:欢迎大佬指正! 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;
}