6月16日省中集训题解

最新推荐文章于 2019-09-20 20:35:02 发布
最新推荐文章于 2019-09-20 20:35:02 发布
阅读量252 收藏
点赞数 1
分类专栏： 学习总结&其他 数学 树状数组 暴力
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/even_bao/article/details/80719413
版权

Problem A Pailindrome

【题意】

            求在k进制下n位回文数的和，答案对1000001取模，结果用十进制输出

            n <= 10^18

            k <= 10^6

            对于20%的数据，n <= 16,k <= 20

            对于50%的数据，n <= 10^4,k <= 10^3

            对于100%的数据,n <= 10^18,k <= 10^6

【题解】

          满分算法： 

                   解法1 ：

                   我们知道对于一个k进制下的n位回文数S(假设n为偶数），一定有 ： S1 = Sn, S2 = Sn-1, .. Sn/2 = Sn/2+1

                   那么，我们只要知道前n/2位数，就能够构造出一个k进制下的n位回文数

                   假设我们知道前n/2位数，那么其对应的回文数在10进制下的值为 :

                   S1k^(n-1) + S2k^(n-2) + ... +  S(n/2)k^(n/2) + S(n/2)k^(n/2-1) + S(n/2-1)k^(n/2-2) + ... + S1

                   = S1(k^(n-1) + 1) + S2(k^(n-2) + k) + ... + S(n/2)(k^(n/2-1) + k^(n/2-2))

                   考虑第一项对答案的"贡献" :

                   当S1在1..k-1范围内时才会对答案产生"贡献"，而S2..S(n/2)，每个数都有0..k-1这k种取值，因此，它对答案

                   的贡献为 : (1 + 2 + 3 + ... + (k - 1) ) * k^(n/2-1) * (k^(n-1) + 1)

                                   = k(k-1)/2k^(n/2-1)(k^(n-1)+1) = (k - 1)k^(n/2)/2(k^(n-1)+1)

                   考虑第2,3,4...n/2项对答案的贡献 ：

                   S1有(k-1)中取值，其它的数有k种取值，

                   所以，它对答案的贡献为 ：

                   （k-1)^2k^(n/2-1)/2(k ^ (n - 1) + k^0 + k^(n-2) + k + ... + k^(n/2-1) + k^(n/2-2))

                     = (k-1)^2k^(n/2-1)/2(k^0 + k + k^2 + ... + k^(n-2) - k^(n/2))

                   n为奇数时同理

                   那么，贡献相加就是最后的答案，注意要用快速幂辅助运算

                   k^0 + k + .. + k^(n-2)是一个等比数列，可以用二分法求 :

                   设f(n) = S + S^2 + .. S^n

                   那么，当n = 1时，f(n) = S

                             当n为奇数时, f(n) = f(n-1) + S^n

                             当n为偶数时,f(n) = f(n/2) + S^(n/2)f(n/2)

                   解法2 ：

                   在k进制下的n位回文数共有 : (k - 1)k^((n+1)/2)个

                   根据这个结论，对第一位，第二位，...第(n+1)/2位讨论，即可，笔者不再赘述

【代码】

          

#include<bits/stdc++.h>
using namespace std;
const long long P = 1000001;

long long s,n,k,ans1,ans2,ans3,ans;

template <typename T> inline void read(T &x)
{
    int f = 1; x = 0;
    char c = getchar();
    for (; !isdigit(c); c = getchar()) { if (c == '-') f = -f; }
    for (; isdigit(c); c = getchar()) x = (x << 3) + (x << 1) + c - '0';
    x *= f;
}
template <typename T> inline void write(T x)
{
    if (x < 0)
    {
        putchar('-');
        x = -x;
    }
    if (x > 9) write(x/10);
    putchar(x%10+'0');
}
template <typename T> inline void writeln(T x)
{
    write(x);
    puts("");
}

inline long long power(long long a,long long n)
{
		long long ans = 1,b = a;
		while (n)
		{
				if (n & 1) ans = (ans * b) % P;
				b = (b * b) % P;
				n >>= 1;	
		}	
		return ans;
}
inline long long calc(long long n)
{
		long long tmp;
		if (n == 1) return k % P;
		else
		{
				if (n % 2 == 1) return (calc(n-1) + power(k,n)) % P;
				else 
				{
						tmp = calc(n/2);
						return (tmp + (power(k,n/2) * tmp) % P) % P;
				}
		}
}
 
int main() 
{
		
		freopen("palindrome.in","r",stdin);
		freopen("palindrome.out","w",stdout);
		
		read(k); read(n);
		if (n & 1)
		{
				s = k * (k - 1) / 2 % P;
				ans1 = ((s * power(k,n/2)) % P * (power(k,n-1) + 1) % P) % P;
				ans2 = ((s * power(k,n/2-1) * (k - 1) % P) % P * ((calc(n-2) - power(k,n/2) + P) % P)) % P ;	
				ans3 = ((s * (k - 1) % P) * power(k,n-2)) % P;
				ans = (((ans1 + ans2) % P) + ans3) % P; 
		} else
		{
				s = k * (k - 1) / 2 % P;
				ans1 = s * power(k,n/2-1) % P * (power(k,n-1) + 1) % P;
				ans2 = (s * (k - 1) % P * power(k,n/2-2)) % P * calc(n-2) % P;
				ans = (ans1 + ans2) % P;	
		}
		writeln(ans);
		
		return 0;
	
}

Problem B wave

            【题意】

                     给定序列a

                     求有多少个子序列满足 ：

                     Ai != Aj

                     i < j < k且Ai < Aj,Ak < Aj或Ai > Aj,Ak > Aj

            【题解】

                     满分算法 ：

                     先将序列离散化

                     对于第一种序列，我们用Cnt[i]表示i出现了多少次，那么，第i个位置对答案的“贡献”为 : i - Cnt[A[i]]

                     对于第二种序列，我们枚举j，根据乘法原理，j这个位置对答案的“贡献”为 ：

                     左边小于Aj数的个数 * 右边小于Aj数的个数 + 左边大于Aj数的个数 * 右边大于Aj数的个数

                     可以用树状数组优化

              【代码】

                        此题常数卡得很紧，写程序时注意 ：

                        1. 能写在一个循环里的不要分到多个循环里写

                        2. long long自带4倍常数，能用int的不要用long long 

                        3. 比较大小cmp函数，用三目运算比if语句快

                        

#include<bits/stdc++.h>
using namespace std;
const int P = 1e9 + 7;
const int MAXN = 1e6 + 5;

int i,n,len;
int tmp[MAXN],a[MAXN],rk[MAXN],s[MAXN];
int la[MAXN],lb[MAXN],ra[MAXN],rb[MAXN];
long long ans;

template <typename T> inline void read(T &x)
{
		int f = 1; x = 0;
		char c = getchar();
		for (; !isdigit(c); c = getchar()) { if (c == '-') f = -f; }
		for (; isdigit(c); c = getchar()) x = (x << 3) + (x << 1) + c - '0';
		x *= f; 
}
template <typename T> inline void write(T x)
{
		if (x < 0)
		{
			putchar('-');
			x = -x;	
		}	
		if (x > 9) write(x/10);
		putchar(x%10+'0');
}
template <typename T> inline void writeln(T x)
{
		write(x);
		puts("");
}
class BinaryIndexedTree
{
		private :
				int c[MAXN];
		public :
			inline int lowbit(int x)
			{
					return x & (-x);
			}
			inline void clear()
			{
					memset(c,0,sizeof(c));
			}
			inline void modify(int pos,int val)
			{
					int i;
					for (i = pos; i <= n; i += lowbit(i)) c[i] = (c[i] + val) % P;
			}	
			inline int query(int x)
			{
					int i;
				  int ans = 0;
					for (i = x; i; i -= lowbit(i)) ans = (ans + c[i]) % P;
					return ans;
			}
			inline int query(int l,int r)
			{
					return (query(r) - query(l-1) + P) % P;
			}
} BIT;

int main()
{
	
		freopen("wave.in","r",stdin);
		freopen("wave.out","w",stdout);
		
		read(n);
		for (i = 1; i <= n; i++) 
		{
				read(a[i]);
				tmp[i] = a[i];
		}
		sort(tmp+1,tmp+n+1);
		len = unique(tmp+1,tmp+n+1) - tmp - 1;
		for (i = 1; i <= n; i++) rk[i] = lower_bound(tmp+1,tmp+len+1,a[i]) - tmp;
		BIT.clear();
		for (i = 1; i <= n; i++)
		{
				s[rk[i]]++;
				ans = (ans + i - s[rk[i]] + P) % P;
				la[i] = BIT.query(rk[i]-1);
				lb[i] = BIT.query(rk[i]+1,len);
				BIT.modify(rk[i],1);	
		}
		BIT.clear();
		for (i = n; i >= 1; i--)
		{
				ra[i] = BIT.query(rk[i]-1);
				BIT.modify(rk[i],1);
				rb[i] = BIT.query(rk[i]+1,len);
		}
		for (i = 1; i <= n; i++) ans = (ans + 1ll * la[i] * ra[i] + 1ll * lb[i] * rb[i]) % P;
		writeln(ans);
	
	return 0;
}

Problem C String

【题意】

           给出一个长度为N的字符串S，只包含A、C、G、U四种大写字母。A可以和U配对，C可以和G配对。定义一个长度        x（1<=x<=N）可以作为S的界，当且仅当S[1..x]能和S[N-x+1..N]每一位都成功配对。求S所有界中的最小值和最大值。

【题解】

          暴力判断x是否S的“界”，即可

【代码】

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 1e5 + 5;

int i,n,ansa,ansb;
char s[MAXN];

inline bool match(int x)
{
		int i;
		for (i = 1; i <= x; i++)
		{
				if (s[i] == 'A' && s[n-x+i] == 'U') continue;
				if (s[i] == 'U' && s[n-x+i] == 'A') continue;
				if (s[i] == 'C' && s[n-x+i] == 'G') continue;
				if (s[i] == 'G' && s[n-x+i] == 'C') continue;
				return false;
		}
		return true;
}
inline void solvea()
{
		int i;
		for (i = 1; i <= n; i++)
		{
			if (match(i))
			{
					ansa = i;
					return;
			}
		}
}
inline void solveb()
{
		int i;
		for (i = n; i >= 1; i--)
		{
				if (match(i))
				{
						ansb = i;
						return;
				}
		}
}

int main(){
	
	freopen("string.in","r",stdin);
	freopen("string.out","w",stdout);
	
	scanf("%d%s",&n,s+1);
	solvea();
	solveb();
	printf("%d %d\n",ansa,ansb);
	
	return 0;
}


          

           

even_bao
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
Jzoj3654 【APIO2014】回文串
扩展の灰(Extended Ash/Cooevjnz)
10-21 440
题意：求出现次数*长度最大的回文串 很水的一道题直接用回文树可以切了 太菜并不会回文树，用了极慢的sam+manacher manacher求出所有回文子串用sam在parent树上面倍增就可以求出一个串的出现次数，manacher是O(n)的，倍增lgn，整体复杂度nlgn （分开的manacher比什么插入特殊符号的好写多了） #include #include #include
P1435 [IOI2000]回文字串
zjsru_Beginner的博客
08-06 396
题目描述 回文词是一种对称的字符串。任意给定一个字符串，通过插入若干字符，都可以变成回文词。此题的任务是，求出将给定字符串变成回文词所需要插入的最少字符数。 比如 “Ab3bd”插入2个字符后可以变成回文词“dAb3bAd”或“Adb3bdA”，但是插入少于2个的字符无法变成回文词。 注：此问题区分大小写 输入格式 一个字符串(0<strlen<=1000) 输出格式 有且只有一个整数，即最少插入字符数 输入输出样例 输入 #1 Ab3bd 输出 #1 2 解题思路
「LOJ6073」「2017 山东一轮集训 Day5」距离-主席树+树链剖分
DSL_HN_2002的博客
12-31 386
Description 给你一棵 nnn 个点的树和一个排列 ppp,边有边权,记 dist(u,v)dist(u, v)dist(u,v) 表示 uuu 到 vvv 的距离,path(u,v)path(u, v)path(u,v) 表示 uuu 到 vvv 路径上所有点组成的集合,现在有 qqq 次询问,每次给出 uiu_iui​ , viv_ivi​ , kik_iki​ ,问: ∑i∈pat...
【题解】sdoj3738[2018.8.8集训]A.海星突击队 构造
翻过这座山，他们就会听到你的故事
08-09 179
问题描述 海星突击队在一次执行任务中，有 n 个人在敌占区被围困，他们最后在一个隐蔽地 方安营扎寨，等待友军支援。 从安全角度出发，每天晚上他们必须安排人来值班，以防敌人发现袭击。他们的长 官决定每天安排至少 3 人值班才能保证基本安全。但是因为晚上值班过于无聊，值班的人 可能发生冲突，所以不能让任何 2 个人一起值班超过 3 次（不一定连续 3 次）。 请问突击队最多能度过多少天，并...
菲波那契数列
kt1776133839的博客
02-26 369
#include&lt;bits/stdc++.h&gt;using namespace std;int n,a[1000001]={1,1},m;int main(){for(int i=2;i&lt;=1000000;i++){ a[i]=(a[i-1]+a[i-2])%1000;}cin&gt;&gt;n;while(n--){cin&gt;&gt;m;cout&lt;...
css_02 | CSS——CSS 选择器详解
weixin_34246551的博客
04-13 409
本文推荐 PC 端阅读~ 本文版权归 “公众号 | 前端一万小时” 所有，未经授权，请勿转载！ 复制代码 css_02 复制代码 1. 伪类选择器有哪些？ 2. 伪元素和伪类的区别？ 3. 关于 css 选择器，以下说法正确的是？（不定项） ❌ header .p 是选择 id 为 header 的类型为 p 的所有后代。 ✅ header > .p 是选择标签为 header 内...
LOJ#6075. 「2017 山东一轮集训 Day6」重建
weixin_30346033的博客
07-13 146
题目描述： 　　给定一个 n个点m条边的带权无向连通图，以及一个大小为k的关键点集合S。有个人要从点s走到点t,现在可以对所有边加上一个非负整数a,问最大的a，使得加上a后，满足：s到t的最短路长度=s到t且只能经过S中的点的最短路长度。 题目分析： 　　暴力 　　记x为只经过关键点的最短路长度，其路径条数为n 　　记y为可经过任意点的最短路长度，其路径条数为m 　　　...
linux命令大全（完善中）
Rio520的博客
09-20 2170
bash作为大多数的linux系统的默认字符解析器，所以必须了解bash的优势 1.默认保存历史命令（可以用上下键翻看） 2.命令输入前几位可用tab键补全 3.强大的批量处理脚本 4.实用的环境变量 一条完整命令的格式 命令 字符 【命令参数】 【操作对象】 命令的参数可选用长格式（完整的选项名称）也可选用短格式（单个字母的缩写），分别用“--”；“-”做前缀...
【NOIP2017Day1T1】 小凯的疑惑
even_bao的博客
03-28 677
【题目链接】          点击打开链接【算法】          px + qy不能表示的最大整数为 pq - p - q          证明见这篇博客，过程很详细，推荐阅读 ：          https://blog.csdn.net/qwerty1125/article/details/78661916【代码】          var a,b : int64; begi...
【扬中集训DAY1T1】 微信群
even_bao的博客
02-22 431
【题目链接】          点击打开链接【算法】         对问题稍加分析后，发现其实要求的就是 ：         C(N,K) + C(N,K+1) + C(N,K+2) + ... + C(N,N)         因为N最大10^9,K最大10^5，很显然，我们不可能通过暴力或杨辉三角等来解决此题         我们发现 ：         C(N,K) + C(N,K+1) ...
【NOIP16提高组】换教室
even_bao的博客
02-15 420
 【题目链接】          点击打开链接 【算法】          概率DP          先跑一遍floyed，求出每个教室之间的最短路径，存在数组dist[][]中，时间复杂度O(V^3)          设计状态，f[i][j][k]表示当前选到第i个教室，已经选了j个教室，当前这个教室选不选（0..1)          那么，状态转移方程是什么呢？          假设当...
【Codeforces 582A】 GCD Table
even_bao的博客
03-13 406
【题目链接】          点击打开链接【算法】          G中最大的数一定也是a中最大的数。          G中次大的数一定也是a中次大的数。          第三、第四可能是由最大和次大的gcd产生的         那么就不难想到下面的算法：         1. 令p为G中最大的数。在G中删除p，a中加入p。         2 . 对于a中的所有其他数(设为q)，在G中...
【SDOI2009】SuperGCD
even_bao的博客
03-18 362
【题目链接】           点击打开链接【算法】          1.关于求最大公约数的算法          若使用辗转相除法，那么显然会超时          不妨这样思考 ：          要求gcd(a,b),          若a为偶数，b为偶数，则gcd(a,b) = 2 * gcd(a/2,b/2)          若a为偶数，b为奇数，则gcd(a,b) = gcd...
【HDU 5698】 瞬间移动
even_bao的博客
05-27 348
【题目链接】           点击打开链接【算法】          用f[i][j]表示走到(i,j)这个位置有多少种方案，因为走到(i,j)这个位置，上一步一定在它左上角的矩形中，所以，          f(i,j) = sigma( f(x,y) ) ( (x,y)在左上角的矩形中）           我们尝试将它画出来，发现是斜着的杨辉三角          然后，通过找规律，我们...
【扬中集训Day6T1】 白日梦
even_bao的博客
02-14 333
【题目描述】          白日梦（daydream.c/cpp/pas）时间限制: 1 s  空间限制: 256 MB题目描述    SR需要相当大的睡眠量    某日，他做了一个奇怪的梦，他梦见自己成为了怪物猎人，为了狩猎，他要去一个岛上住N+1天（编号为0到N）。这个岛位于太平洋中心，每天要么是晴天，要么刮台风。他到达岛的第0天是晴天。然后对于第i天，假如是晴天，那么有P(0&lt;p&...
【POJ 1061】 青蛙的约会
even_bao的博客
02-14 291
【题目链接】           点击打开链接【算法】         列出同余方程，然后用exgcd求解【代码】          #include &lt;algorithm&gt; #include &lt;bitset&gt; #include &lt;cctype&gt; #include &lt;cerrno&gt; #include &lt;clocale&gt; #include ...
【HAOI 2008】 硬币购物
even_bao的博客
05-06 286
【题目链接】          点击打开链接【算法】          此题是一道好题！          首先，我们发现 : 付款方法数 = 不受限制的方法数 - 受限制的方法数          那么，我们怎么求呢？          我们用dp求出不受限制的方法数（f[i]表示买i元的东西，不受硬币限制，有多少种方案），只需用01背包的          方法就可以了，实现非常简单     ...
【2017省中集训DAY1T1】 小X的质数
even_bao的博客
02-15 276
【题目链接】          点击打开链接【算法】         如果一个数是小X喜欢的数，那么有两种可能：         1.这个数是质数         2.这个数除以它的最小质因子是一个质数         所以我们可以用线性筛+前缀和的方式预处理，询问的时候O（1）计算就可以了【代码】        #include&lt;bits/stdc++.h&gt; using namesp...
"湖南省计算机等级考试VB程序设计题解
这份文档包含了湖南省计算机等级考试中涉及的VB程序设计题目，每道题目均会进行严格的分析和解答。 二、题目解答 1. 题目一 题目一要求使用VB程序设计实现一个简单的计算器功能。首先，需要设计一个界面，包括数字...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值