FZU - 2042 The Mad Mathematician 数位dp + 算贡献

传送门

文章目录

  • 题意:
  • 思路:

题意:

sumsumsum
在这里插入图片描述
a,b,c,d,e≤1e18a,b,c,d,e\le1e18a,b,c,d,e1e18

思路:

这是一篇无从考究的题解,因为fzu现在进不去。
看到这种题直接考虑数位dpdpdp,对于[A,B],[C,D][A,B],[C,D][A,B],[C,D]两个区间我们不太好直接处理,但是我们可以容斥一下,即[0,B][0,D]−[0,A−1][0,D]−[0,B][0,C−1]+[0,A−1][0,C−1][0,B][0,D]-[0,A-1][0,D]-[0,B][0,C-1]+[0,A-1][0,C-1][0,B][0,D][0,A1][0,D][0,B][0,C1]+[0,A1][0,C1]
假设我们求的[0,B][0,D][0,B][0,D][0,B][0,D],我们设计状态f[pos][flag1][flag2][flag3]f[pos][flag1][flag2][flag3]f[pos][flag1][flag2][flag3]
(1)flag1(1)flag1(1)flag1表示i<Bi<Bi<B
(2)flag2(2)flag2(2)flag2表示j<Dj<Dj<D
(3)flag3(3)flag3(3)flag3表示(ixorj)>E(i \ \ xor \ \ j) >E(i  xor  j)>E
到目前为止,如果这个题改成求方案数,就可以直接秒了,但是改不得。
dp1(pos)dp1(pos)dp1(pos)为从pospospos位开始的方案数,dp2(pos)dp2(pos)dp2(pos)为从pospospos位开始的贡献,考虑如何求贡献,我们直接按位来考虑,假设当前到了第pospospos位,且当前位(ixorj)==1(i \ \ xor \ \ j)==1(i  xor  j)==1,那么当前位的贡献就是(1ll<<pos)∗dp1(pos−1)+dp2(pos−1)(1ll<<pos)*dp1(pos-1)+dp2(pos-1)(1ll<<pos)dp1(pos1)+dp2(pos1),这个显然可以维护一个pairpairpair进行转移。
所以还是一个裸的数位dpdpdp
代码也比较好写辣,一定要注意取模,比如代码686868行,1ll<<pos1ll<<pos1ll<<pos一定要取模,不然炸LLLLLL了,因为这个对拍的时候一直过不了大样例。。。

// Problem: The Mad Mathematician
// Contest: Virtual Judge - FZU
// URL: https://vjudge.net/problem/FZU-2042
// Memory Limit: 32 MB
// Time Limit: 1000 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>
#include<random>
#include<cassert>
#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<LL,LL> PII;const int N=110,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;LL a,b,c,d,e;
int A[N],B[N],E[N];
PII f[70][2][2][2];PII dp(int pos,int flag1,int flag2,int flag3) {if(pos==-1) return {flag3,0};if(f[pos][flag1][flag2][flag3].X!=-1) return f[pos][flag1][flag2][flag3];int x=flag1? 1:A[pos];int y=flag2? 1:B[pos];int z=flag3? 0:E[pos];PII ans={0,0};for(int i=0;i<=x;i++) {for(int j=0;j<=y;j++) {if((i^j)<z) continue;PII now=dp(pos-1,flag1||i<x,flag2||j<y,flag3||((i^j)>z));ans.X+=now.X; ans.X%=mod;ans.Y+=now.Y;if((i^j)==1) ans.Y+=(1ll*(1ll<<pos)%mod)*now.X%mod;ans.Y%=mod;}}return f[pos][flag1][flag2][flag3]=ans;
}LL solve(LL x,LL y) {if(x<=0||y<=0) return 0;for(int i=62;i>=0;i--) {A[i]=x>>i&1;B[i]=y>>i&1;}memset(f,-1,sizeof(f));return dp(62,0,0,0).Y%mod;
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);//rd_wa();int _; scanf("%d",&_);for(int __=1;__<=_;__++) {scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&e);for(int i=62;i>=0;i--) E[i]=e>>i&1;LL ans=solve(b,d)%mod; ans-=solve(a-1,d); ans%=mod; ans+=mod; ans%=mod;ans-=solve(b,c-1); ans%=mod; ans+=mod; ans%=mod;ans+=solve(a-1,c-1); ans%=mod; ans+=mod; ans%=mod;printf("Case %d: %lld\n",__,ans%mod);}return 0;
}
/*
1
0 200000300000002 0 200000300002322 3
6860360
53360
*/

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

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

相关文章

「PowerBI」分析服务多维数据结构重回关系数据库的一大绝招

在过往Excel催化剂定位的轻量级ExcelPowerBIDesktop的解决方案中&#xff0c;已经做过了近乎完美的PowerBIDesktop数据模型数据导出到Excel工作表的应用&#xff0c;这也是个人版数据应用的重要的刚需功能&#xff0c;由Excel催化剂免费奉献给PowerBI社区作完结。此篇再次带领大…

【杭电多校2020】Fibonacci Sum【斐波拉契通项】【推式子】

题意&#xff1a;设FiF_iFi​为斐波拉契数列&#xff0c;求 ∑i0N(FiC)k\sum_{i0}^N(F_{iC})^ki0∑N​(FiC​)k 模109910^991099 N,C≤1018,k≤105N,C\leq10^{18},k\leq10^5N,C≤1018,k≤105 把斐波拉契暴力拆开 FiC15[(152)iC−(1−52)iC]F_{iC}\frac 1 {\sqrt{5}}[(\frac {…

Codeforces Round #257 (Div. 1) D. Jzzhu and Numbers 高维前缀和 + 容斥

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 完全想不到容斥啊&#xff0c;看了半天也没看懂渍渍渍。 定义f[i]f[i]f[i]表示iii的超集个数&#xff0c;那么选择的方案就是2f[i]−12^{f[i]}-12f[i]−1了&#xff0c;因为不能一个不选所…

花5分钟时间来了解一下高性能网关Kong会有意外收获

前言前几天开源发布了 Kong.Net 项目&#xff0c;收到了大量园友的反馈&#xff0c;开源当天就突破了 100 个star &#xff0c;可喜可贺&#xff0c;但是从侧面也说明&#xff0c;我们 .NetCore 阵营真的非常需要拥抱开源&#xff0c;应该敞开心扉&#xff0c;集众家之长&#…

【杭电多校2020】Total Eclipse【贪心】【并查集】

题意&#xff1a;nnn个点mmm条边的无向图&#xff0c;每个点有一个正点权&#xff0c;每次选择一个连通子图&#xff0c;将里面的权值都减111。求所有点权为000的最小步数。 T≤10,n≤105,m≤2105T\leq 10,n\leq 10^5,m\leq2\times10^5T≤10,n≤105,m≤2105 考虑一个贪心&…

Codeforces Round #725 (Div. 3) G. Gift Set 二分

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 有两种物品分别有x,yx,yx,y个&#xff0c;每次可以从一个拿出aaa个&#xff0c;另一个拿出bbb个分成一组&#xff0c;问最多能分成多少组。 思路&#xff1a; 这个题有一个显然的单调性&#xff0c;所以二…

长沙.NET技术社区·设计到实现

方案最终夭折回到一年前的四月初&#xff0c;也是这样一个气象更新&#xff0c;春意盎然的一年中最好的时光。 在长沙.NET技术圈中&#xff0c;有一个不知名的QQ群&#xff0c;叫做放肆.NET技术交流群&#xff08;群号&#xff1a;44329330&#xff09;&#xff0c;在这个群里面…

【杭电多校2020】Lead of Wisdom【搜索】【复杂度证明】

题意&#xff1a;nnn件物品&#xff0c;每种物品有一个种类tit_iti​&#xff0c;四个属性ai,bi,ci,dia_i,b_i,c_i,d_iai​,bi​,ci​,di​&#xff0c;每个种类最多选一件物品&#xff0c;求 (100∑ai)(100∑bi)(100∑ci)(100∑di)(100\sum a_i)(100\sum b_i)(100\sum c_i)(10…

Loj #6274. 数字 数位dp + 去重

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 考虑数位dpdpdp&#xff0c;设计状态f[pos][flag1][flag2][flag3][flag4]f[pos][flag1][flag2][flag3][flag4]f[pos][flag1][flag2][flag3][flag4]&#xff0c;其中flag1:x≥Lxflag2:y≤Rx…

微软发布 VS Code 容器化开发工具,大大简化物联网设备开发

近日&#xff08;北京时间 2019 年 7 月 9 日&#xff09;&#xff0c;微软宣布了 Azure IoT Tools for VS Code 的全新功能&#xff1a;通过在容器中封装编译器、设备 SDK 以及基本库&#xff0c;大大简化了设备交叉编译的工作量&#xff0c;使得嵌入式 Linux 设备&#xff08…

【杭电多校2020】Go Running【几何】【最小点覆盖】

题意&#xff1a;有一个数轴&#xff0c;若干人在某个时刻开始从某个点朝某个方向按111的速度走若干时间。已知nnn个条件&#xff0c;形如tit_iti​时刻xix_ixi​的位置有人&#xff0c;求最少可能的人数。 ∑n≤5105\sum n\leq 5\times 10^5∑n≤5105 画出 x−tx-tx−t二维平…

Codeforces Round #636 (Div. 3) F. Restore the Permutation by Sorted Segments 思维 + 暴力

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n≤200n\le200n≤200 思路&#xff1a; 首先关注到rrr从[2,n][2,n][2,n]都出现一次&#xff0c;所以很明显最后一个位置只出现一次&#xff0c;但是这样倒着来不是很好做考虑正着来。 我们可以枚举111位置填…

.NET中生成动态验证码

NET中生成动态验证码验证码是图片上写上几个字&#xff0c;然后对这几个字做特殊处理&#xff0c;如扭曲、旋转、修改文字位置&#xff0c;然后加入一些线条&#xff0c;或加入一些特殊效果&#xff0c;使这些在人类能正常识别的同时&#xff0c;机器却很难识别出来&#xff0c…

【NOI2019】弹跳【二维线段树】【dijkstra】

题意&#xff1a;一个whw\times hwh的二维平面上有nnn个城市&#xff0c;有mmm个弹跳装置&#xff0c;第iii个可以花费tit_iti​的时间从城市pip_ipi​跳到矩形x∈[l,r],y∈[u,d]x\in [l,r],y\in[u,d]x∈[l,r],y∈[u,d]中的任意一个城市。求从111到其他每个城市的最小时间。 w,…

Codeforces Beta Round #51 D. Beautiful numbers 数位dp + 状态优化

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 数位dpdpdp挺经典的一个题辣&#xff0c;有一个很明显的状态就是f[pos][num][lcm]f[pos][num][lcm]f[pos][num][lcm]表示到了第pospospos位&#xff0c;数是numnumnum&#xff0c;个位数最…

在C#中使用Irony实现SQL语句的解析

在上一篇博文 LogoSharp&#xff1a;Logo语言的C#实现中&#xff0c;我介绍了LOGO语言的C#实现&#xff0c;在该项目中&#xff0c;我使用了.NET下的语法解析框架&#xff1a;Irony。Irony框架最开始的时候是由Roman Ivantsov发布在codeplex上的一个开源项目&#xff0c;它使用…

【CF1045A】A Last chance【贪心】【线段树优化建图】【网络流构造方案】

题意&#xff1a;有nnn个武器和mmm个飞船&#xff0c;武器有下面三种 从给定的集合SSS中击破一个。在给定的区间[L,R][L,R][L,R]中击破一个。对于给定的a,b,ca,b,ca,b,c&#xff0c;选择000个或222个击破。特殊地&#xff0c;每个飞船最多被该操作的a,b,ca,b,ca,b,c指定一次。…

P4127 [AHOI2009]同类分布 数位dp + 对状态剪枝

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 比较套路的题&#xff0c;首先也有个明显的状态f[pos][num][sum]f[pos][num][sum]f[pos][num][sum]表示到了pospospos位&#xff0c;当前数为numnumnum&#xff0c;各位数字之和为sumsumsu…

.NET开发框架(六)-架构设计之IIS负载均衡(视频)

前面有关注我们公众号文章的朋友应该都知道&#xff0c;我们的分布式应用服务可以通过Ocelot网关进行负载均衡&#xff0c;这种方式属于应用级别的实现。而今天我们给大家介绍的是平台级别的实现&#xff0c;并且我们首次使用视频方式进行讲解。本文通过视频方式为大家演示&…

【NOI2019】斗主地【期望】【组合数学】【下降幂】【插值】

题意&#xff1a; 有nnn张牌&#xff0c;从上往下的第iii张牌权值是f(i)itypef(i)i^{type}f(i)itype。 进行mmm次洗牌操作&#xff0c;每次操作给定aia_iai​&#xff0c;将从上往下的前aia_iai​张牌分成一堆&#xff0c;后面bib_ibi​张牌分成第二堆。如果两堆的牌数分别为…