D.Modulo Nine (DP)

思路：若要a1*a2*a3....*an能mod9等于0，则该算式中必有两个因子3相乘或者含有0。对于一个右端点r，有多个[l  ,r]的时候，只要最小的那个区间满足，其他肯定都满足了。这样对于一个右端点的限制条件是否满足，取决于前面出现的最后两个3出现的位置（9和0算作两个3）
由此我们设dp[i][j][k]表示前i个数，最后一个因子3出现的位置为j，倒数第二个为k。由于因子中含有3的数有3、6、9，其中9中含有两个因子3。因此，只需要一个0或9，或者两个3、6即可满足条件。

则方案数的计算分为三种情况：

1. 不含因子3
2. 含有一个因子3（3，6）
3. 含有两个因子3（0，9）

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll long long
using namespace std;
ll n,m,dp[55][55][55];
ll l[55];
const ll mod=1e9+7;
int main(){
    while(~scanf("%lld%lld",&n,&m)){
        ll a,b;
        memset(dp, 0, sizeof(dp));
        memset(l, 0, sizeof(l));
        while(m--){
            cin>>a>>b;
            l[b]=max(l[b],a);
        }
        dp[0][0][0]=1;
        for(int i=0;i<n;i++){
            for(int j=0;j<=i;j++){
                for(int k=0;k<=j;k++){
                    if(l[i]<=k){//使得右端点对应的最小的区间满足条件
                        dp[i+1][j][k]=(dp[i][j][k]*6+dp[i+1][j][k])%mod;//不含因子3的方法数
                        dp[i+1][i+1][j]=(dp[i][j][k]*2+dp[i+1][i+1][j])%mod;//含有一个因子3的方法数
                        dp[i+1][i+1][i+1]=(dp[i][j][k]*2+dp[i+1][i+1][i+1])%mod;//含有两个因子3的方法数
                    }
                }
            }
        }
        ll ans=0;
        for(int i=0;i<=n;i++){
            for(int j=0;j<=i;j++){
                if(l[n]<=j)
                ans=(ans+dp[n][i][j])%mod;
            }
        }
        cout<<ans<<endl;
    }
    
    
    return 0;
}