题目链接:Construct a Matrix

题意:构造一个矩阵,要求矩阵的每行每列的和都不相同。矩阵的边长是前n项斐波那契的和。

思路:由sn = 2*(fn-1)+(fn-2)-1,只要知道第n-1和第n-2项即可,n的范围是10^9,可由矩阵快速幂求出第n项。然后,构造矩阵,上三角为1,下三角全为-1,对角线1和0交替。【真是个天才...!!!】矩阵快速幂求第n项时,构造的矩阵是a[0][0] = f2, a[1][0] = f1, a[0][1] = 1, a[1][1] = 0........【ACMer的脑洞真大...可怕....当然不包括我辣...】

知道思路当然就好实现了。附代码:

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
using namespace std; struct Mat{
int a[2][2];
void init1() { // 斐波那契初始矩阵
a[0][0] = a[0][1] = a[1][0] = 1;
a[1][1] = 0;
}
void init2() { //单位矩阵
a[0][0] = a[1][1] = 1;
a[0][1] = a[1][0] = 0;
}
}; Mat mul(Mat aa, Mat b, int m) { // 矩阵乘法mod m
Mat ans;
for (int i=0; i<2; ++i) {
for (int j=0; j<2; ++j) {
ans.a[i][j] = 0;
for (int k=0; k<2; ++k) {
ans.a[i][j] += ((aa.a[i][k]%m)*(b.a[k][j]%m))%m;
ans.a[i][j] %= m;
//ans.a[i][j] += aa.a[i][k]*b.a[k][j];
}
//cout << ans.a[i][j] << "....\n";
}
}
return ans;
} Mat mul_mat_quick(Mat a, int n, int m) { // 矩阵快速幂%m
Mat ans;
ans.init2();
while(n) {
if (n%2) ans = mul(ans, a, m);
a = mul(a, a, m);
n /= 2;
}
return ans;
} int num[210][210]; int main() {
int t;
int cas = 0;
scanf("%d", &t);
while(t--) {
int n, m;
scanf("%d%d", &n, &m);
Mat a;
a.init1();
Mat ans = mul_mat_quick(a, n-1, m);
int fn1 = ans.a[0][0]; // fn-1
int fn2 = ans.a[1][0]; // fn-2
//cout << fn1 << "++++=" << fn2 << endl;
int sn = (2*fn1 + fn2 - 1)%m; //矩阵边长
// cout << sn << "===\n"; if (sn%2 || sn==0) {
printf("Case %d: No\n", ++cas);
continue;
} printf("Case %d: Yes\n", ++cas);
memset(num, 0, sizeof(num));
for (int i=1; i<=sn; ++i) {
for (int j=1; j<=sn; ++j) {
if (i<j) num[i][j] = 1;
else num[i][j] = -1;
}
}
int pre = 1;
for (int i=1; i<=sn; ++i) {
num[i][i] = (pre^1);
pre ^= 1;
}
for (int i=1; i<=sn; ++i) {
for (int j=1; j<=sn; ++j) {
if (j==1) printf("%d", num[i][j]);
else printf(" %d", num[i][j]);
}
printf("\n");
}
}
return 0;
}

  

FZU 1911 Construct a Matrix的更多相关文章

  1. fzu 1911 Construct a Matrix(矩阵快速幂+规律)

    题目链接:fzu 1911 Construct a Matrix 题目大意:给出n和m,f[i]为斐波那契数列,s[i]为斐波那契数列前i项的和.r = s[n] % m.构造一个r * r的矩阵,只 ...

  2. fzu 1911 C. Construct a Matrix

    C. Construct a Matrix Time Limit: 1000ms Case Time Limit: 1000ms Memory Limit: 32768KB Special Judge ...

  3. Construct a Matrix (矩阵快速幂+构造)

    There is a set of matrixes that are constructed subject to the following constraints: 1. The matrix ...

  4. KUANGBIN带你飞

    KUANGBIN带你飞 全专题整理 https://www.cnblogs.com/slzk/articles/7402292.html 专题一 简单搜索 POJ 1321 棋盘问题    //201 ...

  5. [kuangbin带你飞]专题1-23题目清单总结

    [kuangbin带你飞]专题1-23 专题一 简单搜索 POJ 1321 棋盘问题POJ 2251 Dungeon MasterPOJ 3278 Catch That CowPOJ 3279 Fli ...

  6. ACM--[kuangbin带你飞]--专题1-23

    专题一 简单搜索 POJ 1321 棋盘问题POJ 2251 Dungeon MasterPOJ 3278 Catch That CowPOJ 3279 FliptilePOJ 1426 Find T ...

  7. <转载> OpenGL Projection Matrix

    原文 OpenGL Projection Matrix Related Topics: OpenGL Transformation Overview Perspective Projection Or ...

  8. Palindromic Matrix

    Palindromic Matrix time limit per test 2 seconds memory limit per test 256 megabytes input standard ...

  9. Codeforces Round #540 (Div. 3) C. Palindromic Matrix 【暴力】

    任意门:http://codeforces.com/contest/1118/problem/C C. Palindromic Matrix time limit per test 2 seconds ...

随机推荐

  1. POJ3009 Curling

    题目链接:http://poj.org/problem?id=3009 题意:从2出发,要到达3, 0可以通过,碰到1要停止,并且1处要变成0, 并且从起点开始沿着一个方向要一直前进,直至碰到1(或者 ...

  2. 树--四分树(UVa297)

    郑重声明: 数据结构这部分内容, 由于博主才学很少(且很浅)的内容, 所以现在所写的(大都是抄的)一些典型例题, 再加上一些自己想法和理解而已, 等博主勤加修炼, 以后会大有补充和改进. 粗浅之处, ...

  3. 容易混淆的url src href

    新手刚学习的时候会分不清 url  src  href这些,不知道什么情况下应该用哪个.现在让我来理一理. url 统一资源定位符是对可以从互联网上得到的资源的位置和访问方法的一种简洁的表示,是互联网 ...

  4. jQuery 中$(this).parent().parent().remove()无效。

    在写文章系统的删除功能.需要删除一行数据.在删除的页面,需要jQuery 删除一hang. 局部刷新数据. $(".del").click(function(){ var id = ...

  5. VNC的安装与配置

    一,安装tigervnc-server VNC软件包 [root@localhost ~]# yum install tigervnc-server 设置开机自启动 [root@localhost ~ ...

  6. [Selenium]点击Calendar控件后,Calendar dialog很快消失

    有的日历控件使用了“opacity:0”透明度加上“display:none”来控制其是否显示,使用moveToElement方法不能移动到dialog上,可以把这两个属性的值进行修改,使其可见,使用 ...

  7. Druid 数据库用户密码加密 代码实现

    druid-1.0.16.jar 阿里巴巴的开源数据连接池 jar包 明文密码+私钥(privateKey)加密=加密密码 加密密码+公钥(publicKey)解密=明文密码 程序代码如下: pack ...

  8. mysql使用笔记(四)

    一.选择合适的数据类型 1. CHAR vs VCHAR     char是固定长度的字符类型,而varchar是可变长度的字符类型.char(M)的数据列中,每个值都占用M个字节,如果某个长度小于M ...

  9. Http简析

    HTTP协议 属于应用层的面向对象的协议 HTTP协议的主要特点 支持C/S(客户/服务器)模式. 简单快速:客户向服务器请求服务时,只需传送请求方法和路径.请求方法常用的有GET.HEAD.POST ...

  10. maven和jdk版本不匹配

    解决方法:http://blog.csdn.net/mafan121/article/details/51944346