#include<cstdio>
int a[51];
using namespace std;
int dp[51][3];
int main(){
	char s[51];
	{
	scanf("%s",s+1);
	int len=strlen(s);
      for (int i=0;i<len;i++) 
	   a[i]=s[i]-'0';
	}
	for (int i=1;i<=len;i++)
	  {
	  if (a[i]%3==0) {
	      dp[i][0]=dp[i-1][0]+1;
          dp[i][1]=dp[i-1][1]+1;
          dp[i][2]=dp[i-1][2]+1;
	   }
	  if (a[i]&&a[i]%3==1) {
	      dp[i][2]=dp[i-1][1];
	      dp[i][0]=dp[i-1][2]+1;
	      dp[i][1]=dp[i-1][0]+1; 
      }
	  if (a[i]&&a[i]%3==2) dp[i][1]=dp[i-1][1]+1;
	  }
	  printf("%d",dp[len][0]);
return 0;
}

写了一半右崩了,本来是想按照区间长度和余数组合进行dp,结果gg,写不下去

根据题解又错了
gg,字符串处理和dp的边界没有处理到
#include
using namespace std;
char s[1000];
int main(){
scanf("%d",&s);
for (int i=strlen(s)-1;i>=0;i–)
{
for (int j=i;j<strlen(s)-1;j++)
if (s[j+1]%3==0){
dp[i][j+1][0]=dp[i][j]
}
}
}
正解

int mod = 1e9+7;
 
int main()
{
    int dp[55][55][5];//dp[i][j][k]表示在区间(i,j)中被3整除余数为k的个数
    memset(dp,0,sizeof(dp));
    string s;
    cin>>s;
    for(int i=s.length();i>=1;i--) s[i] = s[i-1];
    for(int i=1;i<=s.length();i++) dp[i][i][(s[i]-'0')%3] = 1;
 
    for(int i=s.length();i>=1;i--){
        for(int j=i;j<s.length();j++){
            //s[j+1]为新加入的值
            if(s[j+1]%3 == 0){
                //如果新加入的值是3的倍数,
                //则(i,j)模3余k的个数除了之前的个数;
                //每个都可以加上s[j+1],即dp[i][j][k]*2;
                //而余0的情况还要加上s[j+1]本身。
                dp[i][j+1][0] = 2*dp[i][j][0]+1;
                dp[i][j+1][1] = 2*dp[i][j][1];
                dp[i][j+1][2] = 2*dp[i][j][2];
            }
            else if((s[j+1]-'0')%3 == 1){
                //如果新加入的值是模3余1,
                //则(i,j)模3余k的个数除了之前的个数;
                //还要加上之前模3余k-1的个数;
                //而余1的情况还要加上s[j+1]本身。
                dp[i][j+1][0] = dp[i][j][0] + dp[i][j][2];
                dp[i][j+1][1] = dp[i][j][1] + dp[i][j][0] + 1;
                dp[i][j+1][2] = dp[i][j][2] + dp[i][j][1];
            }
            else{
                //如果新加入的值是模3余2,
                //则(i,j)模3余k的个数除了之前的个数;
                //还要加上之前模3余k-2的个数;
                //而余2的情况还要加上s[j+1]本身。
                dp[i][j+1][0] = dp[i][j][0] + dp[i][j][1];
                dp[i][j+1][1] = dp[i][j][1] + dp[i][j][2];
                dp[i][j+1][2] = dp[i][j][2] + dp[i][j][0] + 1;
            }
            for(int k=0;k<3;k++) dp[i][j+1][k] = dp[i][j+1][k] % mod;
        }
    }
    cout<<dp[1][s.length()][0]<<endl;
    return 0;
}