ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Ant Counting POJ 3046(计数dp)

2019-08-25 10:52:10  阅读:191  来源: 互联网

标签:now Counting const int LL Ant POJ include dp


原题

题目链接

题目分析

计数dp题,感觉其实也可以用组合数学做,但我太菜了,推不出通用公式.dp可以定义dp[i][j]为前i种选j个蚂蚁有多少种选法,然后递推公式如下,其中c[i]表示第i种的数量,

dp[i][j]=Σ(min(j,c[i]),k=0)dp[i-1][j-k].可以化简一下,dp[i][j]=∑(c[i],k=0)dp[i-1][j-1-k]+dp[i-1][j]-dp[i-1][j-1-c[i]],最后答案就是dp[t][s]+dp[t][s+1]+...+dp[t][b].然后这道题要状压,也就是i只有i和i-1,所以第一维可以压缩成i=1和i=0.

代码

 1 #include <iostream>
 2 #include <algorithm>
 3 #include <utility>
 4 #include <cstdio>
 5 #include <cmath>
 6 #include <cstring>
 7 #include <string>
 8 #include <vector>
 9 #include <stack>
10 #include <queue>
11 #include <map>
12 #include <set>
13 
14 using namespace std;
15 typedef long long LL;
16 const int INF_INT=0x3f3f3f3f;
17 const LL INF_LL=0x3f3f3f3f3f3f3f3f;
18 
19 const int mod=1e6;
20 int num[2000];
21 LL dp[2][200000];
22 
23 int main()
24 {
25 //    freopen("black.in","r",stdin);
26 //    freopen("black.out","w",stdout);
27     int t,a,s,b;
28     cin>>t>>a>>s>>b;
29     for(int i=0;i<a;i++)
30     {
31         int x;
32         scanf("%d",&x);
33         num[x]++;
34     }
35     int total=0;
36     for(int i=1;i<=t;i++)
37     {
38         total+=num[i];
39         int now=i&1,last=!now;
40         dp[now][0]=dp[last][0]=1;
41         for(int j=1;j<=total;j++)
42         {
43             if(j-1-num[i]>=0) dp[now][j]=(dp[now][j-1]+dp[last][j]-dp[last][j-1-num[i]])%mod;
44             else dp[now][j]=(dp[now][j-1]+dp[last][j])%mod;
45         }
46     }
47     LL ans=0;
48     int res=t&1;
49     for(int i=s;i<=b;i++) ans=(ans+dp[res][i])%mod;
50     cout<<(ans+mod)%mod<<endl;
51     return 0;
52 }

 

 

标签:now,Counting,const,int,LL,Ant,POJ,include,dp
来源: https://www.cnblogs.com/VBEL/p/11407147.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有