結果

問題 No.685 Logical Operations
ユーザー maroon_kuri
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#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__<<" -- ")
#else
class 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 CAPITAL
void Yes(bool ex=true){
#ifdef CAPITAL
cout<<"YES"<<endl;
#else
cout<<"Yes"<<endl;
#endif
if(ex)exit(0);
}
void No(bool ex=true){
#ifdef CAPITAL
cout<<"NO"<<endl;
#else
cout<<"No"<<endl;
#endif
if(ex)exit(0);
}
const ll infLL=LLONG_MAX/3;
#ifdef int
const int inf=infLL;
#else
const int inf=INT_MAX/2-100;
#endif
const 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);
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0