4fd800b3a8 2011-02-23        kinaba: #include <iostream>
4fd800b3a8 2011-02-23        kinaba: #include <sstream>
4fd800b3a8 2011-02-23        kinaba: #include <iomanip>
4fd800b3a8 2011-02-23        kinaba: #include <vector>
4fd800b3a8 2011-02-23        kinaba: #include <string>
4fd800b3a8 2011-02-23        kinaba: #include <map>
4fd800b3a8 2011-02-23        kinaba: #include <set>
4fd800b3a8 2011-02-23        kinaba: #include <algorithm>
4fd800b3a8 2011-02-23        kinaba: #include <numeric>
4fd800b3a8 2011-02-23        kinaba: #include <iterator>
4fd800b3a8 2011-02-23        kinaba: #include <complex>
4fd800b3a8 2011-02-23        kinaba: #include <queue>
4fd800b3a8 2011-02-23        kinaba: #include <stack>
4fd800b3a8 2011-02-23        kinaba: #include <cmath>
4fd800b3a8 2011-02-23        kinaba: using namespace std;
4fd800b3a8 2011-02-23        kinaba: typedef long long LL;
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: static const LL MODVAL = 1234567891;
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: LL ADD(LL x, LL y) { return (x+y)%MODVAL; }
4fd800b3a8 2011-02-23        kinaba: LL MUL(LL x, LL y) { return (x*y)%MODVAL; }
4fd800b3a8 2011-02-23        kinaba: LL FACT(LL x) { return x<=1 ? 1 : MUL(x, FACT(x-1)); }
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: class PSequence
4fd800b3a8 2011-02-23        kinaba: {
4fd800b3a8 2011-02-23        kinaba: public:
4fd800b3a8 2011-02-23        kinaba: 	map<pair<vector<LL>,int>, LL> memo;
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: 	int count(vector <int> S, int p)
4fd800b3a8 2011-02-23        kinaba: 	{
4fd800b3a8 2011-02-23        kinaba: 		memo.clear();
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: 		vector<LL> cs; // congruences mod p
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: 		bool used[30] = {};
4fd800b3a8 2011-02-23        kinaba: 		for(int i=0; i<S.size(); ++i)
4fd800b3a8 2011-02-23        kinaba: 			if( !used[i] )
4fd800b3a8 2011-02-23        kinaba: 			{
4fd800b3a8 2011-02-23        kinaba: 				int cnt = 0;
4fd800b3a8 2011-02-23        kinaba: 				for(int j=i; j<S.size(); ++j)
4fd800b3a8 2011-02-23        kinaba: 					if( abs(S[i]-S[j]) % p == 0 )
4fd800b3a8 2011-02-23        kinaba: 						{used[j]=true; ++cnt;}
4fd800b3a8 2011-02-23        kinaba: 				cs.push_back(cnt);
4fd800b3a8 2011-02-23        kinaba: 			}
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: 		LL f = num_perm( cs, -1, S.size() );
4fd800b3a8 2011-02-23        kinaba: 		for(int i=0; i<cs.size(); ++i)
4fd800b3a8 2011-02-23        kinaba: 			f = MUL(f, FACT(cs[i]));
4fd800b3a8 2011-02-23        kinaba: 		return f;
4fd800b3a8 2011-02-23        kinaba: 	}
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: 	LL num_perm( vector<LL>& cs, int forbidden, int n )
4fd800b3a8 2011-02-23        kinaba: 	{
4fd800b3a8 2011-02-23        kinaba: 		if( n == 0 )
4fd800b3a8 2011-02-23        kinaba: 			return 1;
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: 		int v = (forbidden >= 0 ? cs[forbidden] : -1);
4fd800b3a8 2011-02-23        kinaba: 		pair<vector<LL>,int> key(cs, forbidden);
4fd800b3a8 2011-02-23        kinaba: 		sort(key.first.begin(), key.first.end());
4fd800b3a8 2011-02-23        kinaba: 		key.second = find(key.first.begin(), key.first.end(), v) - key.first.begin();
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: 		if( memo.count(key) )
4fd800b3a8 2011-02-23        kinaba: 			return memo[key];
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: 		LL np = 0;
4fd800b3a8 2011-02-23        kinaba: 		for(int i=0; i<cs.size(); ++i)
4fd800b3a8 2011-02-23        kinaba: 			if( cs[i] && i!=forbidden ) {
4fd800b3a8 2011-02-23        kinaba: 				cs[i]--;
4fd800b3a8 2011-02-23        kinaba: 				np = ADD(np, num_perm( cs, i, n-1 ));
4fd800b3a8 2011-02-23        kinaba: 				cs[i]++;
4fd800b3a8 2011-02-23        kinaba: 			}
4fd800b3a8 2011-02-23        kinaba: 		return memo[key] = np;
4fd800b3a8 2011-02-23        kinaba: 	}
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: // BEGIN CUT HERE
4fd800b3a8 2011-02-23        kinaba: 	public:
4fd800b3a8 2011-02-23        kinaba: 	void run_test(int Case) { if ((Case == -1) || (Case == 0)) test_case_0(); if ((Case == -1) || (Case == 1)) test_case_1(); if ((Case == -1) || (Case == 2)) test_case_2(); if ((Case == -1) || (Case == 3)) test_case_3(); if ((Case == -1) || (Case == 4)) test_case_4(); }
4fd800b3a8 2011-02-23        kinaba: 	private:
4fd800b3a8 2011-02-23        kinaba: 	template <typename T> string print_array(const vector<T> &V) { ostringstream os; os << "{ "; for (typename vector<T>::const_iterator iter = V.begin(); iter != V.end(); ++iter) os << '\"' << *iter << "\","; os << " }"; return os.str(); }
4fd800b3a8 2011-02-23        kinaba: 	void verify_case(int Case, const int &Expected, const int &Received) { cerr << "Test Case #" << Case << "..."; if (Expected == Received) cerr << "PASSED" << endl; else { cerr << "FAILED" << endl; cerr << "\tExpected: \"" << Expected << '\"' << endl; cerr << "\tReceived: \"" << Received << '\"' << endl; } }
4fd800b3a8 2011-02-23        kinaba: 	void test_case_0() { int Arr0[] = {-1,0,1,2,3}; vector <int> Arg0(Arr0, Arr0 + (sizeof(Arr0) / sizeof(Arr0[0]))); int Arg1 = 10; int Arg2 = 120; verify_case(0, Arg2, count(Arg0, Arg1)); }
4fd800b3a8 2011-02-23        kinaba: 	void test_case_1() { int Arr0[] = {6,2}; vector <int> Arg0(Arr0, Arr0 + (sizeof(Arr0) / sizeof(Arr0[0]))); int Arg1 = 4; int Arg2 = 0; verify_case(1, Arg2, count(Arg0, Arg1)); }
4fd800b3a8 2011-02-23        kinaba: 	void test_case_2() { int Arr0[] = {1,2,3,4}; vector <int> Arg0(Arr0, Arr0 + (sizeof(Arr0) / sizeof(Arr0[0]))); int Arg1 = 3; int Arg2 = 12; verify_case(2, Arg2, count(Arg0, Arg1)); }
4fd800b3a8 2011-02-23        kinaba: 	void test_case_3() { int Arr0[] = {4,6,8,-3,7}; vector <int> Arg0(Arr0, Arr0 + (sizeof(Arr0) / sizeof(Arr0[0]))); int Arg1 = 2; int Arg2 = 12; verify_case(3, Arg2, count(Arg0, Arg1)); }
4fd800b3a8 2011-02-23        kinaba: 	void test_case_4() {
4fd800b3a8 2011-02-23        kinaba:  int Arr0[] = {0,1,2,3,4,5,6,7,8,9,
4fd800b3a8 2011-02-23        kinaba: 10,11,12,13,14,15,16,17,18,19,
4fd800b3a8 2011-02-23        kinaba: 20,21,22,23,24,25,26,27,28,29};
4fd800b3a8 2011-02-23        kinaba: vector <int> Arg0(Arr0, Arr0 + (sizeof(Arr0) / sizeof(Arr0[0])));
4fd800b3a8 2011-02-23        kinaba:  int Arg1 = 15; int Arg2 = 12; verify_case(4, Arg2, count(Arg0, Arg1)); }
4fd800b3a8 2011-02-23        kinaba: 
4fd800b3a8 2011-02-23        kinaba: // END CUT HERE
4fd800b3a8 2011-02-23        kinaba: };
4fd800b3a8 2011-02-23        kinaba: // BEGIN CUT HERE
4fd800b3a8 2011-02-23        kinaba: int main() { PSequence().run_test(-1); }
4fd800b3a8 2011-02-23        kinaba: // END CUT HERE