結果

問題 No.685 Logical Operations
ユーザー WA_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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0