CF296B dp\容斥

传送门

文章目录

  • 题意:
  • 思路:

题意:

在这里插入图片描述
n≤1e5n\le1e5n1e5

思路:

求方案数基本就是考虑dpdpdp了,看到nnn这么大可以考虑一下分情况讨论的dpdpdp状态。
f[i][j]f[i][j]f[i][j]表示到了第iii个,状态为jjj的方案数,考虑如何定义jjj才能使方案不重不漏。
(1)j==0(1)j==0(1)j==0表示前面只有sk=wks_k=w_ksk=wk的方案数。
(2)j==1(2)j==1(2)j==1表示前面只有sj<wjs_j<w_jsj<wj但是可以出现sk=wks_k=w_ksk=wk的方案数。
(3)j==2(3)j==2(3)j==2表示前面只有sj>wjs_j>w_jsj>wj但是可以出现sk=wks_k=w_ksk=wk的方案数。
(4)j==3(4)j==3(4)j==3表示前面有sj<wj,sj>wjs_j<w_j,s_j>w_jsj<wj,sj>wj但是可以出现sk=wks_k=w_ksk=wk的方案数。
初始状态f[0][0]=1f[0][0]=1f[0][0]=1,让后递推即可,转移方程:
(1)ai>bi(1)a_i>b_i(1)ai>bif[i][2]+=f[i−1][2]+f[i−1][0]f[i][2]+=f[i-1][2]+f[i-1][0]f[i][2]+=f[i1][2]+f[i1][0]f[i][3]+=f[i−1][3]+f[i−1][1]f[i][3]+=f[i-1][3]+f[i-1][1]f[i][3]+=f[i1][3]+f[i1][1]
(2)ai=bi(2)a_i=b_i(2)ai=bif[i][0]+=f[i−1][0]f[i][0]+=f[i-1][0]f[i][0]+=f[i1][0]f[i][1]+=f[i−1][1]f[i][1]+=f[i-1][1]f[i][1]+=f[i1][1]f[i][2]+=f[i−1][2]f[i][2]+=f[i-1][2]f[i][2]+=f[i1][2]f[i][3]+=f[i−1][3]f[i][3]+=f[i-1][3]f[i][3]+=f[i1][3]
(3)ai<bi(3)a_i<b_i(3)ai<bif[i][1]+=f[i−1][1]+f[i−1][0]f[i][1]+=f[i-1][1]+f[i-1][0]f[i][1]+=f[i1][1]+f[i1][0]f[i][3]+=f[i−1][3]+f[i−1][2]f[i][3]+=f[i-1][3]+f[i-1][2]f[i][3]+=f[i1][3]+f[i1][2]
解释一下转移的含义:
(1)(1)(1)出现了ai>bia_i>b_iai>bi的情况,我们就可以把ai<bia_i<b_iai<bi的情况加到答案中了,让后再让ai=bia_i=b_iai=bi的情况加到ai>bia_i>b_iai>bi的情况中即可。之后这两个状态就变为无效状态了。
(2)(2)(2)出现了ai=bia_i=b_iai=bi的情况,那么对所有情况递推过来即可。
(3)(3)(3)出现了ai<bia_i<b_iai<bi,与上面的(1)(1)(1)情况类似,也是去掉了两个状态。
对于???我们没法比较大小,但是我们直接枚举[0,9][0,9][0,9]转移就好啦。

当然还可以容斥来做,无非就四种情况:
(1)(1)(1)全部情况ansansans
(2)(2)(2)全部大于的情况ans1ans1ans1
(3)(3)(3)全部等于的情况ans2ans2ans2
(4)(4)(4)全部小于的情况ans3ans3ans3
答案即为ans−ans1−ans2−ans3ans-ans1-ans2-ans3ansans1ans2ans3
也比较容易算,最后给出了代码。

dpdpdp

// Problem: B. Yaroslav and Two Strings
// Contest: Codeforces - Codeforces Round #179 (Div. 2)
// URL: https://codeforces.com/problemset/problem/296/B
// Memory Limit: 256 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n;
char a[N],b[N];
LL f[N][5];
//0 前面只有s_j=w_j
//1 前面只有s_j<w_j
//2 前面只有s_j>w_j
//3 前面有 s_j<w_j,s_j>w_jvoid cal(int j) {for(int i=0;i<4;i++) f[j][i]%=mod;
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);cin>>n>>(a+1)>>(b+1);f[0][0]=1;for(int i=1;i<=n;i++) {if(a[i]!='?'&&b[i]!='?') {if(a[i]>b[i]) {f[i][2]+=f[i-1][2]+f[i-1][0];f[i][3]+=f[i-1][3]+f[i-1][1];}else if(a[i]==b[i]) {f[i][0]+=f[i-1][0];f[i][1]+=f[i-1][1];f[i][2]+=f[i-1][2];f[i][3]+=f[i-1][3];}else {f[i][1]+=f[i-1][1]+f[i-1][0];f[i][3]+=f[i-1][3]+f[i-1][2];}cal(i);}else if(a[i]=='?'&&b[i]!='?') {for(int j='0';j<='9';j++) {if(j>b[i]) {f[i][2]+=f[i-1][2]+f[i-1][0];f[i][3]+=f[i-1][3]+f[i-1][1];}else if(j==b[i]) {f[i][0]+=f[i-1][0];f[i][1]+=f[i-1][1];f[i][2]+=f[i-1][2];f[i][3]+=f[i-1][3];}else {f[i][1]+=f[i-1][1]+f[i-1][0];f[i][3]+=f[i-1][3]+f[i-1][2];}cal(i);}}else if(a[i]!='?'&&b[i]=='?') {for(int j='0';j<='9';j++) {if(a[i]>j) {f[i][2]+=f[i-1][2]+f[i-1][0];f[i][3]+=f[i-1][3]+f[i-1][1];}else if(a[i]==j) {f[i][0]+=f[i-1][0];f[i][1]+=f[i-1][1];f[i][2]+=f[i-1][2];f[i][3]+=f[i-1][3];}else {f[i][1]+=f[i-1][1]+f[i-1][0];f[i][3]+=f[i-1][3]+f[i-1][2];}cal(i);}}else {for(int j='0';j<='9';j++) {for(int k='0';k<='9';k++) {if(j>k) {f[i][2]+=f[i-1][2]+f[i-1][0];f[i][3]+=f[i-1][3]+f[i-1][1];}else if(j==k) {f[i][0]+=f[i-1][0];f[i][1]+=f[i-1][1];f[i][2]+=f[i-1][2];f[i][3]+=f[i-1][3];}else {f[i][1]+=f[i-1][1]+f[i-1][0];f[i][3]+=f[i-1][3]+f[i-1][2];}}cal(i);}}cal(i);}for(int i=1;i<=n;i++) {for(int j=0;j<=3;j++) printf("%d ",f[i][j]);puts("");}printf("%lld\n",f[n][3]%mod);return 0;
}
/**/

容斥

// Problem: B. Yaroslav and Two Strings
// Contest: Codeforces - Codeforces Round #179 (Div. 2)
// URL: https://codeforces.com/problemset/problem/296/B
// Memory Limit: 256 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n;
char a[N],b[N];LL getall() {LL ans=1;for(int i=1;i<=n;i++) {if(a[i]=='?') ans*=10;if(b[i]=='?') ans*=10;ans%=mod;}return ans%mod;
}LL get1() {LL ans=1;for(int i=1;i<=n;i++) {if(a[i]=='?'&&b[i]=='?') ans*=55;else if(a[i]=='?') ans*=(b[i]-'0'+1);else if(b[i]=='?') ans*=(10-(a[i]-'0'));else if(a[i]>b[i]) return 0;ans%=mod;}return ans%mod;
}LL get2() {LL ans=1;for(int i=1;i<=n;i++) {if(a[i]=='?'&&b[i]=='?') ans*=10;else if(a[i]=='?'||b[i]=='?') continue;else if(a[i]!=b[i]) return 0ll;ans%=mod;}return ans%mod;
}LL get3() {LL ans=1;for(int i=1;i<=n;i++) {if(a[i]=='?'&&b[i]=='?') ans*=55;else if(b[i]=='?') ans*=(a[i]-'0'+1);else if(a[i]=='?') ans*=(10-(b[i]-'0'));else if(a[i]<b[i]) return 0;ans%=mod;}return ans%mod;
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);/*int cnt=0;for(int i=0;i<=9;i++) for(int j=0;j<=9;j++) if(i>=j) cnt++;cout<<cnt<<endl;*/cin>>n>>(a+1)>>(b+1);LL ans=getall();ans+=get2(); ans%=mod;ans-=get1(); ans+=mod; ans%=mod;ans-=get3(); ans+=mod; ans%=mod;printf("%lld\n",ans%mod);return 0;
}
/**/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315107.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

asp.net core 自定义异常处理中间件

Intro在 asp.net core 中全局异常处理&#xff0c;有时候可能不能满足我们的需要&#xff0c;可能就需要自己自定义一个中间件处理了&#xff0c;最近遇到一个问题&#xff0c;有一些异常&#xff0c;不希望记录错误日志&#xff0c;目前主要是用户请求取消导致的 TaskCanceled…

CF786B Legacy 线段树优化建图

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 实现如下连边后跑最短路。 思路&#xff1a; 优化建图板子题&#xff0c;优化思路就是将区间分割成若干个线段树上的线段&#xff0c;与线段树分治有点类似&#xff0c;由于有点向区间也有区间向点的边&a…

【ZJOI2015】幻想乡 Wi-Fi 搭建计划【几何】【贪心】【dp】

传送门 题意&#xff1a;一个x∈(−∞,∞),y∈[0,R]x\in(-\infin,\infin),y\in[0,R]x∈(−∞,∞),y∈[0,R]的矩形中有nnn个点&#xff0c;矩形外有mmm个半径均为RRR的圆&#xff0c;有独立的代价cic_ici​。求覆盖最多的点所需的最小代价。 n,m≤100n,m\leq100n,m≤100 显然先…

.NET架构开发应知应会

.NET程序是基于.NET framework、.NET Core、Mono、UWP【.NET实现】开发和运行的 &#xff0c;定义以上【.NET实现】的标准规范称为.NET StandardL1&#xff1a;.NET Standard.NET标准是一组API集合&#xff0c;由上层三种【.NET实现】的Basic Class Library实现&#xff0c;更正…

几个冷门字符串算法的学习笔记(最小表示法,exKMP,Lyndon Word)

所有下标均从1开始 最小表示法 给定一个串&#xff0c;求字典序最小的循环同构。 我们把串复制一遍接在后面&#xff0c;然后求出[1,N][1,N][1,N]开始的长为NNN的子串中最小的 先设i1,j2i1,j2i1,j2 然后暴力找出iii和jjj往后匹配的第一个不同的位置&#xff0c;记为ikikik…

.NET Core IdentityServer4实战 第Ⅴ章-单点登录

OiDc可以说是OAuth的改造版&#xff0c;在最初的OAuth中&#xff0c;我们需要先请求一下认证服务器获取下Access_token&#xff0c;然后根据Access_token去Get资源服务器, 况且OAuth1 和 2 完全不兼容&#xff0c;易用性差&#xff0c;而OIDC可以在登陆的时候就把信息返回给你&…

【CF594E】Cutting the Line 【贪心】【Lyndon Word】【扩展kmp】

传送门 题意&#xff1a;给一个字符串SSS和正整数kkk&#xff0c;将SSS分成最多kkk段&#xff0c;每段不变或翻转&#xff0c;使得最后的字典序最小。 ∣S∣≤5106|S|\leq5\times10^6∣S∣≤5106 发现不翻转可以看成拆成若干单字符分别翻转&#xff0c;所以先分析一下必须翻转…

一份好的工作总结才能帮你升职加薪

这里是Z哥的个人公众号每周五早8点 按时送达当然了&#xff0c;也会时不时加个餐&#xff5e;我的第「79」篇原创敬上最近有点忙&#xff0c;搬出之前攒的一篇文章来应急一下。一篇能助你挣更多钱的文章。好了&#xff0c;下面开始。我的读者们大部分是互联网行业的&#xff0c…

腾讯开源软件镜像站上线

腾讯开源软件镜像站(Tencent Open Source Mirror Site)已于近日上线&#xff0c;其官方名称为「腾讯云软件源」&#xff0c;由腾讯云提供支持。官方表示搭建此开源镜像站的目的在于宣传自由软件的价值&#xff0c;提高自由软件社区文化氛围&#xff0c;推广自由软件在国内的应用…

ASP.NET Core on K8S学习初探(2)

“ [LOG] ASP.NET Core on K8S Starting...”在上一篇《单节点环境搭建》中&#xff0c;通过Docker for Windows在Windows开发机中搭建了一个单节点的K8S环境&#xff0c;接下来就是动人心弦的部署ASP.NET Core API到K8S了。但是&#xff0c;在部署之前&#xff0c;我还是把基本…

Educational Codeforces Round 96 E. String Reversa 线段树模拟序列交换

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 与上一篇题解大同小异&#xff0c;无非就是不需要枚举排列了。 // Problem: E. String Reversal // Contest: Codeforces - Educational Codeforces Round 96 (Rated for Div. 2) // URL:…

ASP.NET Core on K8S学习初探(1)

“ [LOG] ASP.NET Core on K8S Starting...”01—写在之前当近期的一个App上线后&#xff0c;发现目前的docker实例&#xff08;应用服务BFF中台服务工具服务&#xff09;已经越来越多了&#xff0c;而我司目前没有专业的运维人员&#xff0c;发现运维的成本逐渐开始上来&#…

.NET Core IdentityServer4实战 第Ⅳ章-集成密码登陆模式

回顾下ClientCredentials模式&#xff0c;在ReSourceApi中定义了我们公开服务&#xff0c;第三方网站想要去访问ReSourceApi则需要在身份验证服务中获取toekn&#xff0c;根据token的内容&#xff0c;硬编码去访问公开服务&#xff08;ResApi&#xff09;,这个还是非常简单的&a…

Codeforces Round #656 (Div. 3) F. Removing Leaves 贪心 + 模拟

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先有一个贪心策略就是每次都找一个叶子节点最多的点&#xff0c;让后删掉他的kkk个叶子节点&#xff0c;现在我们就来考虑如何模拟这个过程。 我们整一个vector<set<int>>ve…

.NetCore中三种注入生命周期的思考

.NetCore彻底诠释了“万物皆可注入”这句话的含义&#xff0c;在.NetCore中到处可见注入的使用。因此core中也提供了三种注入方式的生命周期使用&#xff0c;分别是&#xff1a;AddTransient&#xff1a;每次请求&#xff0c;都获取一个新的实例。即使同一个请求获取多次也会是…

兰州大学第一届 飞马杯 体育课排队 二分 + 最大流 + 输出路径

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 读懂题就会辣&#xff0c;经典模型了&#xff0c;二分时间&#xff0c;让后将其转换成二分图&#xff0c;左边是人&#xff0c;右边是位置&#xff0c;能在规定时间到的连边&#xff0c;跑…

从零开始实现ASP.NET Core MVC的插件式开发(一) - 使用Application Part动态加载控制器和视图...

如果你使用过一些开源CMS的话&#xff0c;肯定会用过其中的的插件化功能&#xff0c;用户可以通过启用或者上传插件包的方式动态添加一些功能&#xff0c;那么在ASP.NET Core MVC中如何实现插件化开发呢&#xff0c;下面我们来探究一下。本系列只是笔者的一些尝试&#xff0c;并…

兰州大学第一届 飞马杯 ★★快乐苹果树★★ 树链剖分 + 懒标记 + 树状数组

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 第一次听说树链剖分能在fa[top[i]]fa[top[i]]fa[top[i]]的地方加懒标记&#xff0c;学到了学到了。 首先不能被题目吓住&#xff0c;这个题目仔细剖析一下不难发现一些性质&#xff1a; 以…

.NET分布式框架 | Orleans 知多少

引言公司物联网项目集成Orleans以支持高并发的分布式业务&#xff0c;对于Orleans也是第一次接触&#xff0c;本文就分享下个人对Orleans的理解。这里先抛出自己的观点&#xff1a;Orleans 是一个支持有状态云生应用/服务水平伸缩的基于Virtual Actor 模型的.NET分布式框架。下…

nowcoder 牛牛的最大兴趣组 质因子 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先nnn很小的话可以暴力连边&#xff0c;让后染个色求一个颜色最多的即可。但是这个题显然不行&#xff0c;由于是三次方&#xff0c;所以考虑质因子入手。 首先很容易就能想到将所有的数…