编程之美挑战赛初赛第一场题目1

时间限制:2000ms

单点时限:1000ms

内存限制:256MB

描述

一般来说，我们采用针孔相机模型，也就是认为它用到的是小孔成像原理。

在相机坐标系下，一般来说，我们用到的单位长度，不是“米”这样的国际单位，而是相邻像素的长度。而焦距在相机坐标系中的大小，是在图像处理领域的一个非常重要的物理量。

假设我们已经根据相机参数，得到镜头的物理焦距大小(focal length)，和相机胶片的宽度(CCD width)，以及照片的横向分辨率(image width)，则具体计算公式为：

Focal length in pixels = (image width in pixels) * (focal length on earth) / (CCD width on earth)

比如说对于Canon PowerShot S100, 带入公式得

Focal length in pixels = 1600 pixels * 5.4mm / 5.27mm = 1639.49 pixels

现在，请您写一段通用的程序，来求解焦距在相机坐标系中的大小。

输入

多组测试数据。首先是一个正整数T，表示测试数据的组数。

每组测试数据占一行，分别为

镜头的物理焦距大小(focal length on earth)

相机胶片的宽度(CCD width on earth)

照片的横向分辨率大小(image width in pixels)，单位为px。

之间用一个空格分隔。

输出

每组数据输出一行，格式为“Case X: Ypx”。 X为测试数据的编号，从1开始；Y为焦距在相机坐标系中的大小(focallength in pixels)，保留小数点后2位有效数字，四舍五入取整。

数据范围

对于小数据：focal length on earth和CCD width on earth单位都是毫米(mm)

对于大数据：长度单位还可能为米(m), 分米(dm), 厘米(cm), 毫米(mm), 微米(um)，纳米(nm)

样例输入

2
5.4mm 5.27mm 1600px
5400um 0.00527m 1600px

样例输出

Case 1: 1639.47px
Case 2: 1639.47px

代码如下：

//source her
#include<vector>
#include<map>
#include<iostream>
#include<algorithm>
#include<string>
using namespace std;
int main()
{
	int T;
	double b;
	double a=0;
	double danwei=0;
	string temp,dw1;
	int dw=0;
	cin>>T;
	vector<double> *p=new vector<double>[T];
	for(int i=0;i<T;++i)
	{
	  for(int j=0;j<3;++j)
	 {
		cin>>temp;
		//cout<<temp;
		dw1.clear();
		dw=0;
		a=0;
		b=1;
		for(string::iterator k=temp.end()-1;!(*k>='0'&&*k<='9');k--)
		{
			dw+=*k;
			dw1+=*k;
		}
		switch(dw)
		{
		case 'm'+'m':danwei=1;break;
		case 'm'+'u':danwei=0.001;break;
		case 'm':danwei=1000;break;
		case 'm'+'d':danwei=100;break;
		case 'm'+'c':danwei=10;break;
		case 'm'+'n':danwei=0.000001;break;
		case 'p'+'x':danwei=1;break;
		default :break;
		}

		string::iterator pos=find(temp.begin(),temp.end(),'.');
		if(pos==temp.end())
		{
		    b=1;
			for(string::iterator pos1=temp.end()-dw1.size();(--pos1)!=temp.begin();)
		   {
		     a+=(*pos1-'0')*b;
		     b*=10;
		   }
		   a+=(*temp.begin()-'0')*b;
		}
		else
		{
		   for(string::iterator pos1=pos;(--pos1)!=temp.begin();)
		  {
		  a+=(*pos1-'0')*b;
		  b*=10;
		  }
		a+=(*temp.begin()-'0')*b;
		b=0.1;
		  for(string::iterator pos2=pos;(++pos2)!=temp.end()-dw1.size();)
		  {
		   a+=((double)(*pos2-'0'))*b;
		   b*=0.1;
		  }
		}
	 a*=danwei;
	  p[i].push_back(a);
	 }
	  double result1=(p[i][2]*p[i][0])/p[i][1];
	  int result2=result1*100;
	  if(((int)(result1*1000)%10)>=5)
		++result2;
	  result1=(result2/100.0);
	  cout<<"Case "<<i+1<<": "<<result1<<"px"<<endl;
	}
	//system("pause"); 
	return 0;
}