失踪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. 2017北京国庆刷题Day4 morning

    期望得分:0+40+30=70 实际得分:0+10+10=20 题目修改:只能由0变1,只能用一次操作 大模拟 #include<cstdio> #include<cstring&g ...

  2. 通过URL传递PDF名称参数显示PDF

    1 <%@ page language="java" import="java.util.*,java.io.*" 2 pageEncoding=&quo ...

  3. 算法题丨3Sum Closest

    描述 Given an array S of n integers, find three integers in S such that the sum is closest to a given ...

  4. Linux入门:vi 和 vim

    vi编辑器是所有Unix及Linux系统下标准的编辑器,它的强大不逊色于任何最新的文本编辑器. 本文介绍了vi (vim)的基本使用方法,但对于普通用户来说基本上够了! 转自:http://www.l ...

  5. C#制作ActiveX插件

    首先新建项目--->类库,取名:ActiveXDemo 右键项目属性:应用属性==>程序集信息=>使程序集Com可见, 生成==>输出==>为com互操作注册 新建接口类 ...

  6. Stanford依存句法关系解释

    ROOT:要处理文本的语句 IP:简单从句 NP:名词短语 VP:动词短语 PU:断句符,通常是句号.问号.感叹号等标点符号 LCP:方位词短语 PP:介词短语 CP:由'的'构成的表示修饰性关系的短 ...

  7. python3.6执行pip3时 Unable to create process using '"'

    问题需求 由于在windows操作系统中已经安装了python2.7,要在安装python3的时候 将python3.6安装在C:\Python36目录下 然后进入C:\Python36目录下执行pi ...

  8. Django REST framework+Vue 打造生鲜超市(五)

    六.商品类别数据展示 6.1. 商品类别数据接口 (1)商品分类有两个接口: 一种是全部分类:一级二级三级 一种是某一类的分类以及商品详细信息: 开始写商品分类的接口 (2)序列化 给分类添加三级分类 ...

  9. [转]pymongo常用操作函数

    pymongo 是 mongodb 的 python Driver Editor.记录下学习过程中感觉以后会常用多一些部分,以做参考. 1. 连接数据库 要使用pymongo最先应该做的事就是先连上运 ...

  10. Oracle SQL Developer 免费的DB2客户端

    问题地址:https://stackoverflow.com/questions/8600735/is-there-any-opensource-db2-client 软件地址:http://www. ...