1307E - Cow and Treats

题意

有一排给定口味的草,并且给m头牛,每个牛都只吃一种口味的草,并且要吃给定数量个。现在可以安排牛从两边出发,方向向另一方向进发,每次路过符合他口味的草就会吃,吃到他需要的数量就会停下开始睡觉,开始睡觉后就是个路障,不能通行。如果一个牛吃到了他需要的数量的草就会Happy,让你安排牛从两边走,求最大数量的happy牛的数量,并且求出满足最大数量情况下的方案数。方案的不同体现在集合的不同,比如从左边出发的牛集合不同,就算做不同方案。

题解

首先我们先来解决最大数量的问题。对于吃每一种口味的牛来说,同口味的牛,不能从同方向出发超过两只,显然可得,那么也就是说同一个口味的牛最多左边右边出发一只,也就是最多2只。对每一种牛求一下即可得最大口味数,本题的难点在于求方案数。这里要分左右n只有5000,一个很自然的想法是枚举从左和从右的分界线。假如从左开始到分界线一共有right[i]个i口味的草,那么需要吃的数量小于等于right[i]的牛都可以从这里出发,左边同理。我们设从左出发满足条件的牛有a个,从右有b个,那么从两边方案数量为a*b-min(a,b) 这里用的是容斥的思想。因为min(a,b),一定是max(a,b)的子集,所以可得。A交B的数量为min(a,b),从一边的方案数为(a+b)。从左到右枚举分界点算一遍即可。这里还要解决的是重复的问题,例如[1,10]区间假如[1,4] [7,10] 分别是两边行进的路径,那么[1,5][5,10] 就会把这个集合再算一遍,就重复了。怎么解决这个问题呢,我们可以考虑枚举每一个有牛休息的点,对于这个点,固定该点口味的牛,其他的点只要符合要求随意停留,这样就能不重不漏计算出方案数了。

这里可以用upper_bound 很方便得找到符合条件的数量,同时用std::binary_search可以方便判断一个数在不在数组中。

代码

#include<bits/stdc++.h>
#define pb push_back
#define mkp make_pair
using namespace std;
typedef long long ll;
const int mod=1e9+7;
const int maxn=5000+5;
int lefts[maxn],rights[maxn],num[maxn],ways[maxn],s[maxn];
int n,m;
vector<int>v[maxn];
void solve1(int id){
int a=upper_bound(v[id].begin(),v[id].end(),rights[id])-v[id].begin();
int b=upper_bound(v[id].begin(),v[id].end(),lefts[id])-v[id].begin();
ll cnt2=a*b-min(a,b);//tow sides
ll cnt1=a+b;//one sides
if(cnt2>0){
ways[id]=cnt2;
num[id]=2;
}
else if(cnt1>0){
ways[id]=cnt1;
num[id]=1;
}
else {
ways[id]=1;
num[id]=0;
}
}
ll ans_tot=-1,ans_ways=1;
ll add(ll a,ll b ){
a=(a+mod)%mod,b=(b+mod)%mod;
return (a+b+mod)%mod;
}
void update(ll tmp_tot,ll tmp_ways){
if(tmp_tot>ans_tot){
ans_tot=tmp_tot;
ans_ways=tmp_ways;
}
else if(tmp_tot==ans_tot){
ans_ways=add(ans_ways,tmp_ways);
}
}
void solve2(int id){
int a=upper_bound(v[id].begin(),v[id].end(),rights[id])-v[id].begin();
if(rights[id]>=lefts[id])//因为同一种口味同1数量的牛不可能有2个 所以固定一个在边界睡觉后 如果right[id]>=left[id] 既这个牛也满足从左边出发,只需要减掉他即可
a--;
if(a>0){
ways[id]=a;
num[id]=2;
}
else {
ways[id]=1;
num[id]=1;
}
}
ll mul(ll a,ll b){
a=(a+mod)%mod;
b=(b+mod)%mod;
return (a%mod)*(b%mod)%mod;
}
ll fpow(ll a,ll b){
ll ans=1;
while(b){
if(b&1)ans=mul(ans,a);
a=mul(a,a);
b>>=1;
}
return ans;
}
ll inv(ll x){
return fpow(x,mod-2);
}
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)scanf("%d",&s[i]),rights[s[i]]++;
for(int i=1;i<=m;i++){
int x,y;
scanf("%d%d",&x,&y);
v[x].push_back(y);
}
for(int i=1;i<=m;i++)sort(v[i].begin(),v[i].end());
for(int i=1;i<=n;i++){
solve1(i);
}
ll tmp_tot=0,tmp_ways=1;
for(int i=1;i<=n;i++){
tmp_tot+=num[i];
tmp_ways=mul(tmp_ways,ways[i]);
}
update(tmp_tot,tmp_ways);
//cout<<ans_tot<<" "<<ans_ways<<endl;
for(int i=1;i<=n;i++){
tmp_tot-=num[s[i]];
tmp_ways=mul(tmp_ways,inv(ways[s[i]]));
//cout<<"??? "<<tmp_ways<<endl;
rights[s[i]]--,lefts[s[i]]++;
if(binary_search(v[s[i]].begin(),v[s[i]].end(),lefts[s[i]])){
solve2(s[i]);
update(tmp_tot+num[s[i]],mul(tmp_ways,ways[s[i]])); //cout<<i<<" "<<lefts[s[i]]<<" "<<tmp_tot+num[s[i]]<<" "<<ways[s[i]]<<endl;
//cout<<ans_tot<<" zzz "<<ans_ways<<endl;
}
solve1(s[i]);
tmp_tot+=num[s[i]];
//cout<<ways[s[i]]<<endl;
tmp_ways=mul(tmp_ways,ways[s[i]]);
}
printf("%lld %lld\n",ans_tot,ans_ways);
return 0;
}

1307E - Cow and Treats 二分 枚举边界 容斥的更多相关文章

  1. hdu-2197 本原串---枚举因子+容斥定理

    题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=2197 题目大意: 由0和1组成的串中,不能表示为由几个相同的较小的串连接成的串,称为本原串,有多少个 ...

  2. 牛客练习赛43 Tachibana Kanade Loves Game (简单容斥)

    链接:https://ac.nowcoder.com/acm/contest/548/F来源:牛客网 题目描述 立华奏是一个天天打比赛的萌新. 省选将至,萌新立华奏深知自己没有希望进入省队,因此开始颓 ...

  3. POJ 3189——Steady Cow Assignment——————【多重匹配、二分枚举区间长度】

     Steady Cow Assignment Time Limit:1000MS     Memory Limit:65536KB     64bit IO Format:%I64d & %I ...

  4. YYHS-分数(二分+容斥)

    题目描述 KJDH是个十分善于探索的孩子,有一天他把分子分母小于等于n的最简分数列在了纸上,他想找到这些分数里第k小的数,这对于KJDH来说当然是非常轻易,但是KJDH最近多了很多妹子,他还要去找妹子 ...

  5. 第k个互质数(二分 + 容斥)

    描述两个数的a,b的gcd为1,即a,b互质,现在给你一个数m,你知道与它互质的第k个数是多少吗?与m互质的数按照升序排列. 输入 输入m ,k (1<=m<=1000000;1<= ...

  6. Codeforces 801C Voltage Keepsake(二分枚举+浮点(模板))

    题目链接:http://codeforces.com/contest/801/problem/C 题目大意:给你一些电器以及他们的功率,还有一个功率一定的充电器可以给这些电器中的任意一个充电,并且不计 ...

  7. nyoj762——分解质因数+容斥+二分

    第k个互质数 时间限制:1000 ms  |  内存限制:65535 KB 难度:4   描述 两个数的a,b的gcd为1,即a,b互质,现在给你一个数m,你知道与它互质的第k个数是多少吗?与m互质的 ...

  8. 【周期性/容斥+二分】POJ2773-HAPPY 2006

    [题目大意] 求与n互质的第k个数. [思路] 先求出小于k且与n互质的数,再利用gcd(bt+a,b)=gcd(a,b)的性质求解,效率低.枚举与n互质的数的效率是O(nlogn),求解第k个数的效 ...

  9. POJ 2112—— Optimal Milking——————【多重匹配、二分枚举答案、floyd预处理】

    Optimal Milking Time Limit:2000MS     Memory Limit:30000KB     64bit IO Format:%I64d & %I64u Sub ...

随机推荐

  1. ES6 - 基础学习(6): 对象扩展

    对象对于JavaScript至关重要,在ES6中对象又加了很多新特性. 对象字面量:属性的简洁表示法 ES6允许对象的属性直接写变量,这时候属性名是变量名,属性值是变量值. let attr1 = & ...

  2. Java虚拟机——JVM

    一.JVM整体架构 1.JVM(Java虚拟机):指以软件的方式模拟具有完整硬件系统功能.运行在一个完全隔离环境中的完整计算机系统,是物理机的软件实现.常用的虚拟机有VMWare.Virtual Bo ...

  3. 基于Jupyter Notebooks的C# .NET Interactive安装与使用

    .NET Interactive发布预览版了,可以像Python那样用jupyter notebooks来编辑C#代码.具体可以在GitHub上查看dotnet/interactive项目. 安装步骤 ...

  4. Python 编程入门(2):复杂数据类型(列表,字典)

    以下所有例子都基于最新版本的 Python,为了便于消化,每一篇都尽量短小精悍,希望你能尽力去掌握 Python 编程的「概念」,可以的话去动手试一下这些例子(就算目前还没完全搞懂),加深理解. 在 ...

  5. (未完成)【Android】MVP模式初见(一)

    最近在阅读郭霖大神的公众号时,分类中架构引起了我的注意. 虽然是个人开发(水平很菜的那种),但最终都要向企业正式项目开发靠近.因此接下来一段时间,主要学习一下MVP架构.Retrofit以及RxJav ...

  6. Linux /etc/network/interfaces

    Linux下/etc/network/interfaces文件用来配置网络接口. 1. 使用动态IP地址 auto eth0 iface eth0 inet dhcp 2. 使用静态IP地址 auto ...

  7. JavaScript自学笔记(3)--- 用JS来实现网页浮窗

    最近做个小项目,给网页加个浮窗,考验了基础的css,js技术,还是蛮有意思的,代码如下(部分代码来源于引用,见底部) <!DOCTYPE html> <html> <he ...

  8. 安装Gitlab到CentOS(YUM)

    运行环境 系统版本:CentOS Linux release 7.3.1611 (Core) 软件版本:Gitlab-ce-11.10.1 硬件要求:最低2核4GB,建议4核8GB 安装过程 1.安装 ...

  9. C# Socket的方式发送中文,接收方收到乱码

    场景: 使用 Socket的方式向下位机发送中文信息,下位机收到的中文是乱码 原因: 了解到的原因是上位机与下位机的发送与接收的编码与解码方式不一致 比如上位机采用 Encoding.UTF8.Get ...

  10. Maven快速入门使用

    1. Maven 介绍 1.1 为什么使用 Maven 由于 Java 的生态非常丰富,无论你想实现什么功能,都能找到对应的工具类,这些工具类都是以 jar 包的形式出现的,例如 Spring,Spr ...