USACO section 1.4.3 Arithmetic Progressions
来源:互联网 发布:郑州财经学院网络管 编辑:程序博客网 时间:2024/06/06 05:11
1. 枚举可以算出来,但是肯定超时,所以要剪枝,以等差数列a+nb说明;
2. 几个可以剪枝的地方:b的范围是从1~2*m*m/(n-1);
a的范围是0~2*m*m - (n-1) * b
然后循环判断是不是pq了
剪枝不是很充分,但是最慢3秒多给过了,哈哈
3.以下是代码:
/*ID: dollar4PROG: ariprogLANG: C++*/#include <iostream>#include <fstream>#include <string>#include <algorithm>#include <cstring>using namespace std;int main(){ ofstream fout ("ariprog.out"); ifstream fin ("ariprog.in"); int is[125010], sum = 0, i, j, k, m, n; memset(is, 0, sizeof(is)); fin >> n >> m; for (i = 0; i <= m; i++) for (j = 0; j <= m; j++) { is[i*i+j*j] = 1; } for (i = 1; i <= (2*m*m)/(n-1); i++)//等差数列a+nb中的b { int e = (2 * m * m) - (n - 1) * i; for (j = 0; j <= e; j++) //a { int flag = 0; for (k = j; k <= j + (n - 1) * i; k += i) { if (!is[k]) { flag = 1; break; } } if (flag == 0) { fout << j << ' ' << i << endl; sum++; } } } if (sum == 0) fout << "NONE" << endl; return 0;}
4. 官方参考代码:
This can be done by brute force, enumerating over all possible sequences. It has to be done a little carefully in order to get it to run in time, however.
Precalculate the bisquares, first of all. Calculate both a sorted list of the bisquares, along with a boolean array saying whether each number between 1 and 125000 (the maximum bisquare possible, for p,q < 250) is a bisquare.
Go through the skips in increasing order, starting at 1, and continuing along as the sequence starting at 1 with that skip doesn't exceed the maximum bisquare. For each bisquare, determine if the sequence starting at that location and with the current skip value consists of all of bisquares. If so, output it.
Here is the solution of Felix Arends from Germany (modified by Iran's Saber Fadaee):
#include <stdio.h>#include <assert.h>#include <string>using namespace std;// open filesFILE *fin = fopen ("ariprog.in", "r");FILE *fout = fopen ("ariprog.out", "w");// global variablesunsigned int N, M, maxMM;unsigned int numbers [65000];unsigned int number_size = 0;unsigned char num_available [125001];unsigned char dist_available [125001];int have_res = 0;int skipstep = 1;// read the inputint read_input () { fscanf (fin, "%d %d", &N, &M); return 0;}int cmp_int (const void *a, const void *b) { return (*(int *)a - *(int *)b);}void asm_num (int a, int b) { for (unsigned int n = 1; n < N; n++) if (num_available [a + n * b] == 0) return; fprintf (fout, "%d %d\n", a, b); have_res ++; if (have_res==1) skipstep = b;}void asm_num () { for (unsigned int b = 1; b < maxMM; b+=skipstep) { if (dist_available [b]) { for (unsigned int p = 0; p < number_size && numbers [p] + (N -1) * b <= maxMM; p++) asm_num (numbers [p], b); } }}int process () { memset (num_available, 0, sizeof (unsigned char) * 125001); memset (dist_available, 0, sizeof (unsigned char) * 125001); for (unsigned int m1 = 0; m1 <= M; m1++) { for (unsigned int m2 = m1; m2 <= M; m2++) { int n = m1 * m1 + m2 * m2; if (!num_available [n]) { num_available [n] = 1; numbers [number_size++] = n; } } } qsort (numbers, number_size, sizeof (unsigned int), cmp_int); maxMM = M * M + M * M; for (unsigned int n1 = 0; n1 < number_size; n1++) { unsigned int _n1 = numbers [n1]; for (unsigned int n2 = n1 + 1; n2 < number_size && _n1 + (numbers[n2] - _n1) * (N - 1) <= maxMM; n2++) { assert (numbers [n2] - _n1 >= 0 && numbers [n2] - _n1 < 125001); if (num_available [_n1 + (numbers [n2] - _n1) * (N - 1)] && num_available [_n1 + (numbers [n2] - _n1) * (N - 2)]) dist_available [numbers [n2] - _n1] = true; } } asm_num (); if (!have_res) fprintf (fout, "NONE\n"); return 0;}int main () { read_input (); process (); fclose (fin); fclose (fout); return 0;}
Here is an even faster solution of "UaE ProGrammeR":
#include <fstream>#include <iostream>using namespace std;void quicksort (int[], int ,int);int pivotlist (int[], int,int);ofstream out ("ariprog.out"); int n;int main () { ifstream in ("ariprog.in"); bool array[125001] = {false}, noneF; int m, upper, upperdef, def, p; int places[300000], pl = 0; noneF = true; in>>n>>m; for (int i = 0; i <= m; i++) for (int j = 0; j <= m; j++) { if (!array[i*i+j*j]) { places[pl] = i*i+j*j; //Saving generated numbers pl++; } array[i*i+j*j] = true; } upper = 2*m*m; upperdef = upper/(n-1); quicksort (places, 0, pl-1); for ( def = 1; def<=upperdef; def++) // Loop to check for solutions // It looks for solutions in // correct order so you // print the solution directly // without sorting first, thnx to who said: // Trade Memory for Speed !! { for ( p = 0; places[p]<=(upper-((n-2)*def)); p++) { bool is; is = true; int where; for (int c = (n-1); c>=0 ; c--) if (!array[places[p]+c*def]) { is = false; where = (p+c*def); break; } if (is) { noneF = false; out<<places[p]<<" "<<def<<endl; } } } if (noneF) out<<"NONE"<<endl; return 0;}void quicksort (int array[], int start, int last) { int pivot; if (start < last) { pivot = pivotlist(array, start,last); quicksort (array, start,pivot-1); quicksort (array, pivot+1,last); }}int pivotlist(int array[], int f, int l) { int pivotpoint; int pivotvalue, temp; pivotvalue = array[f]; pivotpoint = f; for (int i = f+1;i<=l; i++) { if (array[i]<pivotvalue) { pivotpoint++; temp = array[i]; array[i] = array[pivotpoint]; array[pivotpoint] = temp; } } temp = array[f]; array[f] = array[pivotpoint]; array[pivotpoint] = temp; return pivotpoint;}
- USACO section 1.4.3 Arithmetic Progressions
- USACO Section 1.4 Arithmetic Progressions
- USACO-Section 1.4 Arithmetic Progressions
- USACO Section 1.4 Arithmetic Progressions
- USACO Section 1.4 Arithmetic Progressions (ariprog)
- USACO Section 1.4 - Arithmetic Progressions (暴力)
- USACO-Section 1.4 Arithmetic Progressions(枚举)
- [USACO 1.4.3] Arithmetic Progressions
- Section 1.4 Arithmetic Progressions
- USACO 1.4-Arithmetic Progressions
- [usaco]1.4 Arithmetic Progressions
- usaco 1.4 Arithmetic Progressions
- USACO 1.4 Arithmetic Progressions
- USACO 1.4 Arithmetic Progressions
- USACO 1.4-Arithmetic Progressions
- USACO 1.4 Arithmetic Progressions
- Usaco 1.4.3 等差数列(Arithmetic Progressions)
- USACO:1.4.3 Arithmetic Progressions 等差数列 解析
- 图片二进制互相转换C#
- Java中final、finally 和 finalize
- Windows Phone 7的XML操作详解:读取,过滤以及数据绑定
- 怎样判断ios app 第一次启动
- MPICH安装后jumpshot未安装的解决方法
- USACO section 1.4.3 Arithmetic Progressions
- iPhone程序包结构
- StringBuffer、StringBuilder与String剖析
- TCO2012Round1B-2-FoxAndDoraemon
- 无刷新对表格的一列排序
- 使用CSS3让文字按指定路径显示
- PHP Mysql Apache 配置
- Qt中Ui名字空间以及setupUi函数的原理和实现
- HDU-1279 模拟