結果

問題 No.685 Logical Operations
ユーザー WA_TLEWA_TLE
提出日時 2018-05-11 23:07:38
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 2,823 bytes
コンパイル時間 1,335 ms
コンパイル使用メモリ 126,764 KB
最終ジャッジ日時 2025-01-05 10:31:59
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<deque>
#include<queue>
#include<vector>
#include<algorithm>
#include<iostream>
#include<set>
#include<cmath>
#include<tuple>
#include<string>
#include<chrono>
#include<functional>
#include<iterator>
#include<random>
#include<unordered_set>
#include<array>
#include<map>
#include<iomanip>
#include<assert.h>
using namespace std;
typedef long long int llint;
typedef long double lldo;
#define mp make_pair
#define mt make_tuple
#define pub push_back
#define puf push_front
#define pob pop_back
#define pof pop_front
#define fir first
#define sec second
#define res resize
#define ins insert
#define era erase
//cout<<setprecision(20)
//cin.tie(0);
//ios::sync_with_stdio(false);
const llint mod=1000000007;
const llint big=4.19e18+1;
const long double pai=3.141592653589793238462643383279502884197;
const long double eps=1e-15;
template <class T,class U>void mineq(T& a,U b){if(a>b){a=b;}}
template <class T,class U>void maxeq(T& a,U b){if(a<b){a=b;}}
llint gcd(llint a,llint b){if(a%b==0){return b;}else return gcd(b,a%b);}
llint lcm(llint a,llint b){return a/gcd(a,b)*b;}
template<class T> void SO(T& ve){sort(ve.begin(),ve.end());}
template<class T> void REV(T& ve){reverse(ve.begin(),ve.end());}
int LBI(vector<int>&ar,llint in){return lower_bound(ar.begin(),ar.end(),in)-ar.begin();}
int UBI(vector<int>&ar,llint in){return upper_bound(ar.begin(),ar.end(),in)-ar.begin();}
int main(void){
	llint n,i;cin>>n;
	//0 0-> 0 0 0
	//0 1-> 0 1 1
	//1 1-> 1 0 1
	//めちゃ複雑な桁DP
	//orがでかいのはわかる xor<orで
	//xor==orのとき、x==y==0 かand<xorになる
	//だから別々に計算
	llint dp[8]={0};
	llint ans=0;
	dp[0]=1;
	//x<yかどうか(=になってない),y<nが確定したかどうかand<xorかどうか
	//4,2,1
	for(int h=60;h>=0;h--){
		llint ddp[8]={0};
		if(n&(1LL<<h)){
			for(i=0;i<8;i++){
				ddp[i|2]+=dp[i];//0 0
				ddp[i|5]+=dp[i];//0 1
				if(i&4){ddp[i|3]+=dp[i];}//1 0
				if(i&1){ddp[i]+=dp[i];}//1 1
			}
		}else{
			for(i=0;i<8;i++){
				ddp[i]+=dp[i];//0 0
				if(i&2){ddp[i|5]+=dp[i];}//0 1
				if(i&4){ddp[i|1]+=dp[i];}//1 0
				if((i&1)&&(i&2)){ddp[i]+=dp[i];}//1 1
			}
		}
		for(i=0;i<8;i++){dp[i]=ddp[i]%mod;}
	}
	ans=dp[7]+dp[5]+dp[3]+dp[1];
	//xor==orの判定 これを引く
	for(i=0;i<8;i++){dp[i]=0;}
	dp[0]=1;
	//x<yかどうか(=になってない),y<nが確定したかどうか
	//4,2 
	//1のbitはないです
	for(int h=60;h>=0;h--){
		llint ddp[8]={0};
		if(n&(1LL<<h)){
			for(i=0;i<8;i++){
				ddp[i|2]+=dp[i];//0 0
				ddp[i|4]+=dp[i];//0 1
				if(i&4){ddp[i|2]+=dp[i];}//1 0
			}
		}else{
			for(i=0;i<8;i++){
				ddp[i]+=dp[i];//0 0
				if(i&2){ddp[i|4]+=dp[i];}//0 1
				if(i&4){ddp[i]+=dp[i];}//1 0
			}
		}
		for(i=0;i<8;i++){dp[i]=ddp[i]%mod;}
	}
	ans-=dp[6]+dp[4]+dp[2]+dp[0];
	ans++;ans+=mod;ans%=mod;
	cout<<ans<<endl;
	return 0;
}
0