algorithm is this integration - andstudy/forge GitHub Wiki

[http://www.programming-challenges.com/pg.php?page=downloadproblem&probid=111307&format=html 문제보기] [http://www.programming-challenges.com/pg.php?page=studenthome&threshold=0&perpage=8&pagenum=12&orderby=12 13장 Summit Pages]

Mastojun

문제 풀이

백투더 대학교 1학년!
오랫만에 적분표꺼내들고 풀었습니다
Programming-Challenges에서는 WA, UVa에서는 AC뜨네요.
소스는 민망합니다.
대부분의 풀이는 종이와 계산기로.

소스 코드

#include <stdio.h>

    int main()
    {
    	double in;
    
    	double s, d, r;
    
    	while( scanf("%lf", &in) >= 1 )
    	{
    		in *= in;
    		s = in*0.31514674363;
    		d = in*0.51129916633;
    		r = in*0.17355409004;
    
    		printf("%.3lf %.3lf %.3lf\n", s, d, r);
    	}
    
    	return 0;
    }

CynicJJ

wrong answer. 왜?
- UVa에서는 Accepted뜨네요! :D PC쪽이 이상한듯...
  - 우왕ㅋ굳ㅋ 고맙습니다 ^^
여기저기 선 그어서 풀이

적분인가? 라는 질문에 대한 대답은 아니오

  #include <iostream>
  #include <math.h>
  
  using namespace std;
  
  const double PI = 3.14159265359;
  
  double getStripedHeight(double radius)
  {
  	double height = radius * (sqrt(3.0)-1) / 2;
  	return height;
  }
  
  double getStriped(double radius)
  {
  	double circlePart = (PI*radius*radius) / 12;
  	double triangle = radius/2 * getStripedHeight(radius) / 2;
  	double striped = 4 * (circlePart - 2*triangle);
  	return striped;
  }
  
  double getDotted(double radius)
  {
  	double circlePart = (PI*radius*radius)/4 - (radius*radius)/2;
  	double dotted = 4 * (circlePart - getStriped(radius)/2);
  	return dotted;
  }
  
  int main()
  {
  	double length;
  	while (cin >> length)
  	{
  		double striped = getStriped(length);
  		double dotted = getDotted(length);
  		double rest = length*length - striped - dotted;
  		printf("%.3f %.3f %.3f\n", striped, dotted, rest);
  	}
  	return 0;
  }

Outbreak

문제 요약

원과 정사각형이 교차 할 때, 각 부분 영역의 넓이 구하기

문제 해결

삼각형과 부채꼴의 볼록부분 넓이를 구해서 해결

소감

답에 차이가 나는 이유를 모르겠음;
- 아 잘못 풀었네요.. 올리다 문제점 발견 흐..
- 음, 제대로 고친거 같은데.. Wrong Answer.. 흐..;
  - UVa결과 Accepted! Programming-challenges 홈페이지에 문제가 있는듯^^;

풀이

// "Is This Integration" Solution By Outbreak 

#include <iostream>
#include <iomanip>
#include <cmath>

using namespace std;

//#define _TDD

//---------------------------------------------------------------------------

namespace Solver
{
	static const double PI = 3.14159265358979;
	
	class ABCD
	{
	public:
		ABCD(double r) : r(r) {}

	public:
		// 각 영역 넓이
		double GetStripedArea()
		{
			return GetRectangleArea() + GetChordArea() * 4;
		}

		double GetDottedArea()
		{
			return  GetTriangleArea() * 4 + GetChordArea() * 4;
		}

		double GetRestArea()
		{
			return r * r - GetStripedArea() - GetDottedArea();
		}

	private:
		// 부채꼴 ABD, CDB
		double GetBigFanArea()
		{
			return (r * r * PI) / 4;
		}

		// 삼각형 ABD, CBD
		double GetBigTriangleArea()
		{
			return r * r / 2;
		}

		// striped 영역 내부 정사각형 넓이
		double GetRectangleArea()
		{			
			return GetChordLength() * GetChordLength();
		}

		// dotted 영역 내부 정삼각형 넓이
		double GetTriangleArea()
		{
			return sqrt((double)3)*GetChordLength()*GetChordLength()/4;
		}
	
		// striped 사각형 또는 dotted 삼각형의 볼록 부분 넓이
		double GetChordArea()
		{		
			return (GetBigFanArea() - GetBigTriangleArea() - GetRectangleArea()/2 - GetTriangleArea() ) / 3;
		}

		double GetChordLength()
		{
			return  (r/2) / cos(PI/12);
		}

	private:
		double r;
	};

	void Calculate(const double r, double& striped, double& dotted, double& rest)
	{
		ABCD abcd(r);

		// striped		
		striped = abcd.GetStripedArea();

		// dotted
		dotted = abcd.GetDottedArea();

		// rest
		rest = abcd.GetRestArea();
	}
}

//---------------------------------------------------------------------------

namespace Runner
{
	void Execute(istream& in, ostream& out)
	{
		double a;

		while( in >> a )
		{
			if( in.fail() )
				break;

			double striped = 0, dotted = 0, rest = 0;

			// 답 구하기
			Solver::Calculate(a, striped, dotted, rest);

			out << setprecision(3) << fixed << striped << " " << dotted << " " << rest << endl;	
		}
	}
}

//---------------------------------------------------------------------------

#ifdef _TDD

#include "unittest++.h"

TEST(Output1)
{
	stringstream input;
    stringstream output;

    input << "0.1";

	Runner::Execute(input, output);

	CHECK_EQUAL("0.003 0.005 0.002\n", output.str());           
}

TEST(Output)
{
	stringstream input;
    stringstream output;

    input << "0.1\n0.2\n0.3\n";

	Runner::Execute(input, output);

	CHECK_EQUAL("0.003 0.005 0.002\n0.013 0.020 0.007\n0.028 0.046 0.016\n", output.str());           
}

#endif

//---------------------------------------------------------------------------

int main()
{

#ifdef _TDD	
	UnitTest::RunAllTests();
#else
	Runner::Execute(cin, cout);
#endif // _TDD
	return 0;

}

itmentor

참고 사이트

적분계산기

소스(C++)

    #include <iostream>
    #include <cmath>
    
    using namespace std;
    
    float centerTriangleArea(float area)
    {
    	//sqrt( 1-x^2 ) --> (x*sqrt(1 - x^2) + arcsin(x))/2;
    	//x=[0.5.. sqrt(3) / 2]
    	float center = 0.26179938090693;
    	float baseArea = 0.5 * ( sqrt(3.0) / 2.0  - 0.5 );
    
    	return 4*area*(center - baseArea);
    }
    
    float topTriangleArea(float area)
    {
    	//sqrt( 1-x^2 ) --> (x*sqrt(1 - x^2) + arcsin(x))/2;
    	//x=[0..0.5]
    	float top = 0.47830573874525906;
    	return area * 8*(0.5 - top);	//8조각
    }
    
    float midTriangleArea(float area)
    {
    	return area  - centerTriangleArea(area) - topTriangleArea(area);
    }
    
    int main()
    {
    	float scale;
    	while( std::cin >> scale )
        {
                if( std::cin.fail() )
                        break;
    
    			scale *= scale;	//area
    
     			float center = centerTriangleArea(scale);
    			float top = topTriangleArea(scale);
    			float mid = midTriangleArea(scale);
    
    			printf("%.3f %.3f %.3f\n", center, mid, top );
        }
    
    
    	return 0;
    }