提出 #44348463

---

ソースコード 拡げる

#include <stdio.h>
 
void x87_sincos(double* sin_res, double* cos_res, double angle) {
	__asm__ __volatile__(
		"fldl %2\n\t"
		"fsincos\n\t"
		"fstpl %1\n\t"
		"fstpl %0\n\t"
	: "=m"(*sin_res), "=m"(*cos_res) : "m"(angle));
}
 
double x87_sqrt(double value) {
	double res;
	__asm__ __volatile__(
		"fldl %1\n\t"
		"fsqrt\n\t"
		"fstpl %0\n\t"
	: "=m"(res) : "m"(value));
	return res;
}
 
double x87_atan2(double y, double x) {
	double res;
	__asm__ __volatile__(
		"fldl %1\n\t"
		"fldl %2\n\t"
		"fpatan\n\t"
		"fstpl %0\n\t"
	: "=m"(res) : "m"(y), "m"(x));
	return res;
}
 
int main(void) {
	double PI = x87_atan2(1, 1) * 4;
	double T;
	double L, X, Y;
	int Q;
	int i;
	if (scanf("%lf", &T) != 1) return 1;
	if (scanf("%lf%lf%lf", &L, &X, &Y) != 3) return 1;
	if (scanf("%d", &Q) != 1) return 1;
	for (i = 0; i < Q; i++) {
		double E;
		double s, c;
		double y, z;
		double teihen, takasa;
		if(scanf("%lf", &E) != 1) return 1;
		x87_sincos(&s, &c, 2.0 * PI * E / T);
		y = -(L / 2) * s;
		z = (L / 2) - (L / 2) * c;
		teihen = x87_sqrt(X * X + (Y - y) * (Y - y));
		takasa = z;
		printf("%.15f\n", x87_atan2(takasa, teihen) * 180.0 / PI);
	}
	return 0;
}

#include <stdio.h>

void x87_sincos(double* sin_res, double* cos_res, double angle) {
	__asm__ __volatile__(
		"fldl %2\n\t"
		"fsincos\n\t"
		"fstpl %1\n\t"
		"fstpl %0\n\t"
	: "=m"(*sin_res), "=m"(*cos_res) : "m"(angle));
}

double x87_sqrt(double value) {
	double res;
	__asm__ __volatile__(
		"fldl %1\n\t"
		"fsqrt\n\t"
		"fstpl %0\n\t"
	: "=m"(res) : "m"(value));
	return res;
}

double x87_atan2(double y, double x) {
	double res;
	__asm__ __volatile__(
		"fldl %1\n\t"
		"fldl %2\n\t"
		"fpatan\n\t"
		"fstpl %0\n\t"
	: "=m"(res) : "m"(y), "m"(x));
	return res;
}

int main(void) {
	double PI = x87_atan2(1, 1) * 4;
	double T;
	double L, X, Y;
	int Q;
	int i;
	if (scanf("%lf", &T) != 1) return 1;
	if (scanf("%lf%lf%lf", &L, &X, &Y) != 3) return 1;
	if (scanf("%d", &Q) != 1) return 1;
	for (i = 0; i < Q; i++) {
		double E;
		double s, c;
		double y, z;
		double teihen, takasa;
		if(scanf("%lf", &E) != 1) return 1;
		x87_sincos(&s, &c, 2.0 * PI * E / T);
		y = -(L / 2) * s;
		z = (L / 2) - (L / 2) * c;
		teihen = x87_sqrt(X * X + (Y - y) * (Y - y));
		takasa = z;
		printf("%.15f\n", x87_atan2(takasa, teihen) * 180.0 / PI);
	}
	return 0;
}

提出情報

ジャッジ結果