当前位置: 代码迷 >> 综合 >> Groundhog and 2-Power Representation
  详细解决方案

Groundhog and 2-Power Representation

热度:80   发布时间:2024-02-08 14:09:08.0

题目描述

在这里插入图片描述

分析

从题目的备注
在这里插入图片描述
中可以看出
这道题的答案要用高精度计算和存储
于是模板先打好

struct BigInt
{ll a[B];BigInt(){}BigInt(ll x){a[0]=x;for(int i=0;i<B-1;i++){a[i+1]=a[i]/M;a[i]%=M;}}
};
BigInt add(BigInt x,BigInt y)
{BigInt ret;for(int i=0;i<B;i++)ret.a[i]=x.a[i]+y.a[i];for(int i=0;i<B-1;i++){ret.a[i+1]+=ret.a[i]/M;ret.a[i]%=M;}return ret;
}
BigInt mul(BigInt x,BigInt y)
{BigInt ret;memset(ret.a,0,sizeof(ret.a));for(int i=0;i<B;i++)for(int j=0;i+j<B;j++)ret.a[i+j]+=x.a[i]*y.a[j];for(int i=0;i<B-1;i++){ret.a[i+1]+=ret.a[i]/M;ret.a[i]%=M;}return ret;
}
void print(BigInt x)
{int f=-1;for(int i=B-1;~i;i--)if(x.a[i]){f=i;break;}if(f>=1){if(x.a[f]>0){for(int i=f;i;i--)x.a[i]--,x.a[i-1]+=M;}else{for(int i=f;i;i--)x.a[i]++,x.a[i-1]-=M;}for(int i=0;i<B-1;i++)x.a[i+1]+=x.a[i]/M,x.a[i]%=M;if(!x.a[f]) f--;printf("%lld",x.a[f]);for(int i=f-1;~i;i--)printf("%04lld",abs(x.a[i]));}else printf("%lld",x.a[0]);
}

然后因为题目中有二的幂,所以要再打一个高精度版本的快速幂

BigInt ksm_last_order(int b){BigInt a=BigInt(2);BigInt ret=BigInt(1);while(b){if(b&1) ret=mul(ret,a);a=mul(a,a),b>>=1;}return ret;
}

接下来分析一下给定的输入
发现符号^
被括号给代替了
所以要打括号匹配,并进行递归

int sol(int l,int r)
{int p=0;for(int i=l;i<=r;i++){if(s[i]=='(') p+=ksm(sol(i+1,mat[i]-1)),i=mat[i]+1;if(s[i]=='2') p+=2;}return p;
}

括号匹配

    for(int i=1;i<=len;i++){if(s[i]=='(') st.push(i);if(s[i]==')'){int tp=st.top();st.pop();mat[tp]=i;if(tp==pre)ans=add(ans,ksm_last_order(sol(tp+1,i-1))),pre=i+2,ed=i;//因为调试时最后的2^1和2^0没加,所以用ed标记}}

然后输出就AC了

#include <stack>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define ll long long
using namespace std;
const int B=200,M=10000;
struct BigInt
{ll a[B];BigInt(){}BigInt(ll x){a[0]=x;for(int i=0;i<B-1;i++){a[i+1]=a[i]/M;a[i]%=M;}}
};
BigInt add(BigInt x,BigInt y)
{BigInt ret;for(int i=0;i<B;i++)ret.a[i]=x.a[i]+y.a[i];for(int i=0;i<B-1;i++){ret.a[i+1]+=ret.a[i]/M;ret.a[i]%=M;}return ret;
}
BigInt mul(BigInt x,BigInt y)
{BigInt ret;memset(ret.a,0,sizeof(ret.a));for(int i=0;i<B;i++)for(int j=0;i+j<B;j++)ret.a[i+j]+=x.a[i]*y.a[j];for(int i=0;i<B-1;i++){ret.a[i+1]+=ret.a[i]/M;ret.a[i]%=M;}return ret;
}
void print(BigInt x)
{int f=-1;for(int i=B-1;~i;i--)if(x.a[i]){f=i;break;}if(f>=1){if(x.a[f]>0){for(int i=f;i;i--)x.a[i]--,x.a[i-1]+=M;}else{for(int i=f;i;i--)x.a[i]++,x.a[i-1]-=M;}for(int i=0;i<B-1;i++)x.a[i+1]+=x.a[i]/M,x.a[i]%=M;if(!x.a[f]) f--;printf("%lld",x.a[f]);for(int i=f-1;~i;i--)printf("%04lld",abs(x.a[i]));}else printf("%lld",x.a[0]);
}
const int N=20010;
char s[N],ss[N];
int mat[N];
stack<int> st;
BigInt ans=BigInt(0);
BigInt ksm_last_order(int b){BigInt a=BigInt(2);BigInt ret=BigInt(1);while(b){if(b&1) ret=mul(ret,a);a=mul(a,a),b>>=1;}return ret;
}
int ksm(int b){if(!b) return 1;int a=2,ret=1;while(b){if(b&1) ret*=a;a*=a,b>>=1;}return ret;
}
int sol(int l,int r){int p=0;for(int i=l;i<=r;i++){if(s[i]=='(') p+=ksm(sol(i+1,mat[i]-1)),i=mat[i]+1;if(s[i]=='2') p+=2;}return p;
}
int main(){scanf("%s",ss+1);int lens=strlen(ss+1),len=0;for(int i=1;i<=lens;i++) if(ss[i]=='2'&&ss[i+1]=='(') ss[i]='?';for(int i=1;i<=lens;i++) if(ss[i]!='?') s[++len]=ss[i];int p=0,pre=1,ed;for(int i=1,val;i<=len;i++){if(s[i]=='(') st.push(i);if(s[i]==')'){int tp=st.top();st.pop();mat[tp]=i;if(tp==pre)ans=add(ans,ksm_last_order(val=sol(tp+1,i-1))),pre=i+2,ed=i;}}ans=add(ans,sol(ed+1,len));print(ans),putchar(10);
}
  相关解决方案