http://acm.nyist.net/JudgeOnline/problem.php?pid=301递推求值

2024-01-10 07:48

本文主要是介绍http://acm.nyist.net/JudgeOnline/problem.php?pid=301递推求值,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

矩阵运算。。。这一题让我明白了一些事,在做题的时候一定要考虑数的取值范围。。。否者会多吃WA的这一题我就是因为这wa了好几次。。。

#include<iostream>
#include<string.h>
#include<cstdio>
#include<algorithm>
#define M 1000007
typedef long long L;
typedef struct
{   L s[3][3];
}Node;
Node a,b;
Node ceil(Node p,Node q)
{   Node c;memset(c.s,0,sizeof(c.s));for(int i=0;i<3;++i)for(int j=0;j<3;++j)for(int t=0;t<3;++t)c.s[i][j]=((c.s[i][j]+p.s[i][t]*q.s[t][j])%M+M)%M;return c;
}
Node doit(int k)
{   Node p=b,q=a;while(k){  if(k&1) p=ceil(p,q);q=ceil(q,q);k=k>>1;}return p;
}
int main()
{  int T;scanf("%d",&T);while(T--){  int f1,f2,a1,b1,c1,n;scanf("%d%d%d%d%d%d",&f1,&f2,&a1,&b1,&c1,&n);if(n==1) printf("%d\n",f1);else if(n==2) printf("%d\n",f2);else if(n>=3){  a.s[0][0]=0;a.s[0][1]=a1;a.s[0][2]=0;a.s[1][0]=1;a.s[1][1]=b1;a.s[1][2]=0;a.s[2][0]=0;a.s[2][1]=1;a.s[2][2]=1;b.s[0][0]=1;b.s[0][1]=0;b.s[0][2]=0;b.s[1][0]=0; b.s[1][1]=1;b.s[1][2]=0;b.s[2][0]=0;b.s[2][1]=0;b.s[2][2]=1;n-=2;Node c=doit(n);L k=((f1*c.s[0][1]+f2*c.s[1][1]+c1*c.s[2][1])%M+M)%M;printf("%lld\n",k);}}return 0;
}

一般一个题这样算下来不超过1亿的话,一秒左右都能过;


嗯,大致数量级就是亿级别
十亿就很可能会挂
不过,这都只是一般情况
;
编译器不同,机器配置不同,操作系统不同;

这篇关于http://acm.nyist.net/JudgeOnline/problem.php?pid=301递推求值的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/590073

相关文章

认识、理解、分类——acm之搜索

普通搜索方法有两种:1、广度优先搜索;2、深度优先搜索; 更多搜索方法: 3、双向广度优先搜索; 4、启发式搜索(包括A*算法等); 搜索通常会用到的知识点:状态压缩(位压缩,利用hash思想压缩)。

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n

uva 568 Just the Facts(n!打表递推)

题意是求n!的末尾第一个不为0的数字。 不用大数,特别的处理。 代码: #include <stdio.h>const int maxn = 10000 + 1;int f[maxn];int main(){#ifdef LOCALfreopen("in.txt", "r", stdin);#endif // LOCALf[0] = 1;for (int i = 1; i <=

uva 10025 The ? 1 ? 2 ? ... ? n = k problem(数学)

题意是    ?  1  ?  2  ?  ...  ?  n = k 式子中给k,? 处可以填 + 也可以填 - ,问最小满足条件的n。 e.g k = 12  - 1 + 2 + 3 + 4 + 5 + 6 - 7 = 12 with n = 7。 先给证明,令 S(n) = 1 + 2 + 3 + 4 + 5 + .... + n 暴搜n,搜出当 S(n) >=

如何在Visual Studio中调试.NET源码

今天偶然在看别人代码时,发现在他的代码里使用了Any判断List<T>是否为空。 我一般的做法是先判断是否为null,再判断Count。 看了一下Count的源码如下: 1 [__DynamicallyInvokable]2 public int Count3 {4 [__DynamicallyInvokable]5 get

2、PF-Net点云补全

2、PF-Net 点云补全 PF-Net论文链接:PF-Net PF-Net (Point Fractal Network for 3D Point Cloud Completion)是一种专门为三维点云补全设计的深度学习模型。点云补全实际上和图片补全是一个逻辑,都是采用GAN模型的思想来进行补全,在图片补全中,将部分像素点删除并且标记,然后卷积特征提取预测、判别器判别,来训练模型,生成的像

BUUCTF靶场[web][极客大挑战 2019]Http、[HCTF 2018]admin

目录   [web][极客大挑战 2019]Http 考点:Referer协议、UA协议、X-Forwarded-For协议 [web][HCTF 2018]admin 考点:弱密码字典爆破 四种方法:   [web][极客大挑战 2019]Http 考点:Referer协议、UA协议、X-Forwarded-For协议 访问环境 老规矩,我们先查看源代码

【Linux】应用层http协议

一、HTTP协议 1.1 简要介绍一下HTTP        我们在网络的应用层中可以自己定义协议,但是,已经有大佬定义了一些现成的,非常好用的应用层协议,供我们直接使用,HTTP(超文本传输协议)就是其中之一。        在互联网世界中,HTTP(超文本传输协议)是一个至关重要的协议,他定义了客户端(如浏览器)与服务器之间如何进行通信,以交换或者传输超文本(比如HTML文档)。

如何确定 Go 语言中 HTTP 连接池的最佳参数?

确定 Go 语言中 HTTP 连接池的最佳参数可以通过以下几种方式: 一、分析应用场景和需求 并发请求量: 确定应用程序在特定时间段内可能同时发起的 HTTP 请求数量。如果并发请求量很高,需要设置较大的连接池参数以满足需求。例如,对于一个高并发的 Web 服务,可能同时有数百个请求在处理,此时需要较大的连接池大小。可以通过压力测试工具模拟高并发场景,观察系统在不同并发请求下的性能表现,从而