HDU 6145 Arithmetic of Bomb II（数学+矩阵快速幂）

Description

众所周知，度度熊非常喜欢数字。

它最近在学习小学算术，沉迷于计算A+B中不能自拔。

但是这回的算术题可并不简单，由于含有表示bomb的#号，度度熊称之为 Arithmetic of Bomb。

Arithmetic of Bomb的目的与普通算术一样，就是计算一些Bomb Expression的结果。比如，”1-2+3”的结果为2。然而，bomb，也就是#号，会展开一些普通表达式，这会导致需要计算的式子超出了度度熊所能理解的范畴。比如”(1-2+3)#(3)”表示”1-2+3”出现了3次，将会被展开为”1-2+31-2+31-2+3”。

为了方便理解，下面给出了Bomb Expression的BNF表示。

<bomb expression> := <bomb term> | <bomb expression> <bomb term>
<bomb term> := <bomb statement> | '(' <bomb statement> ')' '#' '(' <number> ')'
<bomb statement> := <bomb element> | <bomb statement> <bomb element>
<bomb element> := <digit> | '+' | '-' | '*'
<normal expression> := <norm term> | <normal expression> '+' <norm term> | <normal expression> '-' <norm term>
<norm term> := <number> | <norm term> '*' <number>
<number> := <digit> | <non-zero-digit> <number>
<digit> := '0' | <non-zero-digit>
<non-zero-digit> := '1' | '2' | '3' | '4' | '5' | '6' | '7' | '8' | '9'

请先将Bomb Expression中所有的#号展开，使其成为Normal Expression（题目的输入保证展开后是一个合法的Normal Expression），再来计算这个表达式的结果。

Input

第一行为T，表示输入数据组数。

每组数据包含一个Bomb Expression。

●1≤T≤50

●1≤length(Bomb Statement)≤10

●1≤length(Number in Bomb term)≤10

●1≤length(Bomb Expression)≤300 000

Output

对每组数据输出表达式的结果，结果对 1 000 000 007 取模。

Sample Input

6
1-2+3
(1-2+3)#(3)
(1)#(3)
(1+)#(2)1
(2*3+1)#(2)
(2)#(2)1+1(2)#(2)

Sample Output

2
60
111
3
43
343

Solution

将表达式展开不现实，考虑用矩阵简化运算，由于加法和减法类似，所以下面只讨论数字、加法和乘法
记

$$M_{a,b,c,d,e}^i=a*b+c+d*e= \left( \begin{array}{ccc} 1&0&0&0 \\ 10^i*b&0&0& 0\\ a*b&0&0&0\\ c&d*e&d&1 \end{array} \right)$$
1.在 $M_{a,b,c,d,e}^i$两端加位

若要通过左乘一个矩阵实现在$a*b+c+d*e$左端的$a$加上一个最高位$f$，即要把$M_{a,b,c,d,e}^i$变成$M_{fa,b,c,d,e}^{i+1}$，只需把矩阵中的$(10^i*b,a*b)^T$变成$(10^{i+1}*b,fa*b)^T$，由于$fa*b=a*b+f*10^i*b$，故有

$$\left( \begin{matrix} 10^{i+1}*b\\ fa*b \end{matrix}\right)= \left( \begin{matrix} 10&0\\ f&1\\ \end{matrix} \right) \left( \begin{matrix} 10^i*b\\ a*b\\ \end{matrix} \right)$$
若要通过右乘一个矩阵实现在 $a*b+c+d*e$右端的 $e$加上一个最低位$f$，即要把 $M_{a,b,c,d,e}^i$变成 $M_{a,b,c,d,ef}^{i}$，只需把矩阵中的 $(d*e,d)$变成 $(d*ef,d)$，由于 $d*ef=10*d*e+d*f$，故有
$$\left( \begin{array}{ccc} d*ef&d\\ \end{array} \right)= \left( \begin{array}{ccc} d*e&d\\ \end{array} \right) \left( \begin{array}{ccc} 10&0\\ f&1\\ \end{array} \right)$$
故对于每个一位数 $f$，令
$$M_{f}= \left( \begin{array}{ccc} 1&0&0&0 \\ 0&10&0& 0\\ 0&f&1&0\\ 0&0&0&1 \end{array} \right)$$
那么通过对 $M_{a,b,c,d,e}^i$左乘右乘 $M_f$即可实现在左边右边加一位 $f$

2.在$M_{a,b,c,d,e}^i$两端加数字

若要通过左乘一个矩阵实现在$a*b+c+d*e$左端加一个数字$f$，即要把$M_{a,b,c,d,e}^{i}$变成$M_{f,1,a*b+c,d,e}^{1}$，只需把矩阵中的$(1,10^i*b,a*b,c)^T$变成$(1,10*1,f*1,a*b+c)^T$，故有

$$\left( \begin{array}{ccc} 1\\ 10*1\\ f*1\\ a*b+c \end{array} \right)= \left( \begin{array}{ccc} 1&0&0&0\\ 10&0&0&0\\ f&0&0&0\\ 0&0&1&1 \end{array} \right) \left( \begin{array}{ccc} 1\\ 10^i*b\\ a*b\\ c \end{array} \right)\ = M_f \left( \begin{array}{ccc} 1&0&0&0\\ 1&0&0&0\\ 0&0&0&0\\ 0&0&1&1 \end{array} \right) \left( \begin{array}{ccc} 1\\ 10^i*b\\ a*b\\ c \end{array} \right)$$
若要通过右乘一个矩阵实现在 $a*b+c+d*e$右端加一个数字 $f$，即要把$M_{a,b,c,d,e}^{i}$变成 $M_{a,b,c+d*e,1,f}^{i}$，只需把矩阵中的 $(c,d*e,d,1)$变成 $(c+d*e,f,1,1)$，故有
$$\left( \begin{array}{ccc} c+d*e&f&1&1 \end{array} \right)= \left( \begin{array}{ccc} c&d*e&d&1 \end{array} \right) \left( \begin{array}{ccc} 1&0&0&0\\ 1&0&0&0\\ 0&0&0&0\\ 0&f&1&1 \end{array} \right)= \left( \begin{array}{ccc} c&d*e&d&1 \end{array} \right) \left( \begin{array}{ccc} 1&0&0&0\\ 1&0&0&0\\ 0&0&0&0\\ 0&0&1&1 \end{array} \right) M_f$$
故对于加号，令
$$M_{+}= \left( \begin{array}{ccc} 1&0&0&0 \\ 1&0&0& 0\\ 0&0&0&0\\ 0&0&1&1 \end{array} \right)$$
那么通过对 $M_{a,b,c,d,e}^i$左乘右乘 $M_+$即可实现在左边右边带上加号

注：类似的可以定义减号

$$M_{-}= \left( \begin{array}{ccc} 1&0&0&0 \\ -1&0&0& 0\\ 0&0&0&0\\ 0&0&1&1 \end{array} \right)$$

3.在$M_{a,b,c,d,e}^i$两端乘数字

若要通过左乘一个矩阵实现在$a*b+c+d*e$左端乘一个数字$f$，即要把$M_{a,b,c,d,e}^{i}$变成$M_{f,a*b,c,d,e}^{1}$，只需把矩阵中的$(10^i*b,a*b)^T$变成$(10*a*b,f*a*b)^T$，故有

$$\left( \begin{array}{ccc} 10*a*b\\ f*a*b\\ \end{array} \right)= \left( \begin{array}{ccc} 0&10\\ 0&f \end{array} \right) \left( \begin{array}{ccc} 10^i*b\\ a*b\\ \end{array} \right)= \left( \begin{array}{ccc} 10&0\\ f&1 \end{array} \right) \left( \begin{array}{ccc} 0&1\\ 0&0 \end{array} \right) \left( \begin{array}{ccc} 10^i*b\\ a*b\\ \end{array} \right)$$
若要通过右乘一个矩阵实现在 $a*b+c+d*e$右端乘一个数字 $f$，即要把$M_{a,b,c,d,e}^{i}$变成 $M_{a,b,c,d*e,f}^{i}$，只需把矩阵中的 $(d*e,d)$变成 $(d*e*f,d*e)$，故有
$$\left( \begin{array}{ccc} d*e*f&d*e \end{array} \right)= \left( \begin{array}{ccc} d*e&d \end{array} \right) \left( \begin{array}{ccc} f&1\\ 0&0 \end{array} \right)= \left( \begin{array}{ccc} d*e&d \end{array} \right) \left( \begin{array}{ccc} 0&1\\ 0&0 \end{array} \right) \left( \begin{array}{ccc} 10&0\\ f&1 \end{array} \right)$$
故对于乘号，令
$$M_{*}= \left( \begin{array}{ccc} 1&0&0&0 \\ 0&0&1& 0\\ 0&0&0&0\\ 0&0&0&1 \end{array} \right)$$
那么通过对 $M_{a,b,c,d,e}^i$左乘右乘 $M_*$即可实现在左边右边带上乘号

综上，通过把每一个位$f,+,-,*$分别变成$M_f,M_{+},M_{-},M_{*}$后，原表达式即转化为一些矩阵的乘积，而对于重复的表达式，只需求矩阵幂即可

注：

1.$M_{-}$的意义：当在$a*b+c+d*e$右方加上$-f$时，我们先把$f$用和加法相同的方法放在了第四行第二列，然后用一个符号变量记录这个负号，当后面又出现加法或者减法时，这个$-f$会被并到第四行第四列中的$c+d*e$中，此时由于通过符号变量的记录我们用$c+d*e$减去$f$而不是加上$f$来得到$c+d*e-f$的正确值

2.当处理重复表达式时，我们会用一个矩阵来表示这个重复表达式，然后用矩阵幂去快速计算，但是由于在计算过程中需要用一个符号变量记录$M_{a,b,c,d,e}^{i}$中的$d*e$的真实符号，故对于第一个该表达式，其符号取决于之前表达式继承下来的符号，而后面的表达式的符号是取决于这个表达式中的最后一个运算是加还是减，故第一个表达式要单独计算，后面的重复部分再用幂运算加速

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define maxn 300005
#define mod 1000000007
vector<string>term;
vector<ll>num;
char c[maxn]; 
int T;
void Deal()
{
    term.clear(),num.clear();
    int n=strlen(c);
    for(int i=0;i<n;)
    {
        string s;
        if(c[i]=='(')
        {
            i++;
            while(c[i]!=')')s.push_back(c[i++]);
            i+=3;
            ll t=0;
            while(c[i]!=')')t=t*10+c[i++]-'0';
            if(t)term.push_back(s),num.push_back(t);
            i++;
        }
        else
        {
            while(i<n&&c[i]!='(')s.push_back(c[i++]);
            term.push_back(s),num.push_back(1);
        }
    }
}
typedef ll Mat[4][4];
Mat A,B,C;
void Mul(Mat &a,Mat b)
{
    Mat c;
    for(int i=0;i<4;i++)
        for(int j=0;j<4;j++)
        {
            c[i][j]=0;
            for(int k=0;k<4;k++)
                c[i][j]+=a[i][k]*b[k][j];
        }
    for(int i=0;i<4;i++)
        for(int j=0;j<4;j++)
            a[i][j]=(c[i][j]%mod+mod)%mod;
}
void Pow(Mat &a,ll b)
{
    if(b==1)return ;
    Mat c;
    for(int i=0;i<4;i++)
        for(int j=0;j<4;j++)
            c[i][j]=(i==j);
    while(b)
    {
        if(b&1)Mul(c,a);
        Mul(a,a);
        b>>=1;
    }
    for(int i=0;i<4;i++)
        for(int j=0;j<4;j++)
            a[i][j]=c[i][j];
}
ll Count()
{
    for(int i=0;i<4;i++)
        for(int j=0;j<4;j++)
            A[i][j]=(i==j);
    int sign=1;
    for(int k=0;k<term.size();k++)
    {
        string s=term[k];
        ll t=num[k];
        int n=s.size();
        for(int i=0;i<4;i++)
            for(int j=0;j<4;j++)
                B[i][j]=(i==j);
        for(int i=0;i<n;i++)
        {
            memset(C,0,sizeof(C));
            C[0][0]=C[3][3]=1;
            if(s[i]>='0'&&s[i]<='9')
                C[1][1]=10,C[2][1]=s[i]-'0',C[2][2]=1;
            else if(s[i]=='+')
                C[1][0]=sign,C[3][2]=1,sign=1;
            else if(s[i]=='-')
                C[1][0]=sign,C[3][2]=1,sign=-1;
            else C[1][2]=1;
            Mul(B,C); 
        }
        Mul(A,B);
        if(t>1)
        {
            for(int i=0;i<4;i++)
                for(int j=0;j<4;j++)
                    B[i][j]=(i==j);
            for(int i=0;i<n;i++)
            {
                memset(C,0,sizeof(C));
                C[0][0]=C[3][3]=1;
                if(s[i]>='0'&&s[i]<='9')
                    C[1][1]=10,C[2][1]=s[i]-'0',C[2][2]=1;
                else if(s[i]=='+')
                    C[1][0]=sign,C[3][2]=1,sign=1;
                else if(s[i]=='-')
                    C[1][0]=sign,C[3][2]=1,sign=-1;
                else C[1][2]=1;
                Mul(B,C); 
            }
            Pow(B,t-1);
            Mul(A,B);
        }
    }
    ll ans=(A[2][0]+A[3][0])%mod;
    if(sign==1)ans=(ans+A[2][1]+A[3][1])%mod;
    else ans=(ans-A[2][1]-A[3][1])%mod;
    ans=(ans+mod)%mod;
    return ans;
}
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%s",c);
        Deal();
        printf("%I64d\n",Count()); 
    }
    return 0;
}