oj.uz

Submission #367782

# Submission time Handle Problem Language Result Execution time Memory
367782 2021-02-18T10:04:43 Z arnold518 Semafor (COI20_semafor) C++14
12 / 100
 450 ms 28652 KB 
#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef vector<vector<ll>> Mat;

const int MAXN = 1500;
const ll MOD = 1e9+7;

const int MM[10]={12, 8, 5, 25, 10, 19, 20, 9, 23, 27};

int M, X;
ll N, K;

Mat operator + (const Mat &p, const Mat &q)
{
	Mat ret;
	int sz=p.size();
	ret=vector<vector<ll>>(sz, vector<ll>(sz));
	for(int i=0; i<sz; i++)
	{
		for(int j=0; j<sz; j++)
		{
			ret[i][j]=0;
			for(int k=0; k<sz; k++)
			{
				ret[i][j]+=p[i][k]*q[k][j]%MOD;
				ret[i][j]%=MOD;
				//if(ret[i][j]>=MOD) ret[i][j]-=MOD;
			}
		}
	}
	return ret;
}

Mat mypow(Mat x, ll y)
{
	if(y==0)
	{
		int sz=x.size();
		Mat ret=vector<vector<ll>>(sz, vector<ll>(sz, 0));
		for(int i=0; i<sz; i++) ret[i][i]=1;
		return ret;
	}
	if(y%2) return mypow(x, y-1)+x;
	Mat ret=mypow(x, y/2);
	return ret+ret;
}

Mat A, B, C;
Mat P[MAXN+10], Q, R, S;

ll comb[MAXN+10][MAXN+10];
ll ans[100];

int main()
{
	scanf("%d%lld%lld%d", &M, &N, &K, &X);
	for(int i=1; i<=K; i++)
	{
		for(int j=0; j<=i; j++)
		{
			if(j==0 || j==i) comb[i][j]=1;
			else comb[i][j]=(comb[i-1][j]+comb[i-1][j-1])%MOD;
		}
	}

	A=vector<vector<ll>>(32, vector<ll>(32));
	for(int i=0; i<32; i++)
	{
		for(int j=0; j<32; j++)
		{
			if(__builtin_popcount(i^j)==1) A[i][j]=1;
		}
	}

	P[0]=mypow(A, 0);
	for(int i=1; i<=K; i++)
	{	
		P[i]=P[i-1]+A;
	}


	Q=vector<vector<ll>>(100, vector<ll>(100));
	for(int i=0; i<=K; i++)
	{
		for(int p=0; p<100; p++)
		{
			for(int q=0; q<100; q++)	
			{
				Q[p][q]+=P[i][MM[p/10]][MM[q/10]]*P[K-i][MM[p%10]][MM[q%10]]%MOD*comb[K][i]%MOD;
				Q[p][q]%=MOD;
			}
		}
	}

	R=vector<vector<ll>>(100, vector<ll>(100));
	for(int i=0; i<=N%K; i++)
	{
		for(int p=0; p<100; p++)
		{
			for(int q=0; q<100; q++)	
			{
				R[p][q]+=P[i][MM[p/10]][MM[q/10]]*P[N%K-i][MM[p%10]][MM[q%10]]%MOD*comb[N%K][i]%MOD;
				R[p][q]%=MOD;
			}
		}
	}

	S=mypow(Q, N/K);

	for(int i=0; i<100; i++)
	{
		printf("%lld\n", S[X][i]);
	}
	return 0;
/*
	B=mypow(A, K);
	C=mypow(A, N%K);

	P=vector<vector<ll>>(10, vector<ll>(10));
	Q=vector<vector<ll>>(10, vector<ll>(10));

	for(int i=0; i<10; i++)
	{
		for(int j=0; j<10; j++)
		{
			P[i][j]=B[MM[i]][MM[j]];
			Q[i][j]=C[MM[i]][MM[j]];
		}
	}

	R=mypow(P, N/K)+Q;

	for(int i=0; i<10; i++)
	{
		printf("%lld\n", R[X][i]);
	}
*/
}

Compilation message

semafor.cpp: In function 'int main()':
semafor.cpp:60:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   60 |  scanf("%d%lld%lld%d", &M, &N, &K, &X);
      |  ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Subtask #1
0 / 6.0

# Verdict Execution time Memory Grader output
1 Incorrect 6 ms 876 KB Output isn't correct
2 Halted 0 ms 0 KB -

Subtask #2
0 / 15.0

# Verdict Execution time Memory Grader output
1 Incorrect 6 ms 876 KB Output isn't correct
2 Halted 0 ms 0 KB -

Subtask #3
12.0 / 12.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 1004 KB Output is correct
2 Correct 34 ms 2284 KB Output is correct
3 Correct 303 ms 18156 KB Output is correct
4 Correct 339 ms 21356 KB Output is correct
5 Correct 350 ms 21484 KB Output is correct
6 Correct 386 ms 24692 KB Output is correct
7 Correct 446 ms 28652 KB Output is correct
8 Correct 450 ms 28652 KB Output is correct
9 Correct 324 ms 20460 KB Output is correct

Subtask #4
0 / 12.0

# Verdict Execution time Memory Grader output
1 Incorrect 71 ms 2080 KB Output isn't correct
2 Halted 0 ms 0 KB -

Subtask #5
0 / 15.0

# Verdict Execution time Memory Grader output
1 Incorrect 71 ms 2080 KB Output isn't correct
2 Halted 0 ms 0 KB -

Subtask #6
0 / 40.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 1004 KB Output is correct
2 Correct 34 ms 2284 KB Output is correct
3 Correct 303 ms 18156 KB Output is correct
4 Correct 339 ms 21356 KB Output is correct
5 Correct 350 ms 21484 KB Output is correct
6 Correct 386 ms 24692 KB Output is correct
7 Correct 446 ms 28652 KB Output is correct
8 Correct 450 ms 28652 KB Output is correct
9 Correct 324 ms 20460 KB Output is correct
10 Incorrect 71 ms 2080 KB Output isn't correct
11 Halted 0 ms 0 KB -