AcWing 1214. 波动数列

摆烂了，有问题欢迎评论

#include <iostream>

using namespace std;

int f[1010][1010];
int n,s,a,b;
int N=100000007;

int mod(int x,int y) {
    return (x%y+y)%y;
}
int main() {
    cin >> n >> s >> a >> b;
    f[0][0]=1;
    for(int i=0; i<n; i++)
        for(int j=0; j<n; j++) {
            f[i+1][mod(j+(n-i-1)*a,n)]=(f[i+1][mod(j+(n-i-1)*a,n)]+f[i][j])%N;
            f[i+1][mod(j-(n-i-1)*b,n)]=(f[i+1][mod(j-(n-i-1)*b,n)]+f[i][j])%N;
        }
    cout << f[n-1][mod(s,n)];
    return 0;
}