csps-s模拟测试60嘟嘟噜,天才绅士少女助手克里斯蒂娜,凤凰院凶真题解
题面:https://www.cnblogs.com/Juve/articles/11625190.html
嘟嘟噜:
约瑟夫问题
第一种递归的容易re,但复杂度较有保证
第二种适用与n大于m的情况
第三种O(n)用于n不太大或m大于n时
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define int long long
#define re register
using namespace std;
int t,n,m;
int calc(int n,int m){
if(n==1) return 0;
if(n<m) return (calc(n-1,m)+m)%n;
int s=calc(n-n/m,m)-n%m;
return s<0?s+n:s+s/(m-1);
}
signed main(){
scanf("%lld",&t);
while(t--){
scanf("%lld%lld",&n,&m);
cout<<calc(n,m)+1<<endl;
re int ans=0;
int i=2;
while(i<=n){
int tim=(i-ans-1)/m+1;
if(i-1+tim>=n) tim=n-(i-1);
i+=tim,ans=(ans+tim*m)%(i-1);
}
cout<<ans+1<<endl;
ans=0;
for(re int i=2;i<=n;++i){
ans=(ans+m)%i;
}
printf("%lld\n",ans+1);
}
return 0;
}
天才绅士少女助手克里斯蒂娜:
就是推式子:
$\sum\limits_{i=l}^{r}\sum\limits_{j=i+1}^{r}|v_i*v_j|^2=\sum\limits_{i=l}^{r}\sum\limits_{j=i+1}^{r}(x_i^2*y_j^2+y_i^2*x_j^2-2*x_i*y_i*x_j*y_j)$
我们把所有$x_i^2$和$y_i^2$和$x_i*y_i$提出来,就有了
$ans=\sum\limits_{i=l}^{r}(x_i^2*\sum\limits_{j=i+1}^{r}y_j^2+y_i^2*\sum\limits_{j=i+1}^{r}x_j^2+x_i*y_i*\sum\limits_{j=i+1}^{r}x_j*y_j)$
然后这个式子可以树状数组维护后面的sigma
但是复杂度还是不优
我们可以把i和j看成无序的,最后再除以2
所以就有了:
$\sum\limits_{i=1}^{n}\sum\limits_{j=i+1}^{n}|v_i*v_j|^2=\sum\limits_{j=1}^{n}(x_j^2)*\sum\limits_{j=1}^{n}(y_j^2)-\sum\limits_{j=1}^{n}(x_j*y_j)^2$
然后就可以愉快地树状数组了
#include<iostream>
#include<cstdio>
#include<iostream>
#include<algorithm>
#define int long long
#define re register
using namespace std;
const int MAXN=1e6+5;
const int mod=20170927;
int n,m;
int read(){
int x=0;char ch=getchar();
while(ch<'0'||ch>'9') ch=getchar();
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return x;
}
struct node{
int x,y,s1,s2,s3;
inline friend int operator * (node p,node q){
return ((p.x*q.y%mod-q.x*p.y%mod)%mod+mod)%mod;
}
}a[MAXN];
inline int lowbit(re int x){
return x&(-x);
}
struct BIT{
int c[MAXN];
inline void update(re int pos,re int val){
for(re int i=pos;i<=n;i+=lowbit(i)){
(c[i]+=val)%=mod;
}
}
inline int query(re int pos){
re int res=0;
for(re int i=pos;i>0;i-=lowbit(i)){
(res+=c[i])%=mod;
}
return res;
}
inline int ask(re int l,re int r){
return ((query(r)-query(l-1))%mod+mod)%mod;
}
}tr[3];
signed main(){
n=read(),m=read();
for(re int i=1;i<=n;++i){
a[i].x=read(),a[i].y=read();
a[i].s1=a[i].x*a[i].x%mod,a[i].s2=a[i].y*a[i].y%mod,a[i].s3=a[i].x*a[i].y%mod;
tr[0].update(i,a[i].s1),tr[1].update(i,a[i].s2),tr[2].update(i,a[i].s3);
}
while(m--){
re int opt=read();
if(opt==1){
re int p=read(),x=read(),y=read();
tr[0].update(p,x*x%mod-a[p].s1);
tr[1].update(p,y*y%mod-a[p].s2);
tr[2].update(p,x*y%mod-a[p].s3);
a[p]=(node){x%mod,y%mod,x*x%mod,y*y%mod,x*y%mod};
}else{
re int l=read(),r=read(),ans=0;
ans=((tr[0].ask(l,r)*tr[1].ask(l,r)%mod-tr[2].ask(l,r)*tr[2].ask(l,r)%mod)%mod+mod)%mod;
printf("%lld\n",ans);
}
}
return 0;
}
凤凰院凶真:
求两个数列的最长公共上升子序列并输出路径
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stack>
#define int long long
using namespace std;
const int MAXN=5005;
int n,a[MAXN],m,b[MAXN],f[MAXN][MAXN],pre[MAXN][MAXN],pos=1;
stack<int>sta;
signed main(){
scanf("%lld",&n);
for(int i=1;i<=n;++i) scanf("%lld",&a[i]);
scanf("%lld",&m);
for(int i=1;i<=m;++i) scanf("%lld",&b[i]);
memset(pre,-1,sizeof(pre));
for(int i=1;i<=n;++i){
int v=0,k=0;
for(int j=1;j<=m;++j){
f[i][j]=f[i-1][j];
if(b[j]<a[i]&&v<f[i-1][j]) v=f[i-1][j],k=j;
if(a[i]==b[j]) f[i][j]=v+1,pre[i][j]=k;
}
}
for(int i=1;i<=n;++i){
if(f[n][pos]<f[n][i]) pos=i;
}
printf("%lld\n",f[n][pos]);
if(f[n][pos]!=0){
for(int i=n;i>=1;--i){
if(pre[i][pos]!=-1)
sta.push(a[i]),pos=pre[i][pos];
}
while(!sta.empty()){
printf("%lld ",sta.top());
sta.pop();
}
puts("");
}
return 0;
}
csps-s模拟测试60嘟嘟噜,天才绅士少女助手克里斯蒂娜,凤凰院凶真题解的更多相关文章
- [CSP-S模拟测试]:天才绅士少女助手克里斯蒂娜(数学+树状数组)
题目描述 红莉栖想要弄清楚楼下天王寺大叔的显像管电视对“电话微波炉(暂定)”的影响. 选取显像管的任意一个平面,一开始平面内有个$n$电子,初始速度分别为$v_i$,定义飘升系数为$$\sum \li ...
- 【CSP模拟赛】天才绅士少女助手克里斯蒂娜(线段树&读入优化&输出优化)
题面描述 红莉栖想要弄清楚楼下天王寺大叔的显像管电视对“电话微波炉(暂定)”的影响.选取显像管的任意一个平面,一开始平面内有个n电子,初始速度分别为vi,定义飘升系数为 $$\sum_{1\leqsl ...
- [CSP-S模拟测试]:凤凰院凶真(LCIS)
题目描述 $\alpha$世界线.凤凰院凶真创立了反抗$SERN$统治的组织“瓦尔基里”.为了脱离$\alpha$线,他需要制作一个世界线变动率测量仪.测量一个世界线相对于另一个世界线的变动率,实质上 ...
- 【CSP模拟赛】凤凰院凶真(最长公共上升子序列)
题目描述 α世界线.凤凰院凶真创立了反抗SERN统治的组织“瓦尔基里”.为了脱离α线,他需要制作一个世界线变动率测量仪. 测量一个世界线相对于另一个世界线的变动率,实质上就是要求出这两个世界线的最长公 ...
- csp-s模拟测试60
csp-s模拟测试60 2019-10-05 RT. 又颓又垃圾. 状态低迷,题都交不上去. 交了也是爆零,垃圾玩家没有什么可说的,就是垃圾. A. 嘟嘟噜 $mlogn$的毒瘤做法. 贴 ...
- [CSP-S模拟测试60]题解
回去要补一下命运石之门了…… A.嘟嘟噜 给定报数次数的约瑟夫,递推式为$ans=(ans+m)\% i$. 考虑优化,中间很多次$+m$后是不用取模的,这种情况就可以把加法变乘法了.问题在于如何找到 ...
- [考试反思]1005csp-s模拟测试60:招魂
最近总是好一场烂一场的.没有连续两场排名波动小于20的... 没人管.反正大脸一点脸没有就又AK了. 但是T3爆零这种事情吧... 爆搜不是很难打,但是想优化想了半天剩的时间不够结果赶忙打出来了,然后 ...
- csp-s模拟测试99
csp-s模拟测试99 九九归一直接爆炸. $T1$一眼板子. $T2$一眼语文题(语文的唯一一次$120+$是给模拟出来的可知我的语文能力). $T3$一眼普及题. ?? Hours Later 板 ...
- csp-s模拟测试92
csp-s模拟测试92 关于$T1$:最短路这一定建边最短路. 关于$T2$:傻逼$Dp$这一定线段树优化$Dp$. 关于$T3$:最小生成树+树P+换跟一定是这样. 深入(?)思考$T1$:我是傻逼 ...
随机推荐
- python下pip使用bug汇总
PS:以下操作全部基于win10 64位操作系统 pip安装任何包都出现问题: Cannot unpack file /tmp/pip-KzJgHD-unpack/simple 报错: Cannot ...
- ATC/TC/CF
10.25 去打 CF,然后被 CF 打了. CF EDU 75 A. Broken Keyboard 精神恍惚,WA 了一发. B. Binary Palindromes 比赛中的憨憨做法,考虑一个 ...
- Petrozavodsk Summer-2016. Warsaw U Contest, XVI Open Cup Onsite.
Petrozavodsk Summer-2016. Warsaw U Contest, XVI Open Cup Onsite. Problem A. Gambling Problem B. Colo ...
- 【转】Java程序CPU飙升问题排查方法
windows环境下cpu飙升问题 线上某台runtime机器(windows Server)cpu报警,这种情况初步就是代码里面死循环了,先把机器下线了保证不再有新的任务分配进来,然而cpu使用依然 ...
- css 写一个向右的箭头
经常写移动端页面时会用到向右的箭头,之前都是用图片,突然想到用css应该也能写,于是一搜还确实有 width: 7px; height: 7px; border-top: 2px solid #fff ...
- SQLserver执行命令
方法一:xp_cmdshell exec master..xp_cmdshell "whoami"默认执行是关闭 EXEC sp_configure 'show advanced ...
- JS函数 编程练习 使用javascript代码写出一个函数:实现传入两个整数后弹出较大的整数。
编程练习 使用javascript代码写出一个函数:实现传入两个整数后弹出较大的整数. 任务 第一步: 编写代码完成一个函数的定义吧. 第二步: 我们来补充函数体中的控制语句,完成函数功能吧. 提示: ...
- React和vue的差异和相似地方
React 单向绑定(加插件后,还是可以双向绑定) Vue 双向绑定 组件化 1. React,需要编写render函数, 2. 当React状态的状态state改变是render就会重新被调用, 重 ...
- redis安装配置使用
1.redis官方下载地址:https://redis.io/download 也可以github上下载,redis 64位下载地址:https://github.com/ServiceStack/r ...
- 树莓派3B+ 人脸识别、摄像头安装和使用
最近在学校里折腾树莓派上的人脸识别,折腾了很久才能用 在此记录下使用的过程和遇到的困难 过程基于超有趣!手把手教你使用树莓派实现实时人脸检测完成的.其中前面opencv的安装是文章中的Raspbian ...