牛客手速月赛43 E(dp)

题目
题意: 给定数字1-9的数量cnt[i]。计算出“满意的集合“的个数。
“满意的集合” 定义：选出的数的和为3的倍数。
空集合和为0，视为合法。此处的集合是本质set概念，123和321属于同一集合。
思路: dp.
定义f[i][j]表示只用1-i,将余数为j的集合的数量。
可以考虑暴力，枚举第i个数使用0-cnt[i]次，obviously，是一个不能接受的时间复杂度。
但是%3的结果只会有0、1、2，0-cnt[i]次中有多少个余0、1、2我们可以直接算出来。
选1、4、7个一样;2、5、8个一样;3、6、9个一样
那么我们可以直接O(1)计算出这三种模数，之后枚举前i-1个数的模数取值
时间复杂度: O(9 * 3)
代码:

// Problem: 满意的集合
// Contest: NowCoder
// URL: https://ac.nowcoder.com/acm/contest/11220/E
// Memory Limit: 524288 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<complex>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<unordered_map>
#include<list>
#include<set>
#include<queue>
#include<stack>
#define OldTomato ios::sync_with_stdio(false),cin.tie(nullptr),cout.tie(nullptr)
#define fir(i,a,b) for(int i=a;i<=b;++i) 
#define mem(a,x) memset(a,x,sizeof(a))
#define p_ priority_queue
// round() 四舍五入 ceil() 向上取整 floor() 向下取整
// lower_bound(a.begin(),a.end(),tmp,greater<ll>()) 第一个小于等于的
// #define int long long //QAQ
using namespace std;
typedef complex<double> CP;
typedef pair<int,int> PII;
typedef long long ll;
// typedef __int128 it;
const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const ll inf = 1e18;
const int N = 2e5+10;
const int M = 1e6+10;
const int mod = 1e9+7;
const double eps = 1e-6;
inline int lowbit(int x){ return x&(-x);}
template<typename T>void write(T x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x>9)
    {
        write(x/10);
    }
    putchar(x%10+'0');
}
template<typename T> void read(T &x)
{
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
int n,m,k,T;
int cnt[10]; 
ll f[10][3]; //只用1-i,余数为j 
void solve()
{
   n = 9;
   for(int i=1;i<=n;++i)
   {
   	  read(cnt[i]);
   }
   f[0][0] = 1;
   for(int i=1;i<=n;++i)
   {
   	 int mod1 = i % 3; int mod2 = i*2 % 3; int mod3 = 0; 
   	 ll sum1 = (cnt[i] + 2) / 3; // 1、4、7倍
   	 ll sum2 = (cnt[i] + 1) / 3; // 2、5、8倍
   	 ll sum3 = cnt[i] / 3 + 1; // 0、3、6、9倍
   	 for(int j=0;j<3;++j)
   	 {
   	   f[i][(mod1+j)%3] += sum1 * f[i-1][j]; //选sum1个i，使得余数从j -> j+mod1
   	   f[i][(mod2+j)%3] += sum2 * f[i-1][j];
   	   f[i][(mod3+j)%3] += sum3 * f[i-1][j];
   	   f[i][(mod1+j)%3] %= mod; f[i][(mod2+j)%3] %= mod; f[i][(mod3+j)%3] %= mod;
   	 }
   }
   write(f[n][0]);
}
signed main(void)
{  
   T = 1;
   // OldTomato; cin>>T;
   // read(T);
   while(T--)
   {
   	 solve();
   }
   return 0;
}