三峡大学校OJ 79.数字列表 (15分)

79.数字列表 (15分)

C时间限制：3000 毫秒 | C内存限制：3000 Kb

题目内容：
给出包含M个数字的列表，和列表中所有数字的所有质因数。求出最长的子列表，使得子列表中所有数字的乘积是一个完全平方数.

输入描述
输入文件包含多组测试数据。第一行包含两个整数N , M ( 1 <= N <= 30 , 1 <= M <= 30000 ). N
是质因数的个数。接下来一行有N个整数，给出所有的质因数。然后一行包含M个整数，给出列表。

输出描述
对于每组数据，输出最长子列表的两个位置坐标l r。l是该子列表在列表中的起始位置，r是结束位置
。如果多种情况都满足子列表长度最大，输出l最小的一个。如果不存在这样的子列表输出“None”。

输入样例
3 4
2 3 5
4 9 25 6
3 4
2 3 5
6 6 3 3
0 0

输出样例
1 3
1 4

思路：这道题当然不能做乘法来判断完全平方数，最多30000个数相乘，肯定会爆的，完全平方数其实可以分成许多质因数相乘，每个质因数的指数一定是偶数，我们把列表每个数所有质因数的指数得到，放在结构体数组中，可以知道一个合乎条件的区间，里面所有元素所分解出的单个质因数的指数的和一定是偶数，这就是判断条件。而关于如何找到最大的左右区间，我实在没想到特别巧妙的方法，只能先从最大的区间找起，一点点缩小规模，从下标小的找到大的，暴力枚举出来，可以说用到了贪心吧（笑）。

#include<iostream>
using namespace std; 
int n,m,prm[35];
typedef struct ll
{
	int prm[35] = { 0 };
	int zhi;
}list;
list a[30005];
int main()
{
	while (cin >> n >> m)
	{
		if (m == 0 && n == 0)
		{
			return 0;
		}
		
		for (int i = 0; i < n; i++)
		{
			cin>>prm[i];
		}
		for (int i = 0; i < m; i++)
		{
			cin >> a[i].zhi;
			for (int j = 0; j < n; j++)
			{
				a[i].prm[j]=0;
				while(a[i].zhi%prm[j]==0&&a[i].zhi!=1)
				{
					a[i].prm[j]++;
					a[i].zhi/=prm[j];
				}
			}
		}
		int k = 0,pd=0;
		for (int i = m; i >= 1; i--)
		{
			for (int j = 0; j <= m - i; j++)
			{
				for (int p = 0; p < n; p++)
				{
					k = 0;
					for (int v = j; v < i + j; v++)
					{
						k += a[v].prm[p];
					}
					if (k % 2 != 0)break;
					if (p == n - 1)
					{
						cout << j+1 << " " << j + i << endl;
						pd = 1;
						break;
					}
				}
				if (pd == 1)break;
			}
			if (pd == 1)break;
		}
		if(pd==0)cout<<"None"<<endl; 
	}
}