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

void solve(){
  int n,p;
  cin>>n>>p;
  vector<int> a(n),tmp(n);
  for (int i=0;i<n;i++){
    cin>>a[i];
    a[i]--;
    tmp[a[i]]=i;
  }
  swap(a,tmp);
  vector<vector<int>> g;
  vector<pair<int,int>> g2;
  {
    vector<bool> vis(n);
    vector<int> vec;
    auto dfs=[&](auto dfs,int v)-> void {
      if (vis[v]){
        if ((int)vec.size()==2) g2.push_back({vec[0],vec[1]});
        else g.push_back(vec);
        vec.clear();
        return;
      }
      vis[v]=true;
      vec.push_back(v);
      dfs(dfs,a[v]);
    };
    for (int i=0;i<n;i++) if (!vis[i]) dfs(dfs,i);
  }
  vector<int> ans(n,-1);
  {
    int k=g2.size();
    for (int i=0;i<k/2;i++){
      auto [a,b]=g2[i*2];
      auto [c,d]=g2[i*2+1];
      ans[a]=c;
      ans[b]=d ;
      ans[c]=b;
      ans[d]=a;
    }
    if (k%2==1){
      auto [a,b]=g2.back();
      ans[a]=b;
      ans[b]=b;
    }
  }
  {
    for (auto vec:g){
      int k=vec.size();
      if (k==1){
        ans[vec[0]]=vec[0];
        continue;
      }
      int j=-1;
      for (int i=0;i<k;i++){
        if ((i*2*p)%k==1){
          j=i;
        }
      }
      assert(j!=-1);
      for (int i=0;i<k;i++){
        int a=vec[i];
        int b=vec[(i+j)%k];
        ans[a]=b;
      }
    }
  }
  for (int i=0;i<n;i++){
    assert(ans[i]!=-1);
    cout<<ans[i]+1<<" \n"[i==n-1];
  }
  cout<<flush;
}

int main(){
  int t;
  cin>>t;
  while (t--) solve();
}