9.13——模拟赛

T1 全排列

(permutation.cpp/c/pas)
Description
从 n 个不同元素中任取 m(m≤n)个元素,按照一定的顺序排列起来,叫做从 n
个不同元素中取出 m 个元素的一个排列。当 m=n 时所有的排列情况叫全排列。
你觉得 xxy 会问你全排列的个数吗?Xxy:这个问题能淹死你,我才不问呢。我
要问的是求 n 的全排列中,先递增后递
减、先递减后递增的全排列的个数。由于答案可能很大,对 p 取余
Input
输入包含多组测试数据每组测试
数据一行两个整数 n,p
Output
对于每组测试数据输出一行表示答案
Example
permutation.in permutation.out
3 5 4
2 233 0
Hint
设数据组数为 T
对于 10%的数据,n<=10,p<=1000,T=1
对于另外 10%的数据,n<=12,p<=1000,T<=12
对于另外 10%的数据,n<=100,p<=100000,T<=100
对于另外 10%的数据,n<=100000,p<=1000000,T=1
对于另外 10%的数据,n<=100000,p<=1000000,T<=1000
对于另外 20%的数据,n<=1e9,p<=1e9,T<=1000
对于 100%的数据,n<=1e18,p<=1e18,T<=1000

 1 #include <cstdio>
 2
 3 #define LL long long
 4 inline void read(LL &x)
 5 {
 6     x=0; register char ch=getchar();
 7     for(; ch>‘9‘||ch<‘0‘; ) ch=getchar();
 8     for(; ch>=‘0‘&&ch<=‘9‘; ch=getchar()) x=x*10+ch-‘0‘;
 9 }
10
11 LL n,p;
12 inline LL mul(LL a,LL b,LL mod)
13 {
14     LL ret=0;
15     for(; b; b>>=1)
16     {
17         if(b&1) (ret+=a) %=mod;
18         (a<<=1) %=mod;
19     }
20     return ret;
21 }
22 inline LL Quick_pow(LL b,LL mod)
23 {
24     LL ret=1,base=2;
25     for(; b; b>>=1)
26     {
27 //        if(b&1) ret=((ret%mod)*(base%mod))%mod;
28 //        base=((base%mod)*(base%mod))%mod;
29         if(b&1) ret=mul(ret,base,mod);
30         base=mul(base,base,mod);
31     }
32     return (ret-4+mod)%mod;
33 }
34 int Persist()
35 {
36 //    freopen("permutation.in","r",stdin);
37 //    freopen("permutation.out","w",stdout);
38     for(;~scanf("%I64d%I64d",&n,&p); )
39         if(n==1) puts("0");
40         else printf("%I64d\n",Quick_pow(n,p));
41     return 0;
42 }
43
44 int Aptal=Persist();
45 int main(){;}

打表找规律+快速幂+快速乘

T2  埃及分数(题目要求多余codevs)

(egypt.cpp/c/pas)
Description
对于一个分数 a/b(a!=1),将它表示为 1/x + 1/y + 1/z ……的形式,x,
y,z……互不相同。
多解取加数少的。加数相同时,取最小的分数最大的,最小分数相同时,取次小分
数最大的,以此类推。
输入保证 a<b 且 gcd(a,b)=1
Input
输入包含多组测试数据每组测试数
据包含一行 2 个数 a,b
Output
对于每组测试数据,输出一行表示答案,只输出分母,每个分母之间空 1 格
从小到大输出
Example
egypt.in egypt.out
5 6 2 3
8 9 2 3 18
Hint
对于 10%的数据,a,b<=10
对于另外 10%的数据,a,b<=100
对于另外 20%的数据,a,b<=10000
对于另外 30%的数据,a,b<=100000
对于 100%的数据,a,b<=1000000
由于本题时间复杂度不随数据范围的递增而递增,在此给出 std 耗时:
30% <=0.01s
10% <=0.2s
40% <=0.5s
20% <=0.9s

深度:能组成a/b的分数的个数、

剪枝:用没用过的最小的分母,且保证对于当前使用的分母i,剩余分数个数cnt,满足fz/fm<=cnt/i

 1 #include <cstring>
 2 #include <cstdio>
 3
 4 #define LL long long
 5 inline void read(LL &x)
 6 {
 7     x=0; register char ch=getchar();
 8     for(; ch>‘9‘||ch<‘0‘; ) ch=getchar();
 9     for(; ch>=‘0‘&&ch<=‘9‘; ch=getchar()) x=x*10+ch-‘0‘;
10 }
11
12 const int N(1e5+5);
13 LL a,b,n,tmp[N],ans[N];
14 bool flag;
15
16 #define max(a,b) (a>b?a:b)
17 inline int get_first(int a,int b)
18 {
19     for(int i=1; ; i++)
20         if(1.0/i<=a*1.0/b) return i;
21 }
22 inline bool if_change()
23 {
24     for(int i=1; i<=n; ++i)
25         if(tmp[i]!=ans[i]) return ans[i]>tmp[i];
26     return 0;
27 }
28 LL GCD(LL a,LL b)
29 {
30     return !b? a: GCD(b,a%b);
31 }
32 void DFS(LL fz,LL fm,LL now,LL cnt)
33 {
34     if(fz<0) return ;
35     if(fz==1&&fm>=now)
36     {
37         tmp[cnt]=fm;
38         if(if_change()||!flag)
39         {
40             for(int i=n; i>=1; --i)
41                 ans[i]=tmp[i];
42         }
43         flag=1;
44         return ;
45     }
46     if(cnt==1) return ;
47     for(LL gcd,i=max(now,get_first(fz,fm)); (fz*i)<=(fm*cnt); ++i)
48     {
49         tmp[cnt]=i;
50         gcd=GCD(fz*i-fm,fm*i);
51         DFS((fz*i-fm)/gcd,fm*i/gcd,i+1,cnt-1);
52     }
53 }
54
55 int Persist()
56 {
57     freopen("egypt.in","r",stdin);
58     freopen("egypt.out","w",stdout);
59
60     read(a),read(b);
61     ans[1]=1<<30;
62     for( n=1; !flag; ++n)
63         DFS(a,b,get_first(a,b),n);
64     for( ; --n; ) printf("%I64d ",ans[n]);
65     return 0;
66 }
67
68 int Aptal=Persist();
69 int main(){;}

迭代加深+剪枝

 1 #include <cstdio>
 2
 3 #define LL long long
 4 inline void read(LL &x)
 5 {
 6     x=0; register char ch=getchar();
 7     for(; ch>‘9‘||ch<‘0‘; ) ch=getchar();
 8     for(; ch>=‘0‘&&ch<=‘9‘; ch=getchar()) x=x*10+ch-‘0‘;
 9 }
10
11 const int N(1e5+5);
12 LL a,b,n,tmp[N],ans[N];
13 bool flag;
14
15 LL GCD(LL a,LL b)
16 {
17     return !b? a: GCD(b,a%b);
18 }
19 void DFS(LL fz,LL fm,LL now,LL cnt)
20 {
21     if(fz<0) return ;
22     if(cnt==1)
23     {
24         if(fz==1&&fm>=now)
25             if(fm<ans[1])
26             {
27                 flag=1;
28                 ans[1]=fm;
29                 for(int i=n; i>1; --i)
30                     ans[i]=tmp[i];
31             }
32         return ;
33     }
34     for(LL gcd,i=now; (fz*i)<=(fm*cnt); ++i)
35     {
36         tmp[cnt]=i;
37         gcd=GCD(fz*i-fm,fm*i);
38         DFS((fz*i/gcd-fm/gcd),fm*i/gcd,i+1,cnt-1);
39     }
40 }
41
42 int Persist()
43 {
44     read(a),read(b);
45     ans[1]=1<<30;
46     for( n=1; ; ++n)
47     {
48         DFS(a,b,2,n);
49         if(flag) break;
50     }
51     for( ; n; n--) printf("%lld ",ans[n]);
52     return 0;
53 }
54
55 int Aptal=Persist();
56 int main(){;}

CODEVS AC

T3走楼梯  (COGS 2752. [济南集训 2017] 数列运算

(stair.cpp/c/pas)
Description
由于经典走楼梯问题被 xxy 升级了 2 次,上帝非常满意,所以当 xxy 完成这最
后一次跳跃时,会有一定的几率开启轮回之梯。轮回之梯有 n 阶,每一阶上都会有一个
轮回指数,xxy 用不同的策略走到轮回之梯的第 n 阶,会得到不同的轮回值。
由于轮回之梯不是普通的楼梯,每走一步需要消耗极大的体力,xxy 体力有限。上
帝允许 xxy 在中途休息。所以 xxy 如果在第 i 阶,他可以选择立刻走到第 i+1 阶,也可
以选择在第 i 阶休息一会儿。休息和立刻走看做不同的走楼梯方式。
上帝说:我现在给你提供所有台阶的轮回指数,我需要你回答所有不同的走轮回之
梯的方式获得的轮回值之和。如果某种方式你算不出来,你可以亲自去走一遍。你走的次数
越少,就会有越大的几率开启六道轮回。Xxy 非常想开启六道轮回,但由于之前在楼梯上跳
来非常累,现在不想动弹,连飞也不想飞。所以只能在求助学信息学奥赛的你了。
轮回值计算方式:
轮回值为所有的轮回子值之和。
设第 i 个台阶的轮回指数为 xi,如果 xxy 连续从第 L 阶走到第 R 阶,那么
xxy 可以得到的轮回子值为 xL 一直累乘到 xR。特别的,当 L=R 时,轮回子值为 xL。
注意:xxy 开始在轮回之梯的第 1 阶,可以选择先休息一会儿,也可以立刻走到
第 2 阶。走一遍轮回之梯指的是从第 1 阶走到第 n 阶。
由于答案很大,上帝可不想看到几千几百个数字,所以你只需要告诉他答案对
1000000007 取模。
Input
第一行一个整数 n。
接下来 n 个整数,表示 xi
Output
输出一行表示答案
Example
stair.in stair.out
2 30
1 2 4
Hint
对于 10%的数据,1<=n<=10,1<=xi<=10
对于另外 20%的数据,所有的 Xi=1
对于另外 20%的数据,n<=100,xi<=1e5
对于另外 20%的数据,n<=1000,xi<=1e5
对于 100%的数据,n<=100000,xi<=1e9

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<iostream>
 4 #define mod 1000000007
 5 using namespace std;
 6 int n,a[100001];
 7 long long cf[100001];
 8 long long dp[100001],pre[100001],inv[100001];
 9 long long pow(long long a,long long b)
10 {
11     long long res=1;
12     while(b)
13     {
14         if(b&1) res*=a,res%=mod;
15         b>>=1; a*=a; a%=mod;
16     }
17     return res;
18 }
19 void read(int &x)
20 {
21     x=0; char c=getchar();
22     while(!isdigit(c)) c=getchar();
23     while(isdigit(c)) { x=x*10+c-‘0‘; c=getchar(); }
24 }
25 int main()
26 {
27     read(n);
28     for(int i=1;i<=n;i++)   read(a[i]);
29     cf[0]=1;
30     for(int i=1;i<=n;i++) cf[i]=cf[i-1]*2%mod;
31     long long tot; pre[0]=1;
32     for(int i=1;i<=n;i++) pre[i]=pre[i-1]*a[i]%mod;
33     for(int i=1;i<=n;i++) inv[i]=pow(pre[i],mod-2);
34     long long sum1=0,sum2=1;
35     for(int i=1;i<=n;i++)
36     {
37         dp[i]=(sum1+sum2*pre[i]%mod)%mod;
38         sum1=(sum1+dp[i])%mod;
39         sum2=(sum2+cf[i-1]*inv[i]%mod)%mod;
40     }
41     printf("%I64d",dp[n]);
42 }

qiliao

时间: 2024-10-31 18:57:49

9.13——模拟赛的相关文章

2017/9/13模拟赛

粉饰(decorate) [题目描述] 小D有一块被分为n*m个格子的矩形鱼片.为了装饰鱼片,小D决定给每个格子上色.由于小D很喜欢红白,所以小D给每个格子涂上了红色或白色,第i行第j列的格子颜色记为c[i,j].涂完之后,小D想评估这块鱼片的"XY值".我们定义一个有序无重复三元格子组{(x1,y1),(x2,y2),(x3,y3)}为"XY组"当且仅当: |(x1-x2)*(y1-y2)|+|(x3-x2)*(y3-y2)|=0   (c[x1,y1]-c[x2

9.13模拟赛

全排列 (permutation.cpp/c/pas)Description从 n 个不同元素中任取 m(m≤n)个元素,按照一定的顺序排列起来,叫做从 n个不同元素中取出 m 个元素的一个排列.当 m=n 时所有的排列情况叫全排列.你觉得 xxy 会问你全排列的个数吗?Xxy:这个问题能淹死你,我才不问呢.我要问的是求 n 的全排列中,先递增后递减.先递减后递增的全排列的个数.由于答案可能很大,对 p 取余Input输入包含多组测试数据每组测试数据一行两个整数 n,pOutput对于每组测试数

2014.12.13模拟赛【准考证号】

准考证号 128M 0.1s ticket.cpp escription 蒟蒻hzwer NOIP2014惨跪,他依稀记得他的准考证号是37,现在hzwer又将要面临一场比赛,他希望准考证号不出现37(连续),同时他又十分讨厌4,所以也希望4出现在准考证号中...现在他想知道在A和B之间有多少合法的准考证号 Input 包含两个整数,A B Output 一个整数. Sample Input [输入样例一]1 10[输入样例二]25 50 Sample Output [输出样例一]9[输出样例二

8.13模拟赛 T2

题面 大意:10000组询问,每次查询Σ(0<=i<=n)中(i! 末尾0的个数为偶数的 i 的个数)(n<=100000) 0的个数也就是求5的因子数  但我只会打暴力,这里我参考了xumingyang大佬的做法,用了类数位DP的方法 首先把n转化为5进制  得到这样子的一个数 如 1231424,可知这个数为 5?6*1+5?5*2+5?4*3+5?3*1+5?2*4+5?1*2+5?0*4 f[i][0/1]表示到第i位设一个数有n位 , 原文地址:https://www.cnbl

【BZOJ】【2741】【FOTILE模拟赛】L

可持久化Trie+分块 神题……Orz zyf & lyd 首先我们先将整个序列搞个前缀异或和,那么某一段的异或和,就变成了两个数的异或和,所以我们就将询问[某个区间中最大的区间异或和]改变成[某个区间中 max(两个数的异或和)] 要是我们能将所有[l,r]的答案都预处理出来,那么我们就可以O(1)回答了:然而我们并不能. 一个常见的折中方案:分块! 这里先假设我们实现了一个神奇的函数ask(l,r,x),可以帮我们求出[l,r]这个区间中的数,与x最大的异或值. 我们不预处理所有的左端点,我

10.30 NFLS-NOIP模拟赛 解题报告

总结:今天去了NOIP模拟赛,其实是几道USACO的经典的题目,第一题和最后一题都有思路,第二题是我一开始写了个spfa,写了一半中途发现应该是矩阵乘法,然后没做完,然后就没有然后了!第二题的暴力都没码QAQ 现在我来写解题报告了,有点饿了QAQ.. 第一题 题目 1: 架设电话线 [Jeffrey Wang, 2007] 最近,Farmer John的奶牛们越来越不满于牛棚里一塌糊涂的电话服务,于 是,她们要求FJ把那些老旧的电话线换成性能更好的新电话线.新的电话线架设 在已有的N(2 <=

【题解】PAT团体程序设计天梯赛 - 模拟赛

由于本人愚笨,最后一题实在无力AC,于是只有前14题的题解Orz 总的来说,这次模拟赛的题目不算难,前14题基本上一眼就有思路,但是某些题写起来确实不太容易,编码复杂度有点高~ L1-1 N个数求和 设计一个分数类,重载加法运算符,注意要约分,用欧几里得算法求个最大公约数即可. 1 #include <cstdio> 2 3 long long abs(long long x) 4 { 5 return x < 0 ? -x : x; 6 } 7 8 long long gcd(long

[GRYZ]寒假模拟赛

写在前面 这是首次广饶一中的OIERS自编自导,自出自做(zuo)的模拟赛. 鉴于水平气压比较低,机(wei)智(suo)的WMY/XYD/HYXZC就上网FQ下海找了不少水(fei)题,经过他们优(le)美(se)的文字加工后,有故事有题目有人物有奸情的模拟赛正式呈上. 我是正文 题目名 GRYZ娱乐时刻 GRYZ追击时刻 GRYZ就餐时刻 源文件 hyxzc.cpp/c/pas clikar.cpp/c/pas eat.cpp/c/pas 输入文件 hyxzc.in clikar.in ea

[铁一中OI模拟赛]2017.8.19 Day1

T1 小Z的情书 题目链接 思考: 题目主要难度在于旋转后的位置如何,在手写了样例之后不难发现规律. #include <cstdio> #include <cstring> #define up(a,b,c) for(register int c=a;c<=b;++c) #define down(a,b,c) for(register int c=a;c>=b;--c) const int Maxn=1005; int n; bool Map[Maxn][Maxn],