HDU-5755-Gambler Bo-高斯消元
2016-07-28 11:49
447 查看
题意:高斯消元
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
int a[905][905];
int n,m;
int maxn;
int MOD=3;
int ans[811111];
int loc(int x,int y) {
return x*m-m+y;
}
void swapl(int x,int y) {
int tt;
for(int i=1;i<=m+1;i++) {
tt=a[x][i];
a[x][i]=a[y][i];
a[y][i]=tt;
}
return;
}
inline int gcd(int a,int b) {
return b==0?a:gcd(b,a%b);
}
inline int lcm(int a,int b) {
return a/gcd(a,b)*b;
}
long long inv(long long a,long long m) {
if(a==1) return 1;
return inv(m%a,m)*(m-m/a)%m;
}
int Guass(int equ,int var) {
int max_r,col,k;
for(k=0,col=0;k<equ&&col<var;k++,col++) {
max_r=k;
for(int i=k+1;i<equ;i++)
if(abs(a[i][col])>abs(a[max_r][col]))
max_r=i;
if(a[max_r][col]==0) {
k--;
continue;
}
if(max_r!=k) {
for(int j=col;j<var+1;j++)
swap(a[k][j],a[max_r][j]);
}
for(int i=k+1;i<equ;i++) {
if(a[i][col]!=0) {
int LCM=lcm(abs(a[i][col]),abs(a[k][col]));
int ta=LCM/abs(a[i][col]);
int tb=LCM/abs(a[k][col]);
if(a[i][col]*a[k][col]<0) tb=-tb;
for(int j=col;j<var+1;j++) {
a[i][j]=((a[i][j]*ta-a[k][j]*tb)%MOD+MOD)%MOD;
}
}
}
}
for(int i=var-1;i>=0;i--) {
int temp=a[i][var];
for(int j=i+1;j<var;j++) {
if(a[i][j]!=0) {
temp-=a[i][j]*ans[j];
temp=(temp%MOD+MOD)%MOD;
}
}
ans[i]=a[i][i]*temp%MOD;
}
return 0;
}
void prit(void) {
for(int i=0;i<maxn;i++) {
for(int j=0;j<maxn+1;j++)
cout << a[i][j] << " ";
cout << endl;
}
}
int main()
{
int t;
int kase=0;
int xx;
scanf("%d",&t);
while(t--){
scanf("%d%d",&n,&m);
memset(ans,0,sizeof ans);
memset(a,0,sizeof a);
maxn=n*m;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++) {
scanf("%d",&xx);
a[loc(i,j)-1][maxn+1-1]=(3-xx)%3;
a[loc(i,j)-1][loc(i,j)-1]=2;
if(i>1) a[loc(i,j)-1][loc(i-1,j)-1]=1;
if(j>1) a[loc(i,j)-1][loc(i,j-1)-1]=1;
if(i<n) a[loc(i,j)-1][loc(i+1,j)-1]=1;
if(j<m) a[loc(i,j)-1][loc(i,j+1)-1]=1;
}
// prit();
Guass(maxn,maxn);
// cout << endl;
// prit();
int cnt = 0;
for(int i = 0; i < maxn; i++) {
if(ans[i])
cnt += ans[i];
}
printf("%d\n", cnt);
for(int i = 0; i < maxn; i++) {
while(ans[i]) {
printf("%d %d\n", i/m+1, i%m+1);
ans[i]--;
}
}
//printf("%d\n",tot);
}
return 0;
}
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
int a[905][905];
int n,m;
int maxn;
int MOD=3;
int ans[811111];
int loc(int x,int y) {
return x*m-m+y;
}
void swapl(int x,int y) {
int tt;
for(int i=1;i<=m+1;i++) {
tt=a[x][i];
a[x][i]=a[y][i];
a[y][i]=tt;
}
return;
}
inline int gcd(int a,int b) {
return b==0?a:gcd(b,a%b);
}
inline int lcm(int a,int b) {
return a/gcd(a,b)*b;
}
long long inv(long long a,long long m) {
if(a==1) return 1;
return inv(m%a,m)*(m-m/a)%m;
}
int Guass(int equ,int var) {
int max_r,col,k;
for(k=0,col=0;k<equ&&col<var;k++,col++) {
max_r=k;
for(int i=k+1;i<equ;i++)
if(abs(a[i][col])>abs(a[max_r][col]))
max_r=i;
if(a[max_r][col]==0) {
k--;
continue;
}
if(max_r!=k) {
for(int j=col;j<var+1;j++)
swap(a[k][j],a[max_r][j]);
}
for(int i=k+1;i<equ;i++) {
if(a[i][col]!=0) {
int LCM=lcm(abs(a[i][col]),abs(a[k][col]));
int ta=LCM/abs(a[i][col]);
int tb=LCM/abs(a[k][col]);
if(a[i][col]*a[k][col]<0) tb=-tb;
for(int j=col;j<var+1;j++) {
a[i][j]=((a[i][j]*ta-a[k][j]*tb)%MOD+MOD)%MOD;
}
}
}
}
for(int i=var-1;i>=0;i--) {
int temp=a[i][var];
for(int j=i+1;j<var;j++) {
if(a[i][j]!=0) {
temp-=a[i][j]*ans[j];
temp=(temp%MOD+MOD)%MOD;
}
}
ans[i]=a[i][i]*temp%MOD;
}
return 0;
}
void prit(void) {
for(int i=0;i<maxn;i++) {
for(int j=0;j<maxn+1;j++)
cout << a[i][j] << " ";
cout << endl;
}
}
int main()
{
int t;
int kase=0;
int xx;
scanf("%d",&t);
while(t--){
scanf("%d%d",&n,&m);
memset(ans,0,sizeof ans);
memset(a,0,sizeof a);
maxn=n*m;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++) {
scanf("%d",&xx);
a[loc(i,j)-1][maxn+1-1]=(3-xx)%3;
a[loc(i,j)-1][loc(i,j)-1]=2;
if(i>1) a[loc(i,j)-1][loc(i-1,j)-1]=1;
if(j>1) a[loc(i,j)-1][loc(i,j-1)-1]=1;
if(i<n) a[loc(i,j)-1][loc(i+1,j)-1]=1;
if(j<m) a[loc(i,j)-1][loc(i,j+1)-1]=1;
}
// prit();
Guass(maxn,maxn);
// cout << endl;
// prit();
int cnt = 0;
for(int i = 0; i < maxn; i++) {
if(ans[i])
cnt += ans[i];
}
printf("%d\n", cnt);
for(int i = 0; i < maxn; i++) {
while(ans[i]) {
printf("%d %d\n", i/m+1, i%m+1);
ans[i]--;
}
}
//printf("%d\n",tot);
}
return 0;
}
相关文章推荐
- 使用 Spring 的 AOP 机制来输出 Log
- Solidworks安装界面是英文的解决办法
- 行为识别笔记:Stacked Fisher Vector基本原理
- php学习1-----制作PHP留言板:
- 题目25 A Famous Music Composer
- Sublime ctags 函数跳转插件安装
- Murano Weekly Meeting 2016.07.05
- Django session
- Kotlin Reference(二)-基础部分
- 【NYOJ58】最小步数
- 别猜白块滑动效果实现
- Studio中的项目接入百度地图sdk后,地图不能正常显示及出现 errorcode: 230 uid: -1 appid -1 msg: APP Scode码校验失败 问题的解决方法
- 谈 UIView Animation 编程艺术
- CSS学习笔记--基础(一)
- 关于mysql的join(图)
- 70.打印所有Spring boot载入的bean【从零开始学Spring Boot】
- APP测试要点
- 【OpenJudge】c语言_计算2的幂
- 80个Python经典资料(教程+源码+工具)汇总
- mysql 的一些基本操作