結果
問題 | No.685 Logical Operations |
ユーザー |
![]() |
提出日時 | 2018-05-11 23:19:43 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 3,955 bytes |
コンパイル時間 | 1,683 ms |
コンパイル使用メモリ | 177,140 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-28 09:02:26 |
合計ジャッジ時間 | 2,498 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 27 |
ソースコード
#include <bits/stdc++.h>using namespace std;using ll=int64_t;#define int ll#define FOR(i,a,b) for(int i=int(a);i<int(b);i++)#define REP(i,b) FOR(i,0,b)#define MP make_pair#define PB push_back#define ALL(x) x.begin(),x.end()#ifdef MAROON_LOCAL#define cerr (cerr<<"-- line "<<__LINE__<<" -- ")#elseclass CerrDummy{}cerrDummy;template<class T>CerrDummy& operator<<(CerrDummy&cd,const T&){return cd;}using charTDummy=char;using traitsDummy=char_traits<charTDummy>;CerrDummy& operator<<(CerrDummy&cd,basic_ostream<charTDummy,traitsDummy>&(basic_ostream<charTDummy,traitsDummy>&)){return cd;}#define cerr cerrDummy#endif#define REACH cerr<<"reached"<<endl#define DMP(x) cerr<<#x<<":"<<x<<endl#define ZERO(x) memset(x,0,sizeof(x))#define ONE(x) memset(x,-1,sizeof(x))using pi=pair<int,int>;using vi=vector<int>;using ld=long double;template<class T,class U>ostream& operator<<(ostream& os,const pair<T,U>& p){os<<"("<<p.first<<","<<p.second<<")";return os;}template<class T>ostream& operator <<(ostream& os,const vector<T>& v){os<<"{";REP(i,(int)v.size()){if(i)os<<",";os<<v[i];}os<<"}";return os;}ll read(){ll i;scanf("%" SCNd64,&i);return i;}void printSpace(){printf(" ");}void printEoln(){printf("\n");}void print(ll x,int suc=1){printf("%" PRId64,x);if(suc==1)printEoln();if(suc==2)printSpace();}string readString(){static char buf[3341000];scanf("%s",buf);return string(buf);}char* readCharArray(){static char buf[3341000];static int bufUsed=0;char* ret=buf+bufUsed;scanf("%s",ret);bufUsed+=strlen(ret)+1;return ret;}template<class T,class U>void chmax(T& a,U b){if(a<b)a=b;}template<class T,class U>void chmin(T& a,U b){if(b<a)a=b;}template<class T>T Sq(const T& t){return t*t;}#define CAPITALvoid Yes(bool ex=true){#ifdef CAPITALcout<<"YES"<<endl;#elsecout<<"Yes"<<endl;#endifif(ex)exit(0);}void No(bool ex=true){#ifdef CAPITALcout<<"NO"<<endl;#elsecout<<"No"<<endl;#endifif(ex)exit(0);}const ll infLL=LLONG_MAX/3;#ifdef intconst int inf=infLL;#elseconst int inf=INT_MAX/2-100;#endifconst int mod=1000000007;template<class T,class U>void add(T& a,U b){a=((ll)a+b)%mod;}template<class T,class U>void sub(T& a,U b){a=((ll)a-b%mod+mod)%mod;}template<class T,class U>void mult(T& a,U b){a=((ll)a*b)%mod;}ll modPow(ll a,ll p){ll s=1;while(p){if(p&1)mult(s,a);mult(a,a);p>>=1;}return s;}ll modInv(ll a){return modPow(a,mod-2);}const int Vmax=1000010;int fact[Vmax],factInv[Vmax];void InitFact(){fact[0]=1;FOR(i,1,Vmax){fact[i]=fact[i-1];mult(fact[i],i);}factInv[Vmax-1]=modInv(fact[Vmax-1]);for(int i=Vmax-2;i>=0;i--){factInv[i]=factInv[i+1];mult(factInv[i],i+1);}}int Choose(int n,int k){return ll(fact[n])*factInv[n-k]%mod*factInv[k]%mod;}int Binom(int a,int b){return ll(fact[a+b])*factInv[a]%mod*factInv[b]%mod;}map<pi,int> memo1,memo2,memo3;/*int calc1(int a,int b){return 0;}**/int calc2(int a,int b){if(a==0&&b==0)return 0;if(memo2.count(pi(a,b)))return memo2[pi(a,b)];int res=0;cerr<<a<<" "<<b<<" "<<res<<endl;add(res,calc2(a/2,b/2));cerr<<a<<" "<<b<<" "<<res<<endl;if(a>0){add(res,1);add(res,calc2((a-1)/2,b/2));}cerr<<a<<" "<<b<<" "<<res<<endl;if(b>0){add(res,1);add(res,calc2(a/2,(b-1)/2));}cerr<<a<<" "<<b<<" "<<res<<endl;if(a>0&&b>0)add(res,calc2((a-1)/2,(b-1)/2));cerr<<a<<" "<<b<<" "<<res<<endl;return memo2[pi(a,b)]=res;}int calc3(int a,int b){if(a==0&&b==0)return 0;if(memo3.count(pi(a,b)))return memo3[pi(a,b)];int res=0;add(res,calc3(a/2,b/2));if(a>0)add(res,calc3((a-1)/2,b/2));if(b>0)add(res,calc3(a/2,(b-1)/2));if(a>0&&b>0)add(res,calc2((a-1)/2,(b-1)/2));cerr<<a<<" "<<b<<" "<<res<<endl;return memo3[pi(a,b)]=res;}signed main(){int n=read();int ans=calc3(n,n);mult(ans,modInv(2));print(ans);}