yukicoder : No.1905 PURE PHRASE

ずいぶん変わった問題。
https://yukicoder.me/problems/no/1905

問題

ドレミファソラシド(C4～B4)のいずれかの音をシンセサイザーで鳴らした波形が、44100Hz・16bit・モノラルで1秒分与えられる。
どの音を鳴らしたか答えよ。

解法

変位の配列をAとする。
周波数をfとすると、T=44100/fサンプル毎に同じ値をとるはずなので、A[x]=A[x+T]となる。
実際はTは小数なので、
g(x+T) := A[x+floor(T)]とA[x+ceil(T)]を線形補間した値
を考える。
fを7通り総当たりし、abs(A[x]-g(x+T))の総和が最小なfを選ぼう。

int N;
double A[55555];
double S[]={261.6, 294.3, 327.0, 348.8, 392.4, 436.0, 490.5};
string pat[]={"C4", "D4", "E4", "F4", "G4", "A4", "B4"};

void solve() {
	int i,j,k,l,r,x,y; string s;
	
	cin>>N;
	FOR(i,N) cin>>A[i];

	double dif[7]={};
	FOR(i,7) {
		double span=N/S[i];
		FOR(j,N-1000) {
			double a=A[j+(int)floor(span)];
			double b=A[j+(int)ceil(span)];
			double c;
			double p=a*(1-modf(span,&c))+b*modf(span,&c);
			dif[i]+=abs(p-A[j]);
		}
	}
	
	cout<<pat[min_element(dif,dif+7)-dif]<<endl;
	
}
||

*まとめ