题目背景

没有背景

我写不出来了qwq

题目描述

Chino给定了nn个数a_1...a_na1​...an​,给定常数s,m,她会轮流对这nn个数做k组操作,每组操作包含以下几步:

1.swap(as​,am​)(交换a_s,a_m​)

2.2将n个数都向前平移一位(第11个移动到第n个位置上)

Chino想知道,k组操作后,这nn个数分别是多少?

Orz yky,dyh,wjk,jjy,cxr,gsy,cpy,zcy,tyz,yy,hz,zhr,yg

输入格式

第一行,四个数,n,s,m,k

接下来一行n个数,分别代表a_1,a_2...a_na1​,a2​...an​

输出格式

输出一行,n个数,分别代表a_1,a_2...a_na1​,a2​...an​

输入输出样例

输入 #1复制

4 1 2 3
1 2 3 4
输出 #1复制

1 2 3 4

说明/提示

所有数字均在long long以内

思路:

  矩阵加速递推,构造两个矩阵,一个存交换操作,一个存位移操作。

  假如交换1和2:

      0 1 0 0

      1 0 0 0

      0 0 0 1

      0 0 0 1

  位移操作,所有往前移一位,

         0 1 0 0

      0 0 1 0

      0 0 0 1

      1 0 0 0    

矩阵快速幂求解K次方,即操作次数:

jz ksm(jz a,int b)
{
jz anss;
memset(anss.c,0,sizeof(anss.c));
for(int i=1;i<=n;i++) anss.c[i][i]=1;
for(;b;b>>=1,a=a*a)
{
if(b&1)anss=a*anss;
// a=a*a;当时出错了,因为前面已经计算
}
return anss;
}

重载乘号:

struct jz{
int c[100][100];
}f,base,l1,l2;
int n,m,s,k;
jz operator * (const jz &a,const jz &b)
{
jz lin;
for(int i=1;i<=80;i++)
for(int j=1;j<=80;j++)
{
lin.c[i][j]=0;
for(int k=1;k<=80;k++)
{
lin.c[i][j]+=(a.c[k][j] * b.c[i][k]);
}
}
return lin;
}

 代码:

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cmath>
#include<cstring>
#define int long long
using namespace std;
struct jz{
int c[100][100];
}f,base,l1,l2;
int n,m,s,k;
jz operator * (const jz &a,const jz &b)
{
jz lin;
for(int i=1;i<=80;i++)
for(int j=1;j<=80;j++)
{
lin.c[i][j]=0;
for(int k=1;k<=80;k++)
{
lin.c[i][j]+=(a.c[k][j] * b.c[i][k]);
}
}
return lin;
}
void dy(jz x)//调试用的,可以忽略
{
for(int i=1;i<=n;++i)
{
for(int j=1;j<=n;++j)
cout<<x.c[i][j]<<" ";
cout<<endl;
}
}
jz ksm(jz a,int b)
{
jz anss;
memset(anss.c,0,sizeof(anss.c));
for(int i=1;i<=n;i++) anss.c[i][i]=1;
for(;b;b>>=1,a=a*a)
{
if(b&1)anss=a*anss;
// a=a*a;
}
return anss;
}
signed main()
{
scanf("%lld%lld%lld%lld",&n,&s,&m,&k);
for(int i=1;i<=n;i++) scanf("%lld",&f.c[i][1]);
for(int i=1;i<=n;i++) if(i!=s&&i!=m)l1.c[i][i]=1;
l1.c[s][m]=l1.c[m][s]=1;
for(int i=1;i<=n-1;i++)l2.c[i][i+1]=1;
l2.c[n][1]=1;
base=l1*l2;
base=ksm(base,k);
f=f*base;
for(int i=1;i<=n;i++)printf("%lld ",f.c[i][1]);
cout<<endl;
dy(l1);
cout<<endl;
dy(l2);
return 0;
}

  

感谢wlj_dy 的帮助

【luoguP5550】Chino的数列的更多相关文章

  1. Chino的数列

    题解: 一道练代码能力的题目.. 首先很显然他是一道平衡树裸题 第5个操作是势能分析维护最大值最小值就可以了 另外设置虚点和noip2017队列那题一样(不过我只写过线段树) 具体细节: 1.内存池, ...

  2. C#求斐波那契数列第30项的值(递归和非递归)

    using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.T ...

  3. BZOJ1500[NOI2005]维修数列

    Description Input 输入的第1 行包含两个数N 和M(M ≤20 000),N 表示初始时数列中数的个数,M表示要进行的操作数目.第2行包含N个数字,描述初始时的数列.以下M行,每行一 ...

  4. PAT 1049. 数列的片段和(20)

    给定一个正数数列,我们可以从中截取任意的连续的几个数,称为片段.例如,给定数列{0.1, 0.2, 0.3, 0.4},我们有(0.1) (0.1, 0.2) (0.1, 0.2, 0.3) (0.1 ...

  5. 斐波拉契数列加强版——时间复杂度O(1),空间复杂度O(1)

    对于斐波拉契经典问题,我们都非常熟悉,通过递推公式F(n) = F(n - ) + F(n - ),我们可以在线性时间内求出第n项F(n),现在考虑斐波拉契的加强版,我们要求的项数n的范围为int范围 ...

  6. fibonacci数列(五种)

    自己没动脑子,大部分内容转自:http://www.jb51.net/article/37286.htm 斐波拉契数列,看起来好像谁都会写,不过它写的方式却有好多种,不管用不用的上,先留下来再说. 1 ...

  7. js中的斐波那契数列法

    //斐波那契数列:1,2,3,5,8,13…… //从第3个起的第n个等于前两个之和 //解法1: var n1 = 1,n2 = 2; for(var i=3;i<101;i++){ var ...

  8. 洛谷 P1182 数列分段Section II Label:贪心

    题目描述 对于给定的一个长度为N的正整数数列A[i],现要将其分成M(M≤N)段,并要求每段连续,且每段和的最大值最小. 关于最大值最小: 例如一数列4 2 4 5 1要分成3段 将其如下分段: [4 ...

  9. 剑指Offer面试题:8.斐波那契数列

    一.题目:斐波那契数列 题目:写一个函数,输入n,求斐波那契(Fibonacci)数列的第n项.斐波那契数列的定义如下: 二.效率很低的解法 很多C/C++/C#/Java语言教科书在讲述递归函数的时 ...

随机推荐

  1. 20190728-Python爬取视频&切割视频&视频加水印

    1.视频爬取 1.下载视频的源码如下: import os import requests from bs4 import BeautifulSoup import threading from bj ...

  2. Python【函数使用技巧】

    写成“子函数+主函数”的代码结构,也是因为每个不同的功能封装在单独的函数代码中,方便后续修改.增删 import math # 变量key代表循环运行程序的开关 key = 1 # 采集信息的函数 d ...

  3. 第6章:使用Python监控Linux系统

    1.Python编写的监控工具 1).多功能系统资源统计工具dstat dstat是一个用Python编写的多功能系统资源统计工具,用来取代Linux下的vmstat,iostat,netstat和i ...

  4. 用shell脚本安装MySQL-5.7.22-官方版本

    Install_CentOS7_MySQL57_binary.sh #!/bin/bash MySQL_Package=mysql-5.7.22-linux-glibc2.12-x86_64.tar. ...

  5. STM32之DMA实例

    版权声明:本文为博主原创文章,未经博主允许不得转载. https://blog.csdn.net/zouleideboke/article/details/75092558 DMA简介: DMA(Di ...

  6. SWATS算法剖析(自动切换adam与sgd)

    SWATS算法剖析(自动切换adam与sgd) 战歌指挥官 搬砖.码砖.代查水表.... 27 人赞同了该文章 SWATS是ICLR在2018的高分论文,提出的一种自动由Adam切换为SGD而实现更好 ...

  7. Spring MVC 探讨DispatcherServlet

    先上DispatcherServlet的运行流程图(request processing):

  8. ASP.NET Core 2.0 中读取 Request.Body 的正确姿势

    原文:ASP.NET Core 中读取 Request.Body 的正确姿势 ASP.NET Core 中的 Request.Body 虽然是一个 Stream ,但它是一个与众不同的 Stream ...

  9. js Indexof的用法

    JavaScript中indexOf()函数 JavaScript中indexOf()函数方法返回一个整数值,指出  String 对象内子字符串的开始位置.如果没有找到子字符串,则返回 -1.如果 ...

  10. C#/.net 通过js调用系统相机进行拍照,图片无损压缩后进行二维码识别

    这两天撸了一个需求,通过 JS  调用手机后置相机,进行拍照扫码.前台实现调用手机相机,然后截取图片并上传到后台的功能.后台接收传过来的图片后,通过调用开源二维码识别库 ZXing 进行二维码数据解析 ...