动态规划 01背包 poj 1745

来源:互联网 发布:电脑淘宝商品链接 编辑:程序博客网 时间:2024/05/18 13:43
//初始值 F[0][0]=1;//状态方程  m = N[i] % k ; F[i][j] = F[i-1][ (j+m+k) %k] || F[i-1] [ (j-m+k) %k] //+k是为了不出现负坐标

#include <iostream>#include <stdio.h>#include <cstring>using namespace std;#define MAX_N (10000+1)#define MAX_K (100+1)int N[MAX_N];int K[MAX_K];bool F[MAX_N][MAX_K]; //F[i][j] 定义,前i个数通过+,-组合 %k的值能否为jint n,k;#define printf //void dump(){for(int i=0;i<=n;i++){printf("i=%d\t",i);for(int j=0;j<k;j++)printf("%4d ",F[i][j]);printf("\n");}}int get_mod(int A){//保证结果为正return ((A%k)+k)%k;}//初始值 F[0][0]=1;//状态方程  m = N[i] % k ; F[i][j] = F[i-1][ (j+m+k) %k] || F[i-1] [ (j-m+k) %k] //+k是为了不出现负坐标int main(){while(cin>>n>>k){for(int i=1;i<=n;i++)cin>>N[i];memset(F,0,sizeof(F));//F[0][0]=true;F[1][get_mod(N[1])] = true;for(int i=2;i<=n;i++) // n-1{int m = get_mod(N[i]); //这里注意N[i]%k结果可能为负,封装函数保证结果为正printf("i=%d N[i] = %d k=%d n%%k=%d m=%d\n",i,N[i],N[i]%k,k,m);for(int j=0;j<k;j++){int n1 = get_mod(j-m) ,n2 = get_mod(j+m);F[i][j] = F[i-1][n1]  || F[i-1] [ n2] ;printf("i=%d j=%d n1 =%d,n2=%d F[%d] = %d,F[%d] =%d->F[%d][%d]=%d\n",i,j,n1,n2,n1,F[i-1][n1],n2,F[i-1][n2],i,j,F[i][j]);}}if(F[n][0])cout<<"Divisible"<<endl;elsecout<<"Not divisible"<<endl;dump();}}