失踪OJ回归。

  小C通过这道题mark一下容斥一类的问题。

Description

  硬币购物一共有4种硬币。面值分别为c1,c2,c3,c4。某人去商店买东西,去了tot次。每次带di枚ci硬币,买s的价值的东西。请问每次有多少种付款方法。

Input

  第一行 c1,c2,c3,c4,tot 下面tot行 d1,d2,d3,d4,s。

Output

  每次的方法数。

Sample Input

  1 2 5 10 2
  3 2 3 1 10
  1000 2 2 2 900

Sample Output

  4
  27

HINT

  di,s<=100000,tot<=1000。

Solution

  O(s*tot)的DP算法谁都会写,但是看着这时间复杂度你难道不虚吗?

  你难道甘愿被卡常而就此丢掉10分或是在刷题时打一个这样的暴力草草了事而对其中的精妙不闻不问,你的良心不会痛吗?

  如果你在解题时看到题目中有常数a(a<20)这样的数据,不妨就往O(2a)这样复杂度的算法去想一想。

  因为题目中有一个常数4,而且还都是限制条件,所以我们就往状压、容斥这方面去想。

  显然状压是不可能的,于是我们就只有容斥了。

  首先我们要知道容斥在这道题是干嘛用的:

  容斥就是对多个限制条件下方案的去重工作,也就是你们所熟知的,求多个集合的并集。

  容斥常常伴随的思想是一种逆向思维,就是题目往往要求我们去求多个集合的交,然而我们并没有好的办法,转而跑去求各个集合的补的并,再补回来就是各个集合的交。

  我们需要理解在这个算法中, 交 是可以O(1)求得的,然而 并 需要用容斥求得。

  所以解这样的题目的大致思路就是:

    题目要求我们求A1~An的交,但是我们发现很难求;

    所以我们去求CuA1~CuAn的并,而CuA1~CuAn的并需要我们求CuA1~CuAn的交,但是我们发现CuA1~CuAn的交特别好求,所以我们就圆满地解决了这个问题。

  运用这样的思路,我们就可以很快解出这道题。

  所以我们要求的只剩,在这次购物中,第 i 种硬币用了超过di的方案数。

  脑补一下,我们就知道,我们只要每种硬币先取到它们限制的数量+1,剩下随便取就可以了,这样的方案无论如何都是满足每种硬币都超过限制的。

  所以用一句话概括题解:容斥,求f[s-Σ(ci*(di+1))],f[x]为在没有任何限制下,取硬币得到面值x的方案数。

  时间复杂度O(max(s)*4+tot*16)。注意答案最大为C(s,4)会爆int。

#include <cstdio>
#include <cstring>
#include <algorithm>
#define ll long long
#define MM 100005
using namespace std;
int n,m;
int a[],g[],ys[];
ll f[MM],lt,ans;
bool u; inline int read()
{
int n=,f=; char c=getchar();
while (c<'' || c>'') {if(c=='-')f=-; c=getchar();}
while (c>='' && c<='') {n=n*+c-''; c=getchar();}
return n*f;
} int main()
{
register int i,j;
a[]=read(); a[]=read(); a[]=read(); a[]=read(); n=read();
f[]=ys[]=; ys[]=; ys[]=; ys[]=;
for (i=;i<=;++i)
for (j=a[i];j<MM;++j) f[j]+=f[j-a[i]];
while (n--)
{
g[]=read(); g[]=read(); g[]=read(); g[]=read(); m=read();
ans=;
for (i=;i<;++i)
{
for (u=lt=,j=;j<=;++j)
if (i&ys[j]) u^=,lt+=1LL*a[j]*(g[j]+);
if (lt>m) continue;
ans+=(f[m-lt])*(u?-:);
}
printf("%lld\n",ans);
}
}

Last Word

  大概就是小C关于容斥的一点点想法,当然这样的题目还有很多无法以偏概全。希望这样的思路能对以后有一点帮助吧。

[BZOJ]1042 硬币购物(HAOI2008)的更多相关文章

  1. BZOJ 1042 硬币购物(完全背包+DP)

    题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=1042 题意:给出四种面值的硬币c1,c2,c3,c4.n个询问.每次询问用d1.d2.d ...

  2. BZOJ 1042 硬币购物(背包DP+容斥原理)

    可以看出这是个多重背包,运用单调队列优化可以使每次询问达到O(s).这样总复杂度为O(s*tot). 会TLE. 因为改题的特殊性,每个硬币的币值是不变的,变的只是每次询问的硬币个数. 我们不妨不考虑 ...

  3. BZOJ 1042 硬币购物

    先不考虑限制,那么有dp[i]表示i元钱的方案数. 然后考虑限制,发现可以容斥. 其实整个题就是两个容斥原理.感觉出的蛮好的. #include<iostream> #include< ...

  4. 【BZOJ】【1042】【HAOI2008】硬币购物

    DP+容斥原理 sigh……就差一点…… 四种硬币的数量限制就是四个条件,满足条件1的方案集合为A,满足条件2的方案集合为B……我们要求的就是同时满足四个条件的方案集合$A\bigcap B\bigc ...

  5. Bzoj 1042: [HAOI2008]硬币购物 容斥原理,动态规划,背包dp

    1042: [HAOI2008]硬币购物 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 1747  Solved: 1015[Submit][Stat ...

  6. bzoj 1042: [HAOI2008]硬币购物 dp+容斥原理

    题目链接 1042: [HAOI2008]硬币购物 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 1706  Solved: 985[Submit][ ...

  7. BZOJ 1042: [HAOI2008]硬币购物( 背包dp + 容斥原理 )

    先按完全背包做一次dp, dp(x)表示x元的东西有多少种方案, 然后再容斥一下. ---------------------------------------------------------- ...

  8. BZOJ 1042: [HAOI2008]硬币购物 [容斥原理]

    1042: [HAOI2008]硬币购物 题意:4种硬币.面值分别为c1,c2,c3,c4.1000次询问每种硬币di个,凑出\(s\le 10^5\)的方案数 完全背包方案数? 询问太多了 看了题解 ...

  9. 【BZOJ】1042: [HAOI2008]硬币购物

    1042: [HAOI2008]硬币购物 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 3307  Solved: 2075[Submit][Stat ...

随机推荐

  1. 亚马逊AWS学习——EC2的自定义VPC配置

    1 网络配置 EC2即亚马逊AWS云服务中的虚拟主机.创建EC2实例时如果使用的默认VPC并分配了公有IP是可以上网的.但我们经常需要自定义的网络环境,这时就需要自己定义VPC和子网了. 1.1 配置 ...

  2. 我从业11年来遇到的最奇葩的raid0+1数据恢复经历

    我是一名数据恢复工程师,从事数据恢复行业已经11年了,前几天接到一组4块盘SCSI RAID0+1的数据恢复,客户说做了两组raid1,现在raid状态里显示有3快盘offline.如果两组盘分别作r ...

  3. Python内置函数(34)——filter

    英文文档: filter(function, iterable) Construct an iterator from those elements of iterable for which fun ...

  4. SpringCloud的Config:ConfigServer注册到EurekaServer中,变成一个Eureka服务

    一.概念与定义 1.将SpringCloud ConfigServer注册到 EurekaServer,以便ConfigClient以服务的方式引用ConfigServer 2.客户端不再引用 Con ...

  5. OAuth2.0学习(1-4)授权方式1-授权码模式(authorization code)

    参与者列表: (1) Third-party application:第三方应用程序,又称客户端(client),如:"云冲印".社交应用. (2)HTTP service:HTT ...

  6. 基于DFS的拓扑排序

    传送门:Kahn算法拓扑排序 摘录一段维基百科上的伪码: L ← Empty list that will contain the sorted nodes S ← Set of all nodes ...

  7. MySQL 5.7 基于复制线程SQL_Thread加快恢复的尝试

    1. MySQL 数据恢复常用办法 MySQL恢复的方法一般有三种: 1. 官方推荐的基于全备+binlog , 通常做法是先恢复最近一次的全备,然后通过mysqlbiinlog --start-po ...

  8. Git reset到某一次commit

    下图场景:张三和李四并行开发,张三提交commit1(bc2dd00),李四提交commit2(7f019d2),张三再提交commit3(44d4fc5),如果此时李四revert commit2, ...

  9. js变量的生命周期

    1.在JavaScript中,对于for循环中定义的i变量,其生命周期在循环结束后仍然是有效的. for (var i=0; i < 10; i++){ doSomething(i); } al ...

  10. java 反射(Reflection)

    看了很多关于java 反射的文章,自己把所看到的总结一下.对自己,对他人或多或少有帮助吧. Java Reflection是什么? 首先来看看官方文档Oracle里面对Reflection的描述: R ...