题目大意:x xor 2x=3x(与x xor 3x=2x等价)求满足等式且小于n的x的个数,与满足等式小于2n的数的个数。

因为异或是不进位的二进制加法,那么因为结果正好和加法相同,那么说明x在二进制上没有相邻的1。那么简单的数位DP就可以求出满足这个的答案了。

再看subtask2,根据打表找规律可得,这就是斐波那契数列的第n+2项(以首项是0来说)。那么只需要O(23⋅lgn)的矩阵乘法就可以了。

代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define LL long long unsigned
const LL MOD = 1e9+7;
LL dp[100][2], L, R, cnt;
int n, a[100];
LL DP(int i, int j, int f) {
if(!i) return 1;
if(!f && ~dp[i][j]) return dp[i][j];
LL ans = 0;
int ed = f ? a[i] : 1;
for(int k = 0; k <= ed; ++ k) if(!k||!j) ans += DP(i-1, k, f && k == ed);
if(!f) dp[i][j] = ans;
return ans;
}
LL solve(LL s, int len = 0) {
for(; s; s >>= 1) a[++ len] = s & 1;
return DP(len, 0, 1);
}
struct Mat { LL a[3][3]; } A, B;
Mat Mul(Mat A, Mat B) {
Mat C;
for(int i = 0; i < 2; ++ i)
for(int j = 0; j < 2; ++ j)
C.a[i][j] = 0;
for(int i = 0; i < 2; ++ i)
for(int j = 0; j < 2; ++ j)
for(int k = 0; k < 2; ++ k)
C.a[i][j] = (C.a[i][j] + A.a[i][k] * B.a[k][j]) % MOD;
return C;
}
Mat ksm(Mat A, LL k) {
Mat C;
for(int i = 0; i < 2; ++ i)
for(int j = 0; j < 2; ++ j)
C.a[i][j] = (i == j);
for(; k; k >>= 1) {
if(k & 1) C = Mul(C, A);
A = Mul(A, A);
}
return C;
}
int main() {
memset(dp, -1, sizeof dp);
int T; scanf("%d", &T);
while(T --) {
scanf("%llu", &R);
A.a[0][0] = A.a[0][1] = A.a[1][0] = 1;
A.a[1][1] = 0;
B.a[0][1] = 0; B.a[0][0] = 1;
A = ksm(A, R+1); A = Mul(A, B);
printf("%llu\n%llu\n", solve(R)-1, A.a[0][0]);
}
}

最新文章

  1. NET 强签名
  2. 不让Win7休眠的设置
  3. Segmetation fault你来的真不是时候
  4. [CoreOS 转载] CoreOS实践指南(七):Docker容器管理服务
  5. 常见面试题之ListView的复用及如何优化
  6. Asp.net操作Excel----NPOI!!!!1
  7. 2014多校第十场1004 || HDU 4974 A simple water problem
  8. hdu 1848 Fibonacci again and again (初写SG函数,详解)
  9. Linux Kconfig及Makefile学习
  10. Spring中使用log4j随笔
  11. 创建基于maven的项目模版
  12. [转]UIWebView的Javascript运行时对象
  13. bzoj1070————2016——3——14
  14. C# Lambda表达式运用
  15. HTML5 CSS3 诱人的实例: 3D立方体旋转动画
  16. JS数组Sort方法的使用
  17. qml: 模块定义与使用
  18. svn创建分支(branch/tag)出现“path”already exists
  19. Android ListVeiw控件(转载整理)
  20. iOS多线程编程之多线程简单介绍(转载)

热门文章

  1. IIS7 IIS7.5 IIS8.5 HTTP 错误 500.19 – Internal Server Error解决方案小记
  2. 剑指Offer:面试题25——二叉树中和为某一值的路径(java实现)
  3. Jmeter组件1. CSV Data Set Config
  4. js实现基础运动
  5. Python TCP客户端
  6. js中RegExp类型
  7. mongodb账号安全操作
  8. ecshop的特点,持续加新
  9. SystemVerilog的历史
  10. VC++ chap13 文档与串行化