打了计蒜客上的银川重现赛,总体感觉难度上确实比平时区域赛要低上一些。

这里补一下F题和G题的思路和代码。

upd:I题也补了,理解差不多都在注释里。

  • F题

       做法,玩一下n=10的样例就出来啦!

  解释:显然a^x的反函数为logax,我们先固定外层的求和的a,然后看内层求和的b,b从a开始加到n,注意到对于后半个向上取整的logba,b>=a,所以始终都是1,而对于前半个式子,只有当b经过a^i时才增加,举个例子就是log22~log23向下取整都为1,log24~log27向下取整都为2,log28~log210都为3.

对于n=10的样例,a=2就可以这样玩出来,在玩a=3,同理,b从3~10的过程中,9和10的时候出来2。

而对于4和4以后的值,我们发现后面求和的值都是n-i+1个1,就是从i到n中有多少个数,这部分我们就可以用公式算。咋算呢?

Σi(n-i+1),显然 i 从sqrt(n)+1到n,把这个式子拆开来,就可以得到公式,此处注意 12+22+……+n2=n*(n+1)*(2n+1)/6;

至此,我们就可以发现规律了!n最大到1e12,我们就可以对前sqrt(n)个数暴力的去算。后面部分的完美地去用公式算。

不过比赛的时候,取模写炸了,导致10几分钟就推出来的结论,爆炸搞了2个小时,取模要每步取模,由于这里的n本身就很大,每个数本身都要取模,而更重要的是对于除法,取模意义下要用逆元来算,所以就用到了快速幂来求除法逆元。主要是这里除法经常炸,亏我对拍debug了半天。

代码:

 #include <bits/stdc++.h>
#define debug(x) cout << #x << ": " << x << endl
using namespace std;
typedef long long ll;
const int MAXN=2e5+;
const int INF=0x3f3f3f3f;
const int MOD=; ll quick(ll x,ll n) //快速幂 x^n
{
ll res=;
while(n)
{
if(n&) res=(res*x)%MOD;
x=(x*x)%MOD;
n>>=;
}
return res;
} ll inv(ll a) //逆元 费马小定理,要求 a,mod互素
{
return quick(a,MOD-);
} ll i2sum(ll i)
{
return (((i%MOD)*((i+)%MOD))%MOD*(((*i)%MOD+)%MOD))%MOD*inv()%MOD;
} int main()
{
ios::sync_with_stdio(false);
cin.tie();
ll n;
cin>>n;
ll ans=;
ll k=(ll)sqrt(n)+;
//debug(k);
for(ll i=;i<k;++i)
{
ll t=i;
ll q=;
while(t<=n)
{
q=(q+n-t+)%MOD;
t*=i;
}
ans=(ans+(q%MOD)*(i%MOD)%MOD)%MOD;
//debug(ans);
}
//debug(ans);
ans=(ans+(((((n+)%MOD)*((n-k+)%MOD))%MOD*(((k+n)%MOD)*inv()%MOD)%MOD)%MOD+i2sum(k-)-i2sum(n)+MOD)%MOD)%MOD;
// debug(ans);
//ll t=i2sum(n);
//debug(t);
cout<<ans<<endl;
return ;
}
//200 1366404
//1000 167464908
//10000 36893337
//100000 384927495
//1000000 960529847
//10000000 679483439
//100000000 498142384

debug和对拍数据没删,也可以对照一下。

  • G题

g题一开始看题目看傻了,p是啥完全不知道,后来和队长讨论了一下,结合样例一看,就对上感觉了。

因为操作就是乘2到10之间的一个数,很显然pm | n 且  不pm+1 | n  的m 就是你已经对p乘的操作,需要对照样例感觉一下。

然后维护四棵 区间加法,维护最大值的线段树了。

然后对于2~10的数你可以写if else ,也可以像我这样写,每次差几个log级update的操作问题不大。

然后这两天的数据结构和计算几何题,告诉我这类题一定要scanf,printf,开优化的cin也不行!!!

然后就没啦!(cry!

 #include <bits/stdc++.h>
#define debug(x) cout << #x << ": " << x << endl
#define lson (rt<<1)
#define rson (rt<<1|1) using namespace std;
typedef long long ll;
const int MAXN=2e5+;
const int INF=0x3f3f3f3f;
const int MOD=1e9+;
int prime[]={,,,}; struct tree
{
int seg[MAXN<<];
int lazy[MAXN<<];
inline void pushup(int rt){ seg[rt]=max(seg[lson],seg[rson]); }
inline void pushdown(int rt)
{
if(lazy[rt])
{
lazy[lson]+=lazy[rt];
lazy[rson]+=lazy[rt];
seg[lson]+=lazy[rt];
seg[rson]+=lazy[rt];
lazy[rt]=;
}
}
void build(int l,int r,int rt)
{
seg[rt]=lazy[rt]=;
if(l==r) {seg[rt]=;return;}
int m=(l+r)>>;
build(l,m,lson);
build(m+,r,rson);
pushup(rt);
} void update(int l,int r,int rt,int ql,int qr,ll x)
{
if(ql<=l && r<=qr) {seg[rt]+=x;lazy[rt]+=x;return;}
pushdown(rt);
int m=(l+r)>>;
if(ql<=m) update(l,m,lson,ql,qr,x);
if(qr>m) update(m+,r,rson,ql,qr,x);
pushup(rt);
} int query(int l,int r,int rt,int ql,int qr)
{
if(ql<=l && r<=qr) return seg[rt];
pushdown(rt);
int m=(l+r)>>;
int ans=-INF;
if(ql<=m) ans=max(ans,query(l,m,lson,ql,qr));
if(qr>m) ans=max(ans,query(m+,r,rson,ql,qr));
return ans;
}
}t[]; int main()
{
int n,q;
scanf("%d%d",&n,&q);
//for(int i=0;i<4;++i) t[i].build(1,n,1);
while(q--)
{
char s[];
scanf("%s",s);
int a,b,x;
if(s[]=='U')
{
scanf("%d%d%d",&a,&b,&x);
for(int i=;i<;++i)
{
int tt=;
while(x%prime[i]==)
{
x/=prime[i];
tt++;
}
//debug(tt);
if(tt!=) t[i].update(,n,,a,b,tt);
if(x==) break;
}
}
else
{
scanf("%d%d",&a,&b);
int ans[];
for(int i=;i<;++i)
{
ans[i]=t[i].query(,n,,a,b);
//debug(ans[i]);
}
int a1=max(ans[],ans[]),a2=max(ans[],ans[]);
printf("ANSWER %d\n",max(a1,a2));
}
}
return ;
}
  • I题

2~62任意进制的转换,poj1220

代码来源于大神,看了有些时间才差不多理解。

原来的代码是利用字符ASCII码直接映射的,我自己写的建立了两个互相的映射map,感觉也挺好用的。

转换的核心过程见注释代码。

 #include <bits/stdc++.h>
#define debug(x) cout<< # x <<": "<<x<<endl
using namespace std;
const int MAXN=1e3+; map<char,int> mp;
map<int,char> mp2; void init()
{
int k=;
char c;
for(int i=;i<;++i) //0~9
{
c=i;
mp2[k]=c;
mp[c]=k++;
}
for(int i=;i<+;++i) //A~Z
{
c=i;
mp2[k]=c;
mp[c]=k++;
}
for(int i=;i<+;++i) //a~z
{
c=i;
mp2[k]=c;
mp[c]=k++;
}
}
int n,m;
char s[MAXN],ans[MAXN];
int t[MAXN],a[MAXN]; void solve()
{
int len=strlen(s),k=;
for(int i=len-;i>=;--i)
{
t[len-i-]=mp[s[i]]; //倒置 得到每位的值 下标0表示低位
}
for(k;len;) //模拟 除m取余 的过程 一次得到低位的一个数
{
for(int i=len-;i>=;--i) //
{
t[i-]+=t[i]%m*n;
t[i]/=m;
}
a[k++]=t[]%m;
t[]/=m;
while(len && !t[len-]) //当前高位已为0,可以减少总位数
len--;
}
ans[k]='\0';
for(int i=;i<k;++i)
{
ans[k-i-]=mp2[a[i]]; //数是从下标0开始的,需要倒置
}
}
int main()
{
ios::sync_with_stdio(false);
cin.tie();
init();
cin>>n>>m;
cin>>s;
solve();
cout<<ans<<endl;
return ;
}

2019icpc银川站 复现赛的更多相关文章

  1. 2018 CCPC 桂林站(upc复现赛)补题

    2018 CCPC 桂林站(upc复现赛)补题 G.Greatest Common Divisor(思维) 求相邻数的差值的gcd,对gcd分解素因子,对所有的素因子做一次遍历,找出最小答案. 几个样 ...

  2. 2016ACM/ICPC亚洲区沈阳站-重现赛赛题

    今天做的沈阳站重现赛,自己还是太水,只做出两道签到题,另外两道看懂题意了,但是也没能做出来. 1. Thickest Burger Time Limit: 2000/1000 MS (Java/Oth ...

  3. 2013ACM/ICPC亚洲区南京站现场赛---Poor Warehouse Keeper(贪心)

    题目链接 http://acm.hdu.edu.cn/showproblem.php?pid=4803 Problem Description Jenny is a warehouse keeper. ...

  4. 2016 ACM/ICPC亚洲区青岛站现场赛(部分题解)

    摘要 本文主要列举并求解了2016 ACM/ICPC亚洲区青岛站现场赛的部分真题,着重介绍了各个题目的解题思路,结合详细的AC代码,意在熟悉青岛赛区的出题策略,以备战2018青岛站现场赛. HDU 5 ...

  5. 2016ACM/ICPC亚洲区大连站现场赛题解报告(转)

    http://blog.csdn.net/queuelovestack/article/details/53055418 下午重现了一下大连赛区的比赛,感觉有点神奇,重现时居然改了现场赛的数据范围,原 ...

  6. HDU 6227.Rabbits-规律 (2017ACM/ICPC亚洲区沈阳站-重现赛(感谢东北大学))

    Rabbits Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Others)Total S ...

  7. HDU 6225.Little Boxes-大数加法 (2017ACM/ICPC亚洲区沈阳站-重现赛(感谢东北大学))

    整理代码... Little Boxes Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/O ...

  8. 2019浙师大校赛(浙大命题)(upc复现赛)总结

    2019浙师大校赛(浙大命题)(upc复现赛)总结 早上九点开始.起得迟了,吃了早饭慌慌张张跑过去,刚到比赛就开始了. 开始分别从前往后和从后往前看题,一开始A题,第一发WA,第二次读题发现漏看了还有 ...

  9. 2019-ACM-ICPC-徐州站网络赛- I. query-二维偏序+树状数组

    2019-ACM-ICPC-徐州站网络赛- I. query-二维偏序+树状数组 [Problem Description] ​ 给你一个\([1,n]\)的排列,查询\([l,r]\)区间内有多少对 ...

随机推荐

  1. 【Android - 问题解决】之Android Studio在Refreshing Gradle Project卡住的问题

    今天把Android Studio从2.2版本升级到了2.3版本,但是再打开以前项目的时候,就卡在 “Refreshing XXX Gradle Project” 这里不动了. 网上搜了一些解决方案, ...

  2. MyEclispe启动Tomcat7时出现错误The servlets named [LoginServlet] and [com.liu.control.LoginServlet] are both

    刚开始尝试写Servlet代码,第一天就碰到这个错误,在网上找了很多资料才找到解决办法,在此记录一下. org.apache.catalina.LifecycleException: Failed t ...

  3. 机器学习-Python 01

    机器学习中最常用最流行的语言工具现阶段应该是Python, 这篇文章主要介绍一些常用的Python语法知识.本篇博文适合那些有其他语言基础的程序员们,如果一点基础都没有,我建议先跳过.博主以前是做移动 ...

  4. luogu P2515 [HAOI2010]软件安装 |Tarjan+树上背包

    题目描述 现在我们的手头有N个软件,对于一个软件i,它要占用Wi的磁盘空间,它的价值为Vi.我们希望从中选择一些软件安装到一台磁盘容量为MM计算机上,使得这些软件的价值尽可能大(即Vi的和最大). 但 ...

  5. OSS 对象存储的那些事AmazonS3简单使用

    对象存储是根据AmazonS3来做的封装,主要功能 :文件的上传下载 生成链接 对图片的处理 查看桶内对象等一系列的操作. 本文主要做的是文件的上传下载生成链接以及前期的准备工作 以springboo ...

  6. 2018 ICPC南京网络赛 Set(字典树 + 合并 + lazy更新)

    题解:n个集合,你要进行m个操作.总共有3种操作.第一种,合并两个集合x和y.第二张,把特定的集合里面所有的数字加一.第三种,询问在某个集合里面,对于所有数字对2的k次方取模后,有多少个数字等于x. ...

  7. Ceph 集群搭建

    1.Ceph部署 1.1 Ceph 安装前准备,环境与拓扑图如下: 主机名 外网IP 内网IP 安装组件 磁盘块 c720181 172.22.69.84 192.168.20.181 ceph-de ...

  8. ARTS-S ansible-playbook

    文件a.yml --- - hosts: cluster remote_user: ksotest gather_facts: false tasks: - name: delete dir if e ...

  9. Dubbo学习系列之七(分布式订单ID方案)

    既然选择,就注定风雨兼程! 开始吧! 准备:Idea201902/JDK11/ZK3.5.5/Gradle5.4.1/RabbitMQ3.7.13/Mysql8.0.11/Lombok0.26/Erl ...

  10. 【Mock MVC】使用学习

    [Mock MVC]使用学习 转载:https://www.cnblogs.com/yangchongxing/p/10658311.html