[HNOI2010]MATRIX 矩阵
Description
Input
第一行包含三个正整数N M P表示矩阵的行数列数以及每个数的范围,接下来N行每行包含M个非负整数,其中第i行第j个数表示以格子(i,j)为右下角的2*2子矩阵中的数的和。保证第一行与第一列的数均为0,且每个和都不超过4(P-1)。
Output
包含N行,每行M个整数,描述你求出的矩阵,相邻的整数用空格分开。(行末不要有多余空格)
Sample Input
0 0 0
0 4 5
0 5 3
Sample Output
2 2 1
1 0 0
HINT
1<=N,M<=200
1<P<=10
首先可以推出,一旦A的第一行与第一列已经确定,那么我们就可以按部就班地
算出A的所有元素的值了。这样一来需要决策的元素个数就从$200^{2}$变成了399
首先不看元素在0到P-1范围内的限制,设A(1,k)和A(k,1)都为0,并计算将A补全。
(此时A中可能有负数或者大于等于P的整数)
这样A中有些值不符合条件,调整A矩阵
当我们改变A(1,1)+=k时,
要把A调整成合法矩阵,事实上只需把i+j为偶数的A(i,j)
全部增加 k,i+j为奇数的 A(i,j)全部减去 k 即可
当我们把A(1,i)增加k时,要把A调整成合法,只需要把第i列奇数位加k,偶数位减k即可。
我们只需要搜索A的第一行,每当确定一个元素,就可以更新A(j,1)的取值范围。
根据公式,除了A(1,1)以外,A(i,1)的取值是互不影响的。
不断更新第一列每个元素的可以取值的范围,一旦出现有某个元素下界大于上界的情况就剪枝
显然这样搜索的字典序是最优的
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
int l[][],r[][],a[][],s[][],c[][],n,m,p;
int pd(int x)
{
if (x&) return -;
return ;
}
int get_num(int x,int y)
{
return c[x][y]-pd(x+y)*a[][]-pd(x)*a[][y]-pd(y)*a[x][];
}
bool dfs(int x)
{int k,i,j,ok;
if (x>m) return ;
for (k=;k<p;k++)
{
ok=;
a[][x]=k;
for (i=;i<=n;i++)
{
int rl=(c[i][x]-pd(x+i)*a[][]-pd(i)*a[][x])*(pd(x));
int rr=(c[i][x]-pd(x+i)*a[][]-pd(i)*a[][x]-(p-))*(pd(x));
if (rl>rr) swap(rl,rr);
l[i][x]=max(l[i][x-],rl);
r[i][x]=min(r[i][x-],rr);
if (l[i][x]>r[i][x])
{
ok=;
break;
}
}
if (ok&&dfs(x+)) return ;
}
return ;
}
int main()
{int i,j,k;
cin>>n>>m>>p;
for (i=;i<=n;i++)
{
for (j=;j<=m;j++)
{
scanf("%d",&s[i][j]);
l[i][j]=;r[i][j]=p-;
}
}
for (i=;i<=n;i++)
{
for (j=;j<=m;j++)
{
c[i][j]=s[i][j]-c[i-][j]-c[i][j-]-c[i-][j-];
}
}
for (i=;i<p;i++)
{
a[][]=i;
if (dfs())
{
for (j=;j<=n;j++)
{
a[j][]=l[j][m];
}
for (j=;j<=n;j++)
{
for (k=;k<=m;k++)
a[j][k]=get_num(j,k);
}
break;
}
}
for (i=;i<=n;i++)
{
for (j=;j<m;j++)
{
printf("%d ",a[i][j]);
}
printf("%d",a[i][m]);
printf("\n");
}
}
首先可以推出,一旦 A 的第一行与第一列已经确定,那么我们就可以按部就班地
算出 A 的所有元素的值了。这样一来需要决策的元素个数就从 200 2 变成了 399。
为了方便计算,我们试图建立 A i,j 与第一行或第一列之间的直接等式关系。
以下是一个较易理解的建立等量关系的方法:
首先不看元素在 0 到 P-1 范围内的限制,设 A 1,k 和 A k,1 都为 0,并计算将 A 补全。
(此时 A 中可能有负数或者大于等于 P 的整数)
S:
0 0 0
0 4 5
0 5 3A:
0 0 0
0 4 1
0 1 -3
当我们改变 A 1,1 +=k 时,
要把 A 调整成合法矩阵,事实上只需把 i+j 为偶数的 A i,j
全部增加 k,i+j 为奇数的 A i,j 全部减去 k 即可。
A’: A 1,1 +=1
1 -1 1
-1 5 0
1 0 -2
当我们把 A 1,i 增加 k 时,要把 A 调整成合法,只需要把第 i 列奇数位加 k,偶数位
减 k 即可。
A’: A 1,2 +=1
0 1 0
0 3 1
0 2 -3
修改 A i,1 类似。
设通过令第一行第一列都为 0 得到的矩阵为 C,
那么有公式:
A i,j = C i,j + (-1) i+j-2 A 1,1 + (-1) i-1 A 1,j + (-1) j-1 A i,1
(i>1,j>1)
接下来只需要考虑对第一行第一列进行决策即可。
前 30%的数据可以搜索解决。
另 30%实际上是一个 2-Sat.问题。参见 SGU 307 Cipher。
对于 100%的数据并没有很好的方法。
研究 30%的搜索方法,如果加入较强的剪枝,一般出解是很快的。
我们只需要搜索 A 的第一行,每当确定一个元素,就可以更新 A j,1 的取值范围。
根据公式,除了 A 1,1 以外,A i,1 的取值是互不影响的。
不断更新第一列每个元素的可以取值的范围,一旦出现有某个元素下界大于上界的
情况就剪枝。
显然这样搜索的字典序是最优的。
事实证明这个剪枝效率非常高。
[HNOI2010]MATRIX 矩阵的更多相关文章
- 【bzoj2003】[Hnoi2010]Matrix 矩阵
首先可以知道,如果已知第一行和第一列的数,那么可以很容易的计算出其余的数.进一步的,如果笔算将每个数的表达式写出可以得出如下结论: 第i行的第j个数(i>1,j>1)只与(1,1),(i, ...
- 【CSS3】 理解CSS3 transform中的Matrix(矩阵)
理解CSS3 transform中的Matrix(矩阵) by zhangxinxu from http://www.zhangxinxu.com 本文地址:http://www.zhangxinxu ...
- 理解CSS3 transform中的Matrix(矩阵)
一.哥,我被你吓住了 打架的时候会被块头大的吓住,学习的时候会被奇怪名字吓住(如“拉普拉斯不等式”).这与情感化设计本质一致:界面设计好会让人觉得这个软件好用! 所以,当看到上面“Matrix(矩阵) ...
- HDU5015 233 Matrix(矩阵高速幂)
HDU5015 233 Matrix(矩阵高速幂) 题目链接 题目大意: 给出n∗m矩阵,给出第一行a01, a02, a03 ...a0m (各自是233, 2333, 23333...), 再给定 ...
- 理解CSS3 transform中的Matrix(矩阵)——张鑫旭
by zhangxinxu from http://www.zhangxinxu.com本文地址:http://www.zhangxinxu.com/wordpress/?p=2427 一.哥,我被你 ...
- String数据转Matrix矩阵
String数据转Matrix矩阵 private Matrix String_To_Matrix(string str) { int[] Remove_Num = new int[10]; int ...
- Leetcode 566. Reshape the Matrix 矩阵变形(数组,模拟,矩阵操作)
Leetcode 566. Reshape the Matrix 矩阵变形(数组,模拟,矩阵操作) 题目描述 在MATLAB中,reshape是一个非常有用的函数,它可以将矩阵变为另一种形状且保持数据 ...
- Android中的Matrix(矩阵)
写在前面 看这篇笔记之前先看一下参考文章,这篇笔记没有系统的讲述矩阵和代码的东西,参考文章写的也有错误的地方,要辨证的看. 如何计算矩阵乘法 android matrix 最全方法详解与进阶(完整篇) ...
- 汉高澳大利亚matrix矩阵计算器
我在梦中的超级计算机超级计算机锯,使用大量阵列的cpu记忆,完成并行计算.一个手机制造商由于使用普通机械提供的服务,往往造成停机.是铁道部列车网络售票的事实. 无法使用云服务.上万台计算机并行处理,因 ...
随机推荐
- c语言程序设计第4周编程练习(素数和)
1 素数和(5分) 题目内容: 我们认为2是第一个素数,3是第二个素数,5是第三个素数,依次类推. 现在,给定两个整数n和m,0<n<=m<=200,你的程序要计算第n个素数到第m个 ...
- 20162318 实验二《Java面向对象程序设计》实验报告
北京电子科技学院(BESTI) 实 验 报 告 课程:程序设计与数据结构 班级:1623班 姓名:张泰毓 指导老师:娄老师.王老师 实验日期:2017年4月14日 实验密级:非密级 实验器材:带Lin ...
- The sum of numbers form 0 to n.(20.9.2017)
#include <stdio.h> int main() { int a,b,sum; printf("输入一个数字: "); scanf("%d" ...
- Python 实现队列
操作 Queue() 创建一个空的队列 enqueue(item) 往队列中添加一个item元素 dequeue() 从队列头部删除一个元素 is_empty() 判断一个队列是否为空 size() ...
- Python 远程部署 Fabric
参考文章:http://zmrenwu.com/post/21/ Fabric是一个Python的库,它提供了丰富的同SSH交互的接口,可以用来在本地或远程机器上自动化.流水化地执行Shell命令.因 ...
- 如何用tomcat实现类似weblogic那样的热部署方式
平时weblogic部署程序包时一般是到控制台去部署,不需要重启. 相反之前用tomcat部署应用时,我一般都是把tomcat重启来完成程序包的更新或新包部署.但是这次要部署的应用有点多,大概10几个 ...
- java克隆之深拷贝与浅拷贝
版权声明:本文出自汪磊的博客,转载请务必注明出处. Java深拷贝与浅拷贝实际项目中用的不多,但是对于理解Java中值传递,引用传递十分重要,同时个人认为对于理解内存模型也有帮助,况且面试中也是经常问 ...
- ajax中设置contentType: “application/json”的作用
最近在做项目交互的时候,刚开始向后台传递数据返回415,后来百度添加了 contentType:"application/json"之后返回400,然后把传输的数据格式改为json ...
- SpringCloud的应用发布(四)vmvare+linux,防火墙和selinux
一.vmvare网络配置为nat模式 二.vmvare的网络设置为桥接bridge模式 1.linux 网卡的ip获取方式dhcp 三.关闭linux的防火墙和selinux 1.临时关闭防火墙 sy ...
- 【52ABP实战教程】0.1-- Devops如何用VSTS持续集成到Github仓库!
工欲善其事,必先利其器.在开始正式的教程之前我们先来聊聊准备工作. 管理工具会VSTS. 代码管理会用GITHUB. 服务器会用Azure. 所有的东西都是利用现有服务.不会说自己从虚拟机开始玩.我们 ...