黑客(续) (压位高精+状压dp)

黑客(续)

  • description
  • solution
  • code

description

【问题描述】
在破解了世界首富 Bychaha 的银行账户后,知名黑客 pks 发现,要得到
Bychaha 的全部财产,必须再破解一道密码。

作为客户账户安全的最后一道防线,这一次的密码将由长达 N 位的数码组
成,每一位的数码范围为[1,K]。

pks 想要估算自己破解密码的大致时间,所以他想要你帮他快速计算出,总
共有多少种满足条件的密码,同时,pks 还对每一种密码视为十进制数之后求
和的结果很感兴趣,希望你也能告诉他。

当然,作为世界知名黑客,pks 不会傻傻的枚举,他已经从银行系统中窃取
到了关于密码的 M 个信息,每个信息由两个数字𝑎,𝑏表示,代表数码𝑎不会出
现在数码𝑏之前。

【输入格式】
第一行三个整数 N,M,K,含义如题面所述。
接下来 M 行,每行两个数𝑎,𝑏。

【输出格式】
输出共两行。
第一行一个整数表示有多少种满足条件的密码。
第二行一个整数表示所有密码视为十进制数之后求和的结果。

【样例输入输出】

4 4 3 
1 1 
1 2 
2 2 
3 1 
7 
19020 

【数据规模与约定】
对于 20%的数据, 𝑁 ≤ 6
对于 30%的数据, 𝑁 ≤ 50
对于 40%的数据, 𝑁 ≤ 200
另有 50%的数据,𝑀 = 0
对于 100%的数据,1 ≤ 𝑁 ≤ 500,0 ≤ 𝑀 ≤ 100 , K ≤ 9

solution

当你发现没有模数的那一刻,可能心脏已经跳脱了

这该死的 甜美 大数还是来了

N<=6

直接暴力搜索,每一位选择[1,k][1,k][1,k]最后组合出来,再判断


m=0

没有限制条件,摆明了说个数就是knk^nkn,答案就是每一个数字固定在每一位的数值乘以这种形式的密码个数

直接脱光了说:答案都跟你说了,你写出大数就有一半的分,再加上最原始的暴力,就算是很不错的分了

没想到自己在考场上真的把大数敲出来了


最后就是正解了

很简单,可以预处理出已经出现了sss集合内的数字,且是合法的,能接在后面的所有数字

fi,s/gi,s:f_{i,s}/g_{i,s}:fi,s/gi,s: 在第iii位已经出现了sss集合内的数字的所有密码和/个数

这个状压转移应该都能写吧

但是会发现,答案特别大,就需要——压位高精

这里好巧不巧选择压171717位(1e171e171e17)就能过

之前以为压位很难,结果发现就是将是十七位当成一位,(mod1e17)\pmod {1e17}(mod1e17)就行

最后输出再每一位扩充出来%17lld\text{\%17lld}%17lld171717是表示输出是171717位不够的自动填充000

当然注意最高的一位不能扩充,ta是多少位就是多少位


所以这道题的难点是压位高精??

code

#pragma GCC optimize(2)
#include <cstdio>
#include <vector>
#include <cstring>
using namespace std;
#define int long long
const int mod = 1e17;
int n, m, k;
bool vis[10][10];
vector < int > nxt[1 << 9];struct Int {int len;int c[60];Int(){ memset( c, 0, sizeof( c ) ); }Int operator * ( int x ) {Int ans;ans.len = len;int add = 0;for( int i = 1;i <= len;i ++ ) {int ret = c[i] * x + add;ans.c[i] = ret % mod;add = ret / mod;}if( add ) ans.c[++ ans.len] = add;return ans;}Int operator + ( Int x ) {Int ans;int ip = 1;while( ip <= x.len or ip <= len ) {int ret = ans.c[ip] + x.c[ip] + c[ip];if( ret >= mod ) {ans.c[ip + 1] ++;ans.c[ip] = ret - mod;}elseans.c[ip] = ret;++ ip;}ans.len = ip;if( ! ans.c[ans.len] ) ans.len --;return ans;}void print() {printf( "%lld", c[len] );for( int i = len - 1;i > 0;i -- ) printf( "%017lld", c[i] );printf( "\n" );}void init() {for( int i = 1;i <= len;i ++ ) c[i] = 0;len = 0;}
}f[2][1 << 9], g[2][1 << 9];signed main() {scanf( "%lld %lld %lld", &n, &m, &k );for( int i = 1, a, b;i <= m;i ++ ) {scanf( "%lld %lld", &a, &b );vis[a][b] = 1;}int lim = 1 << k;for( int s = 0;s < lim;s ++ )for( int i = 1;i <= k;i ++ ) {for( int j = 1;j <= k;j ++ )if( ( 1 << j - 1 & s ) and vis[j][i] ) goto opt;nxt[s].push_back( i );opt :;}g[0][0].c[g[0][0].len = 1] = 1;for( int i = 1;i <= n;i ++ ) {int o = i & 1;for( int s = 0;s < lim;s ++ )for( int j : nxt[s] ) {f[o][1 << j - 1 | s] = f[o][1 << j - 1 | s] + f[o ^ 1][s] * 10 + g[o ^ 1][s] * j;g[o][1 << j - 1 | s] = g[o][1 << j - 1 | s] + g[o ^ 1][s];}for( int s = 0;s < lim;s ++ )f[o ^ 1][s].init(), g[o ^ 1][s].init();}g[n & 1][0].c[1] = 0;for( int i = 1;i < lim;i ++ ) {f[n & 1][0] = f[n & 1][i] + f[n & 1][0];g[n & 1][0] = g[n & 1][i] + g[n & 1][0];}g[n & 1][0].print();f[n & 1][0].print();return 0;
}

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

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

相关文章

东莞.NET俱乐部线下技术沙龙-活动报名

自广州.NET技术俱乐部在2018年12月08日线下活动顺利开展后&#xff0c;东莞作为兄弟城市&#xff0c;也想通过线下活动的方式&#xff0c;点燃东莞.NET技术的熊熊之火。现决定先借助广州、深圳兄弟城市的帮助下&#xff0c;开展一场东莞方主办的线下活动&#xff0c;聚集东莞本…

Acwing 1084. 数字游戏 II

Acwing 1084. 数字游戏 II 题意&#xff1a; 指定一个整数闭区间 [a.b]&#xff0c;问这个区间内有多少个取模数。 取模数&#xff1a;这种数字必须满足各位数字之和 mod N 为 0。 题解&#xff1a; 数位dp 这里不细讲数位dp了&#xff0c;可以看看 Acwing 1081. 度的数量&…

ybtoj洛谷P3268:圆的异或并(扫描线)

解析 很神奇的一道题。 关键条件&#xff1a;任意两个圆无交。 把一个圆分成上下两个圆弧&#xff0c;那么所有圆弧的高度关系不会发生变化。 所以可以开一个 set&#xff0c;维护一个从左往右扫的扫描线&#xff0c;按照当前扫描线的横坐标定义比较符号&#xff0c;在圆的最…

如何撰写较受欢迎的技术文章

本来我这篇文章的标题是 “如何撰写受欢迎的技术文章”&#xff0c;但反复斟酌之下&#xff0c;还是加了一个“较”字&#xff0c;这主要是考虑我不是什么知名作者&#xff0c;写的文章大多也谈不上很受欢迎&#xff0c;贸然地谈“受欢迎” 是有点忐忑的&#xff0c;而改成现在…

朝鲜时蔬(分数据点写算法+毒瘤数学)

朝鲜时蔬decriptionsolutioncodedecription 对于一个有穷非空正整数集合S{x1,x2,x3,...,xn}⊂N(n≥1)S\{x_1,x_2,x_3,...,x_n\}\subset N^(n\ge 1)S{x1​,x2​,x3​,...,xn​}⊂N(n≥1)&#xff0c;定义其和sum(S)sum(S)sum(S)为所有元素的和 sum(S)x1...xnsum(S)x_1...x_nsum…

洛谷P7515:矩阵游戏(差分约束)

解析 如果没有元素均要在 [0,1e6][0,1e6][0,1e6] 的条件&#xff0c;可以很容易的构造出一个合法解。 那么我们就要通过调整得到的解&#xff0c;使所有数都在合法范围内。 注意到&#xff0c;每次给某一行/列依次1,-1,1,-1…这样仍然符合要求。 让每一行/列错开&#xff0c;…

Recursive sequence HDU - 5950

Recursive sequence HDU - 5950 题意&#xff1a; 给你一个式子&#xff1a;f[n]2f[n-2]f[n-1]n4 给你f[1]和f[2]&#xff0c;给你一个n&#xff0c;求f[n] f[1],f[2],n<231 题解&#xff1a; 很明显&#xff0c;矩阵快速幂&#xff0c;但是太久没做这种题&#xff0c;我…

使用Http-Repl工具测试ASP.NET Core 2.2中的Web Api项目

今天&#xff0c;Visual Studio中没有内置工具来测试WEB API。使用浏览器&#xff0c;只能测试http GET请求。您需要使用Postman&#xff0c;SoapUI&#xff0c;Fiddler或Swagger等第三方工具来执行WEB API的完整测试。在ASP.NET Core 2.2中&#xff0c;引入了一个名为“http-r…

洛谷P7518:宝石(倍增、可撤销并查集)

解析 算法一 定义 upx,kup_{x,k}upx,k​ 为节点 xxx 从自己的颜色所在位置在返祖链上往后跳 2k2^k2k 个颜色到达的节点。 可以像倍增一样的求解。 这样对于一次询问 (s,t)(s,t)(s,t) 我们就能求出 (s,lca)(s,lca)(s,lca) 这一段能取到哪里了。 对于向下的情况&#xff0c;再处…

E - Counting Cliques HDU - 5952

E - Counting Cliques HDU - 5952 题意&#xff1a; 给你n个点&#xff0c;m个边&#xff0c;还有一个s&#xff0c;问这个图中有多少个等于s的点集可以组成一个完全图 题解&#xff1a; 这题。。直接暴力搜索就行 分析复杂度的时候&#xff0c;应该考虑只有1000条边&#…

Docker最全教程之使用TeamCity来完成内部CI、CD流程(十七)

本篇教程主要讲解基于容器服务搭建TeamCity服务&#xff0c;并且完成内部项目的CI流程配置。教程中也分享了一个简单的CI、CD流程&#xff0c;仅作探讨。不过由于篇幅有限&#xff0c;完整的DevOps&#xff0c;我们后续独立探讨。 为了降低容器的使用门槛以及便于大家将容器技…

2021牛客NOIP提高组第二场T2——方格计数(组合数计数)

方格计数descriptionsolutioncodedescription 在左下角是 (&#x1d7ce;, &#x1d7ce;)&#xff0c;右上角是 (W, H)的网格上&#xff0c;有 (W 1) (H 1) 个格点。 现在要在格点上找 N个不同的点&#xff0c;使得这些点在一条直线上。并且在这条直线上&#xff0c; 相邻…

Acwing 1085. 不要62

Acwing 1085. 不要62 题意&#xff1a; 问[n,m]这些数中有多少数不包含4&#xff0c;且不包含连续的62 题解&#xff1a; 经典数位dp&#xff0c;分析过程以前的数位dp博客有写 Acwing 1082. 数字游戏 代码&#xff1a; #include<bits/stdc.h> #define debug(a,b) …

ybtoj洛谷P4406三角形面积并(扫描线)

解析 暴力求出所有三角形之间的所有交点&#xff0c;提出所有的横坐标。 然后任意两个相邻的横坐标之间的面积都是若干个梯形。 那么就可以求出对于每一个横坐标截得的三角形长度的并的和&#xff0c;然后加在一起乘高除以二即可。 在这里插入代码片#include<bits/stdc.h&…

自动将 NuGet 包的引用方式从 packages.config 升级为 PackageReference

在前段时间我写了一篇迁移 csproj 格式的博客 将 WPF、UWP 以及其他各种类型的旧 csproj 迁移成基于 Microsoft.NET.Sdk 的新 csproj&#xff0c;不过全过程是手工进行的&#xff0c;而且到最后处理 XAML 问题也非常头疼。现在&#xff0c;我们可以利用工具自动地完成这个过程。…

[NowCoder牛客]2021NOIP提高组模拟赛第二场T3——树数树(启发式合并堆)

树数树descriptionsolutioncodedescription 【题目描述】 牛牛有一棵 n 个点的有根树&#xff0c;根为 1。 我们称一个长度为 m 的序列 a 是好的&#xff0c;当且仅当&#xff1a; • ∀&#x1d456;∈(1,&#x1d45a;]∀&#x1d456;∈(1, &#x1d45a;]∀i∈(1,m]&#…

Acwing1086. 恨7不成妻(未解决)

Acwing1086. 恨7不成妻 题意&#xff1a; 问一个区间内与7无关的整数的平方和 与7有关包括&#xff1a; 1.整数中某一位是 7&#xff1b; 2.整数的每一位加起来的和是 7 的整数倍&#xff1b; 3.这个整数是 7 的整数倍。 题解&#xff1a; 数位dp&#xff0c;有点难。。 这…

模板:拓展kmp(Z函数)

所谓拓展kmp&#xff0c;就是拓展的kmp &#xff08;逃&#xff09; 前言 小清新算法&#xff0c;更像一个trick。 exkmp能够在线性复杂度内求出字符串所有后缀与字符串本身的lcp。 个人感觉这个东西和kmp关系不大&#xff0c;反而有些像马拉车的思想。 之前学的东西&#xf…

牛客NOIP2021提高组OI赛前模拟赛第一场T3——与巨(数学)

与巨descriptionsolutioncodedescription 【题目描述】 定义无穷序列f:f11,fnfn−1∗21f:f_11,f_nf_{n-1}*21f:f1​1,fn​fn−1​∗21 定义函数G(x)min⁡fi≥x(fi)G(x)\min_{f_i\ge x}(f_i)G(x)minfi​≥x​(fi​) 定义dpc,00,dpc,imax⁡(dpc,i−1,[((i∗c)&G(i))i]∗i)…

ASP.NET Core 自定义认证方式--请求头认证

Intro最近开始真正的实践了一些网关的东西&#xff0c;最近写几篇文章分享一下我的实践以及遇到的问题。本文主要介绍网关后面的服务如何进行认证。解决思路网关可以做一部分的认证和授权&#xff0c;服务内部有时候也会需要用户的信息&#xff0c;这时该怎么办呢&#xff0c;我…