USACO Camelot 解题报告

2024-03-29 08:32
文章标签 报告 解题 usaco camelot

本文主要是介绍USACO Camelot 解题报告,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

这道题也是看网上的解题报告得出的解法。。。

首先,knight接king的点不需要在整个棋盘枚举,king最多move两步。这个是能够通过所有测试点的。至于为什么这样,网上也搜到过证明,但是无心看了。

其次,floyd是个非常低效的算法,因为它确定无疑的要n^3的时间复杂度,适合于密集图(因为算法和边的个数无关),但是我们遇到的一般是稀疏图。SPFA甚至BFS都更好。之前看johnson算法似乎是适用于稀疏图的。

再次,usaco给的标程没看懂。。。

/*
ID: thestor1
LANG: C++
TASK: camelot
*/
#include <iostream>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <climits>
#include <vector>
#include <cassert>
#include <string>
#include <algorithm>
#include <stack>
#include <set>
#include <queue>using namespace std;struct Position{int r, c;Position(){}Position(int r, int c){this->r = r;this->c = c;}Position(const Position& other):r(other.r), c(other.c){}Position& operator =(const Position& other){r=other.r;c=other.c;return *this;}
};const int MAXR = 26;
const int MAXC = 30;
const int MAXN = MAXR * MAXC;int dis[MAXR][MAXC][MAXR][MAXC];const int directx[] = {-2, -1, 1, 2, 2, 1, -1, -2};
const int directy[] = {1, 2, 2, 1, -1, -2, -2, -1};bool marked[MAXR][MAXC];
int index(int r, int c, int C)
{return r * C + c;
}bool isWithin(int r, int c, int R, int C)
{return r >= 0 && r < R && c >= 0 && c < C;
}void spfa(const Position &s, const int R, const int C)
{if(marked[s.r][s.c]){return;}marked[s.r][s.c] = true;for(int i = 0; i < R; ++i){for(int j = 0; j < C; ++j){dis[s.r][s.c][i][j] = INT_MAX;}}dis[s.r][s.c][s.r][s.c] = 0;deque<Position> que;bool onQueue[MAXR][MAXC];memset(onQueue, false, MAXN * sizeof(bool));//int s = index(r, c, C);que.push_back(s);onQueue[s.r][s.c] = true;while(!que.empty()){/*if(que.size() == 51 && que.front().r == 24 && que.front().c == 18){fprintf(stdout, "debug\n");}fprintf(stdout, "que.size: %d\n", que.size());*/Position u = que.front();que.pop_front();//fprintf(stdout, "u: (%d, %d)\n", u.r, u.c);if(dis[s.r][s.c][u.r][u.c] == INT_MAX){continue;}onQueue[u.r][u.c] = false;for(int d = 0; d < 8; ++d){int vr = u.r + directx[d];int vc = u.c + directy[d];if(!isWithin(vr, vc, R, C)){continue;}Position v(vr, vc);//v.r = vr;//v.c = vc;if(dis[s.r][s.c][u.r][u.c] + 1 < dis[s.r][s.c][v.r][v.c]){dis[s.r][s.c][v.r][v.c] = dis[s.r][s.c][u.r][u.c] + 1;if(!onQueue[v.r][v.c]){if(que.empty() || dis[s.r][s.c][v.r][v.c] < dis[s.r][s.c][que.front().r][que.front().c]){que.push_front(v);}else{que.push_back(v);}onQueue[v.r][v.c] = true;}}}}
}int main()
{FILE *fin  = fopen ("camelot.in", "r");FILE *fout = fopen ("camelot.out", "w");//freopen("log.txt", "w", stdout);int R, C;fscanf(fin, "%d%d\n", &C, &R);//int N = R * C;memset(marked, false, R * C * sizeof(bool));Position king;char r;fscanf(fin, "%c", &r);king.r = r - 'A';int c;fscanf(fin, "%d\n", &c);king.c = c - 1;assert(king.r < R && king.c < C);//fprintf(stdout, "king: \n%d %d\n", king.r, king.c);vector<Position> knights;while(fscanf(fin, "%c", &r) > 0){if(r < 'A' || r > 'Z'){continue;}Position knight;knight.r = r - 'A';fscanf(fin, " %d", &c);knight.c = c - 1;knights.push_back(knight);}//for(int i = 0; i < knights.size(); ++i)//{//assert(knights[i].r < R && knights[i].c < C);//fprintf(stdout, "%c %d\n", knights[i].r + 'A', knights[i].c + 1);//}//fprintf(stdout, "number of knights: %d\n", knights.size());for(int i = 0; i < knights.size(); ++i){spfa(knights[i], R, C);/*for(int r = 0; r < R; ++r){for(int c = 0; c < C; ++c){fprintf(stdout, "%d\t", dis[knights[i].r][knights[i].c][r][c]);}fprintf(stdout, "\n");}*/}int minmoves = -1;for(int gr = 0; gr < R; ++gr)//gathering point{for(int gc = 0; gc < C; ++gc){int knightsdis = 0;for(int i = 0; i < knights.size(); ++i){knightsdis += dis[knights[i].r][knights[i].c][gr][gc];}int kingmove = 0;int pr, pc;int knightsmove = 0;for(int kingmover = -2; kingmover <= 2; ++kingmover){for(int kingmovec = -2; kingmovec <= 2; ++kingmovec){kingmove = max(abs(kingmover), abs(kingmovec));//if(!kingmover && !kingmovec)//{//kingmove = 0;//}//else//{//kingmove = 1;//}pr = king.r + kingmover;pc = king.c + kingmovec;if(!isWithin(pr, pc, R, C)){continue;}spfa(Position(pr, pc), R, C);for(int pk = 0; pk < knights.size(); ++pk) //pickup knight{if(dis[knights[pk].r][knights[pk].c][pr][pc] == INT_MAX || dis[pr][pc][gr][gc] == INT_MAX){continue;}knightsmove = knightsdis - dis[knights[pk].r][knights[pk].c][gr][gc] + dis[knights[pk].r][knights[pk].c][pr][pc] + dis[pr][pc][gr][gc];int moves = knightsmove + kingmove;if(minmoves < 0 || moves < minmoves){minmoves = moves;}}}}}}fprintf(fout, "%d\n", minmoves < 0 ? 0 : minmoves);return 0;
}


这篇关于USACO Camelot 解题报告的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

【专题】2024飞行汽车技术全景报告合集PDF分享(附原数据表)

原文链接: https://tecdat.cn/?p=37628 6月16日,小鹏汇天旅航者X2在北京大兴国际机场临空经济区完成首飞,这也是小鹏汇天的产品在京津冀地区进行的首次飞行。小鹏汇天方面还表示,公司准备量产,并计划今年四季度开启预售小鹏汇天分体式飞行汽车,探索分体式飞行汽车城际通勤。阅读原文,获取专题报告合集全文,解锁文末271份飞行汽车相关行业研究报告。 据悉,业内人士对飞行汽车行业

usaco 1.3 Prime Cryptarithm(简单哈希表暴搜剪枝)

思路: 1. 用一个 hash[ ] 数组存放输入的数字,令 hash[ tmp ]=1 。 2. 一个自定义函数 check( ) ,检查各位是否为输入的数字。 3. 暴搜。第一行数从 100到999,第二行数从 10到99。 4. 剪枝。 代码: /*ID: who jayLANG: C++TASK: crypt1*/#include<stdio.h>bool h

usaco 1.3 Calf Flac(暴搜)

思路是暴搜。 需要注意的地方是输入的方法,以及输出时的换行。 代码: /*ID: who jayLANG: C++TASK: calfflac*/#include<stdio.h>#include<string.h>#include<math.h>int main(){freopen("calfflac.in","r",stdin);freopen("calfflac.ou

usaco 1.3 Barn Repair(贪心)

思路:用上M块木板时有 M-1 个间隙。目标是让总间隙最大。将相邻两个有牛的牛棚之间间隔的牛棚数排序,选取最大的M-1个作为间隙,其余地方用木板盖住。 做法: 1.若,板(M) 的数目大于或等于 牛棚中有牛的数目(C),则 目测 给每个牛牛发一个板就为最小的需求~ 2.否则,先对 牛牛们的门牌号排序,然后 用一个数组 blank[ ] 记录两门牌号之间的距离,然后 用数组 an

usaco 1.3 Mixing Milk (结构体排序 qsort) and hdu 2020(sort)

到了这题学会了结构体排序 于是回去修改了 1.2 milking cows 的算法~ 结构体排序核心: 1.结构体定义 struct Milk{int price;int milks;}milk[5000]; 2.自定义的比较函数,若返回值为正,qsort 函数判定a>b ;为负,a<b;为0,a==b; int milkcmp(const void *va,c

usaco 1.2 Palindromic Squares(进制转化)

考察进制转化 注意一些细节就可以了 直接上代码: /*ID: who jayLANG: C++TASK: palsquare*/#include<stdio.h>int x[20],xlen,y[20],ylen,B;void change(int n){int m;m=n;xlen=0;while(m){x[++xlen]=m%B;m/=B;}m=n*n;ylen=0;whi

usaco 1.2 Name That Number(数字字母转化)

巧妙的利用code[b[0]-'A'] 将字符ABC...Z转换为数字 需要注意的是重新开一个数组 c [ ] 存储字符串 应人为的在末尾附上 ‘ \ 0 ’ 详见代码: /*ID: who jayLANG: C++TASK: namenum*/#include<stdio.h>#include<string.h>int main(){FILE *fin = fopen (

usaco 1.2 Milking Cows(类hash表)

第一种思路被卡了时间 到第二种思路的时候就觉得第一种思路太坑爹了 代码又长又臭还超时!! 第一种思路:我不知道为什么最后一组数据会被卡 超时超了0.2s左右 大概想法是 快排加一个遍历 先将开始时间按升序排好 然后开始遍历比较 1 若 下一个开始beg[i] 小于 tem_end 则说明本组数据与上组数据是在连续的一个区间 取max( ed[i],tem_end ) 2 反之 这个

usaco 1.2 Transformations(模拟)

我的做法就是一个一个情况枚举出来 注意计算公式: ( 变换后的矩阵记为C) 顺时针旋转90°:C[i] [j]=A[n-j-1] [i] (旋转180°和270° 可以多转几个九十度来推) 对称:C[i] [n-j-1]=A[i] [j] 代码有点长 。。。 /*ID: who jayLANG: C++TASK: transform*/#include<

usaco 1.1 Broken Necklace(DP)

直接上代码 接触的第一道dp ps.大概的思路就是 先从左往右用一个数组在每个点记下蓝或黑的个数 再从右到左算一遍 最后取出最大的即可 核心语句在于: 如果 str[i] = 'r'  ,   rl[i]=rl[i-1]+1, bl[i]=0 如果 str[i] = 'b' ,  bl[i]=bl[i-1]+1, rl[i]=0 如果 str[i] = 'w',  bl[i]=b