码迷,mamicode.com
首页 > 其他好文 > 详细

Dollar Dayz poj3181

时间:2015-10-14 18:00:37      阅读:153      评论:0      收藏:0      [点我收藏+]

标签:

http://poj.org/problem?id=3181

这个题目一开始就能看出来是个dp问题,但是我并没有一开始就看出来是一个完全背包为题,只是想着根据以前的方法,这个问题应该是可以找到规律的,但是,还是被坑了,这还是一个大数问题!

首先我膜拜一下hankcs大神的:

///////////////////////////////////////////////////////////

  1. #include <iostream>
  2. using namespace std;
  3.  
  4. unsigned long long dp[100 + 16][1000 + 16]; // dp[i][j] := 用i种价格配出金额j的方案数
  5.  
  6. ///////////////////////////SubMain//////////////////////////////////
  7. int main(int argc, char *argv[])
  8. {
  9. #ifndef ONLINE_JUDGE
  10.     freopen("in.txt", "r", stdin);
  11.     freopen("out.txt", "w", stdout);
  12. #endif
  13. int N, K;
  14. cin >> N >> K;
  15. dp[0][0] = 1;
  16. for (int i = 1; i <= K; ++i)
  17. {
  18. for (int k = 0; k <= N; k += i)
  19. {
  20. for (int j = N; j >= k; --j)
  21. {
  22. dp[i][j] += dp[- 1][- k];
  23. }
  24. }
  25. }
  26. cout << dp[K][N] << endl;
  27. #ifndef ONLINE_JUDGE
  28.     fclose(stdin);
  29.     fclose(stdout);
  30.     system("out.txt");
  31. #endif
  32.     return 0;
  33. }
  34. ////////////////////////////////////////////////////////////

  hancks的这个做法是用完全背包

  dp[i][j] = dp[i – 1][j] + dp[i – 1][j – i] + dp[i – 1][j – 2 * i] + … + dp[i – 1][0]

  由这个公式可以再递推:  

将j换成j – i有

dp[i][j – i] = dp[i – 1][j – i] + dp[i – 1][j – 2 * i] + … + dp[i – 1][0]

得出:if j >= i:

  dp[i][j] = dp[i-1][j] + dp[i][j-i];

  我的做法是一开始就推出了这个公式,因为不小心看出了这个规律

    i:1->4 ,j :1->5 dp[i][j]规律是这样的

    1 1 1 1 1

    1 2 2 3 3

    1 2 3 4 5

    1 2 3 4 6     

    得出了j >= i : dp[i][j] = dp[i-1][j] + dp[i][j-i] 

    然而,这还是个大数问题,即使unsigned long long 都不行,开始一直没想通!

/*************************************************************************
    > File Name: DollarDayz_poj3181.cpp
    > Author: spzhao
    > Mail: spzhaol@163.com 
    > Created Time: 2015年10月14日 星期三 11时13分22秒
 ************************************************************************/

#include<iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

#define mod 10000000000000000
using namespace std;

const int N = 1005;
const int K = 105;
unsigned long long dp[100+16][1000+16][2];
int n,k;
void solve()
{
	for (int i = 1;i <= k;i++)
	{
		for (int j = 1;j <= n;j++)
		{
			if (j >= i)
			{
				dp[i][j][0] = dp[i-1][j][0]+dp[i][j-i][0];
				dp[i][j][1] = dp[i-1][j][1]+dp[i][j-i][1];
				dp[i][j][0] += dp[i][j][1]/mod;
				dp[i][j][1] = dp[i][j][1]%mod;
			}
			else
			{
				dp[i][j][0] = dp[i-1][j][0];
				dp[i][j][1] = dp[i-1][j][1];
			}
		}
	}
	if (dp[k][n][0])
		cout << dp[k][n][0];
	cout << dp[k][n][1] << endl;
}
int main ()
{
		cin >> n >> k;
		memset(dp,0,sizeof(dp));
		dp[1][0][1] = 1;
		for (int i = 1;i <= k;i++)
			dp[i][0][1] = 1;
		solve();
	return 0;
}

  

 

 

 

 

 

Dollar Dayz poj3181

标签:

原文地址:http://www.cnblogs.com/ediszhao/p/4877997.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!