QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#121431#5642. DA-Sortkaitlinz100 ✓0ms1612kbC5.2kb2023-07-08 06:10:272023-07-08 06:10:29

Judging History

你现在查看的是最新测评结果

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-07-08 06:10:29]
  • 评测
  • 测评结果:100
  • 用时:0ms
  • 内存:1612kb
  • [2023-07-08 06:10:27]
  • 提交

answer

/*
 * H - DA-Sort
 * Solution by Fred Pickel
 * ACM ICPC GNYR 2016
 */
#include <stdio.h>
#include <stdlib.h>
#include <ctype.h>
#include <string.h>
#include <math.h>

typedef unsigned long DWORD;

#define MAX_PROBS	1000
#define MAX_VALS	1000
#define VAL_LIMIT	1000000000

typedef struct _dasort_
{
	DWORD val;
	DWORD loc;
	struct _dasort_ *pNext;
} DASORT, *PDASORT;

typedef struct _rsort_
{
	PDASORT pHead;
	PDASORT pTail;
} RSORT;

DWORD prevmin, minval;
int minloc, aftermincnt;

#define RSHIFT	10		// should be chosen so 2^(3*RSHIFT) >= VAL_LIMIT
#define RSIZE	1024		//1 << RSHIFT
#define RMASK	0x3ff	// RSIZE - 1

DASORT values[MAX_VALS+1];
RSORT sorts[2][RSIZE];

void ClearRSort(int index) 
{
	int i;
	for(i = 0; i < RSIZE ; i++) {
		sorts[index][i].pHead = sorts[index][i].pTail = NULL;
	}
}

void RSortInsert(PDASORT pdas, DWORD val, int index)
{
	val &= RMASK;
	if(sorts[index][val].pHead == NULL) {
		sorts[index][val].pHead = sorts[index][val].pTail = pdas;
		pdas->pNext = NULL;
	} else {
		sorts[index][val].pTail->pNext = pdas;
		pdas->pNext = NULL;
		sorts[index][val].pTail = pdas;
	}
}

char inbuf[256];

int main()
{
	int nprob, curprob, index, nvals, count;
	int rem, cursz, i, j, curstart;
	DWORD vals[10], ind, prevloc, firstMoved;
	PDASORT pdas, nxtpdas;
	if(fgets(&(inbuf[0]), 255, stdin) == NULL)
	{
		fprintf(stderr, "Read failed on problem count\n");
		return -1;
	}
	if(sscanf(&(inbuf[0]), "%d", &nprob) != 1)
	{
		fprintf(stderr, "Scan failed on problem count\n");
		return -2;
	}
	if(nprob > MAX_PROBS) {
		fprintf(stderr, "Problem count %d > Max allowed %d\n", nprob, MAX_PROBS);
		return -3;
	}
	for(curprob = 1; curprob <= nprob ; curprob++)
	{
		if(fgets(&(inbuf[0]), 255, stdin) == NULL)
		{
			fprintf(stderr, "Read failed on problem %d header\n", curprob);
			return -4;
		}
		if(sscanf(&(inbuf[0]), "%d %d", &index, &nvals) != 2) {
			fprintf(stderr, "scan failed on problem %d header\n", curprob);
			return -5;
		}
		if(index != curprob) {
			fprintf(stderr, "problem index %d not = expected problem %d\n",
				index, curprob);
			return -6;
		}
		if(nvals > MAX_VALS) {
			fprintf(stderr, "Problem %d val count %d > Max allowed %d\n", curprob, nvals, MAX_VALS);
			return -7;
		}
		// read problem data
		rem = nvals; i= 0;
		curstart = 1;
		prevmin = minval = VAL_LIMIT + 1;
		minloc = -1;
		aftermincnt = 0;
		while(rem > 0) {
			i++;
			cursz = rem;
			if(cursz > 10){
				cursz = 10;
			} 
			rem -= cursz;
			if(fgets(&(inbuf[0]), 255, stdin) == NULL)
			{
				fprintf(stderr, "Read failed on problem %d data line %d\n", curprob, i);
				return -8;
			}
			if(sscanf(&(inbuf[0]), "%lu %lu %lu %lu %lu %lu %lu %lu %lu %lu ",
				&(vals[0]), &(vals[1]), &(vals[2]), &(vals[3]), &(vals[4]), 
				&(vals[5]), &(vals[6]), &(vals[7]), &(vals[8]), &(vals[9])) != cursz) {
				fprintf(stderr, "scan failed on problem %d data line %d\n", curprob, i);
				return -9;
			}
			for(j = 0; j < cursz; j++) {
				if(vals[j] > VAL_LIMIT) {
					fprintf(stderr, "problem %d value %d = %lu too large\n", curprob, 10*(i-1)+j, vals[j]);
					return -10;
				}
				values[curstart + j].val = vals[j];
				values[curstart + j].loc = curstart+j;
				if(vals[j] < minval) {
					prevmin = minval;
					minval = vals[j];
					minloc = curstart+j;
					aftermincnt = 0;
				} else if(vals[j] <= prevmin) {
					aftermincnt++;
				}
			}
			curstart += cursz;
		}
		// radix sort value after minloc that are less than prevmin
		// prevmin has to move to end so anything larger has to move as well
		ClearRSort(0);
		pdas = &(values[minloc]);
		for(i = minloc; i <= nvals ; i++, pdas++) {
			if(pdas->val <= prevmin) {
				ind = pdas->val & RMASK;
				RSortInsert(pdas, ind, 0);
			}
		}
		ClearRSort(1);
		for(i = 0; i < RSIZE; i++) {
			pdas = sorts[0][i].pHead;
			while(pdas != NULL) {
				nxtpdas = pdas->pNext;
				ind = (pdas->val >> RSHIFT) & RMASK;
				RSortInsert(pdas, ind, 1);
				pdas = nxtpdas;
			}
		}
		ClearRSort(0);
		for(i = 0; i < RSIZE; i++) {
			pdas = sorts[1][i].pHead;
			while(pdas != NULL) {
				nxtpdas = pdas->pNext;
				ind = (pdas->val >> (2*RSHIFT)) & RMASK;
				RSortInsert(pdas, ind, 0);
				pdas = nxtpdas;
			}
		}
		// candidates to not move are now sorted
		// scan through until the original position of some element is less than original
		// position of its predecessor When that happens we remeber the value (as firstMoved)
		// any entries with the same value that came after the smaller value do not need to be moved
		// any larger values will need to be moved
		count = prevloc = 0;
		firstMoved = VAL_LIMIT + 1;
		for(i = 0; (i < RSIZE) && (prevloc < MAX_VALS) ; i++) {
			pdas = sorts[0][i].pHead;
			while(pdas != NULL) {
				if((pdas->loc >= prevloc) && (pdas->val <= firstMoved)) {
					count++;
					prevloc = pdas->loc;
				} else if(pdas->val > firstMoved){
					prevloc = MAX_VALS;
					break;
				} else {	// not tooo big but out of order so needs to be moved and all larger need to move
					firstMoved = pdas->val;
				}
				pdas = pdas->pNext;
			}
			pdas = nxtpdas;
		}
		printf("%d %d\n", curprob, nvals - count);
	}
	return 0;
}

详细

Test #1:

score: 100
Accepted
time: 0ms
memory: 1612kb

input:

50
1 3
1 3 2
2 6
1 5 2 4 3 6
3 23
67890 56312 999999999 12345 23456 38927 45632 100345 98765 23456
87654 43278 23456 117654 321899 25432 54326 217435 26845 31782
33456 41234 56213
4 1000
18780 15256 15263 15268 15269 15274 15279 15279 15283 15290
15292 15292 15295 15298 15302 15305 15305 15311 15316...

output:

1 1
2 3
3 15
4 1
5 0
6 0
7 0
8 0
9 999
10 1
11 500
12 2
13 299
14 159
15 36
16 564
17 94
18 0
19 140
20 438
21 182
22 364
23 364
24 45
25 147
26 16
27 66
28 85
29 406
30 325
31 564
32 268
33 9
34 248
35 524
36 335
37 498
38 43
39 392
40 589
41 157
42 148
43 41
44 576
45 177
46 332
47 187
48 22
49 99...

result:

ok 50 lines