波动数列 背包变形 组合问题

观察这个数列：
1 3 0 2 -1 1 -2 …
这个数列中后一项总是比前一项增加2或者减少3，且每一项都为整数。
栋栋对这种数列很好奇，他想知道长度为 n 和为 s 而且后一项总是比前一项增加 a 或者减少 b 的整数数列可能有多少种呢？
输入格式
共一行，包含四个整数 n,s,a,b，含义如前面所述。
输出格式
共一行，包含一个整数，表示满足条件的方案数。
由于这个数很大，请输出方案数除以 100000007 的余数。
数据范围
1≤n≤1000
−10^9≤s≤ 10^9
1≤a,b≤10^6
输入样例：
4 10 2 3
输出样例：
2
样例解释
两个满足条件的数列分别是2 4 1 3和7 4 1 -2。

分析：
1：首项可以是任意的数，但会影响接下来的数，后一项总是比前一项增加a或者减少b
设首项是x，第二项为x+1×d1,第三项为x+2×d2，第n项为x+（n-1)×dn-1；
其中di是 +a和-b的一个随机随机值。
2 ：由总和是s 得：nx+1×d1+2×d2+……（n-1）×dn-1=s
变形后：x=s-（d1+2×d2+……（n-1）×dn-1）/n
由于：x一定是整数
即：s%n=（d1+2*d2+……（n-1）×dn-1）%n
3：s是已知的，故题目是满足所有要求的di值，求%n和s%n相同的数量即是方案数。

状态表示：
设f[i][j]是前i个数，且当前总和%n的余数是j的方案数；
状态转移：
其实di有两种选法，不停的做选择直到最终状态，有2的n次方的总数，类似背包问题，每个物品有选和不选的两种状态。实质是一类组合问题。
对于当前的第i个数，它和前面的数息息相关，而前面第i-1个数要想到第i，有两种路径，一种是+a，一种是-b，第i的方案数是前面两个的和；而s%n=j，故j可代表当前的总数，（j-i×a）%n和（j+i×b）%n是第i-1和%n的余数，是由它们两个+a或者-b得到i的;

即：f[i][j]=f[i-1][（j-i×a）%n]+f[i-1][（j+i×b）%n]

初始化：
f[0][0]=1 一个数都没有只有一种方案。

#include <iostream>
using namespace std;
const int N=1005,MOD=100000007;
int f[N][N]={0};
int get_mod(int a,int b)  //得到a%b的正余数 
{
 return (a%b+b)%b;
}
int main()
{
   int n,s,a,b;
   cin>>n>>s>>a>>b;
   f[0][0]=1;
   for(int i=1;i<n;i++)
   {
      for(int j=0;j<n;j++)
   {
    f[i][j]=(f[i-1][get_mod(j-i*a,n)]+f[i-1][get_mod(j+i*b,n)])%MOD;
   } 
   } 
   cout<<f[n-1][s%n]<<endl;
 return 0;
} 

由于中间取余可能出现负数，要添加以下代码：

int get_mod(int a,int b)  //得到a%b的正余数 
{
 return (a%b+b)%b;
}

还有一点是：
最后的输出不是f[i][s%n],而是f[i-1][s%n],最后一项是第n-1项，默认忽视了第一项。