#include<bits/stdc++.h>
using namespace std;


#define MAXN 100005

long long n,k;
long long a[MAXN];
long long g[MAXN];
long long f[MAXN];

const long long mod=100003;
long long fpow(long long a,long long b,long long p){
	long long ret=1;
	for(;b;b>>=1,a*=a,a%=p){
		if(b&1){
			ret*=a; ret%=p;
		}
	}
	return ret;
}
long long fpow(long long a,long long b){
	return fpow(a,b,mod);
}
long long inv(long long x){
	return fpow(x,mod-2);
}

int main(){
	ios::sync_with_stdio(false);
	
	cin>>n>>k;
	for(long long i=1;i<=n;i++){
		cin>>a[i];
	}
	
	long long ans=0;
	for(long long i=n;i>=1;i--){
		if(a[i]){
			ans++;
			for(long long j=1;j<=sqrt(i);j++){
				if(i%j==0){
					a[j]^=1; 
					if(j!=(i/j)){
						a[i/j]^=1;
					}
				}
			}
		}
	}
	
	g[n]=1;
	for(long long i=n-1;i>=k+1;i--){
		g[i]=(n-i)*inv(i)*g[i+1]+n*inv(i);
		g[i]%=mod;
	}
	for(long long i=1;i<=k;i++){
		f[i]=i;
	}
	for(long long i=k+1;i<=n;i++){
		f[i]=f[i-1]+g[i];
	}
	for(long long i=1;i<=n;i++){
		f[ans]=(f[ans]*i)%mod;
	}
	cout<<f[ans]<<endl;
	
	return 0;
}